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PREFACE 

This  bulletin  contains  a  collection  of  articles  ranging  from  analysis  of  observa¬ 
tional  data  on  the  marginal  ice  zone  to  theoretical  papers  relevant  to  the  physics  of 
air/ice/ocean  dynamics  in  the  MIZ.  As  such  it  represents  a  more  general  format  for 
articles  describing  the  marginal  ice  zone  than  previous  bulletins.  We  will  continue 
this  more  general  format  in  most  future  issues,  with  occasional  special  issues  inter¬ 
spersed.  Authors  are  encouraged  to  submit  articles  on  any  aspect  of  marginal  ice 
zone  studies  in  both  the  northern  and  southern  hemispheres. 
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A  Note  on  Estimating  Melt  Rate  in  the  MIZ 


MILES  G.  McPHEE 

McPhee  Research  Company,  Route  8,  Box  304H,  Yakima,  Washington  98908 


How  heat  and  salt  are  exchanged  at  the  ice/ 
ocean  interface  is  a  central  issue  in  understanding 
the  thermodynamics  and  the  dynamics  of  the  mar¬ 
ginal  ice  zone.  The  intent  of  this  note  is  to  assem¬ 
ble  some  ideas  about  boundary  conditions  at  the 
horizontal  ice/ ocean  interface,  to  investigate  what 
they  imply  for  the  oceanic  boundary  layer,  and  to 
present  an  approximate  method  for  estimating 
melt  rates  of  sea  ice  when  the  predominant  heat 
transfer  is  from  the  ocean.  The  note  is  purposely 
informal,  and  if  it  generates  further  discussion,  1 
will  consider  it  worthwhile. 

In  the  1983  MIZ  issue  of  the  Journal  of  Geo¬ 
physical  Research,  Ed  Josberger  and  I  each  pub¬ 
lished  articles  (Josberger,  1983;  McPhee,  1983)  on 
how  sea  ice  melts  in  above-freezing  water.  Consid¬ 
ering  that  we  used  rather  similar  methods,  it  was 
somewhat  surprising  that  we  came  to  diametrically 
opposed  conclusions  regarding  the  importance  of 
buoyancy  flux  in  the  dynamics  of  the  planetary 
boundary  layer  (PBL).  I  wrote  a  comment  on  Jos- 
berger’s  article  (McPhee,  1984)  to  the  effect  that 
he  had  incorrectly  calculated  the  Obukhov  length, 
and  he  replied  (Josberger,  1984)  that  I  failed  to 
treat  conditions  at  the  ice/ ocean  interface  proper¬ 
ly. 

To  some  extent,  Ed  was  right.  The  fact  is,  I  was 
not  sure  what  to  use  for  a  boundary  condition  at 
the  ice/ocean  interface  (I  was  unconvinced  that  sa¬ 
line  ice  would  melt  under  equilibrium  conditions 
in  the  complex  environment  that  makes  up  the  in¬ 
terface),  so  I  avoided  the  question  by  formulating 
the  problem  in  terms  of  change  in  mean  tempera¬ 
ture  across  the  PBL,  and  by  relating  the  surface 
salinity  flux  to  the  product  of  ice  growth  rate  and 
the  difference  between  mixed-layer  and  ice  salini¬ 
ties.  My  approach  amounts  to  considering  some 
level  a  little  removed  from  the  actual  interface, 
then  assuming  that  the  momentum  flux  (stress) 
and  heat  and  salinity  fluxes  are  the  same  there  as 
they  are  at  the  interface,  a  common  and  useful  ap¬ 
proximation  in  boundary-layer  studies.  If  one  as¬ 
sumes  that  the  average  heat  and  salt  fluxes  are 


proportional  to  the  overall  ice  growth  rate,  then 
the  problem  can  be  solved  without  direct  reference 
to  So  and  To,  the  interfacial  salinity  and  tempera¬ 
ture,  even  though  it  is  somewhat  unsatisfying  to 
deal  with  changes  across  the  PBL  rather  than  ab¬ 
solute  values. 

As  Ed  and  I  were  airing  our  differences,  George 
Mellor,  who  was  visiting  the  Naval  Postgraduate 
School  at  the  same  time  I  was,  took  interest  in  the 
ice/ocean  PBL  problem.  George  used  the  freezing¬ 
line  T0/ S0  relationship  to  put  the  boundary  condi¬ 
tions  for  his  numerical  PBL  model  capped  with  ice 
on  a  rigorous  footing.  We  have  since  applied  the 
Mellor-Yamada  ‘level  2.5’  time-dependent  model 
to  the  regime,  first  results  of  which  will  soon  be 
out  in  a  manuscript  by  Mellor,  McPhee,  and 
Steele.  This  note  discusses  those  boundary  condi¬ 
tions  as  George  formulated  them,  and  their  impli¬ 
cation  for  simpler,  steady-state  PBL  models. 

Consider  a  thermodynamic  control  volume  of 
infinitesimal  extent  surrounding  the  ice/ocean  in¬ 
terface  and  moving  with  it,  as  diagrammed  in  Fig¬ 
ure  1 .  We  assume  that  neither  mass  nor  energy  ac¬ 
cumulates  within  the  volume  and  describe  the  rate 
of  freezing  by  the  ice  growth  rate,  d.  Mass  conser¬ 
vation  and  isostasy  require  the  vertical  velocity  of 
fluid  into  the  control  volume  to  be  w  =  -(e/po)  d. 
(Dynamically,  this  ‘injection  velocity’  is  insignifi¬ 
cant  relative  to  other  fluid  velocities  and  is  ignored 
in  the  equations  of  motion.)  The  change  in  enthal¬ 
py  is  dominated  by  the  phase  change  (about  60 
times  the  sensible  heat  change),  so  the  first-law 
balance  is  approximately: 
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where  ks  ( dT\/dz)\»  is  heat  conduction  through  the 
ice  and  p0  c  <w'T'> 0  is  heat  flux  from  the 
ocean.  Latent  heat,  Ls,  and  thermal  conductivity, 
ks,  are  functions  of  ice  salinity  and  can  be  ap¬ 
preciably  different  from  their  fresh-ice  counter¬ 
parts.  The  ‘kinematic’  heat  flux  is: 
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Fit;.  1.  Schematic  of  thermodynamic  control  volume  surrounding  ice/water  interface. 
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typical  value  of  about  60  °C. 

Salt  conservation  is  treated  similarly.  As  the 
control  volume  moves,  flux  and  advection  at  the 
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lower  surface  are  balanced  by  advection  only  at 
the  upper  surface,  since  we  ignore  salt  diffusion  in 
the  ice: 

Qo  <w'S'>o  +  Qo  W  Sj  =  Qo  W  So* 

The  salinity  flux  is  thus 

<w'S’>o  =  w(So-Sj)  (2) 

When  ice  melts,  the  salinity  flux  is  positive, 
meaning  that  salt  (and  mass)  are  transported  up¬ 
ward,  which  is  to  say  that  relatively  fresh  melt 
water  is  mixed  down  by  turbulent  (and  molecular) 
processes.  At  temperatures  near  freezing,  the  ther¬ 
mal  expansion  coefficient  of  sea  water  is  small, 
and  buoyancy  is  governed  mainly  by  the  salt  flux. 

What  do  the  surface  flux  conditions  imply 
about  the  boundary  layer?  First,  suppose  we  know 
the  mean  salinity  and  temperature  at  some  level 
near  the  interface  in  the  turbulent  flow.  By  choos¬ 
ing  the  level  close  enough,  the  ‘wall-layer’  condi¬ 
tions  of  constant  flux  and  small  buoyancy  effects 
can  be  approximated.  We  further  assume  that  tur¬ 
bulent  diffusivities  for  momentum,  heat,  and  salt 
are  the  same,  so  we  have 


dT 

as 

dz  ’ 

dz 

where  k  is  von  Karman’s  constant,  equal  to  0.4. 

According  to  Gill  (1982),  the  freezing  tempera¬ 
ture  of  sea  water  is 

Tf(S,p)  =  -0.0575S  +  1.710523  x  10” Sv’ 

-2.1 54996  x  10”  S2  -7.53  x  10’p 

where  S  is  in  ppt  and  p  in  bars.  This  is  often  linear¬ 
ized  to 


77  =  -mS 


where  m  =  0.06  °C/(%,) 


If  interfacial  stress  (q0«,2)  is  known,  we  now 
have  a  third  equation  for  the  unknowns  T0,  So, 
and  w  that  is  readily  solved.  We  used  this  ap¬ 
proach  to  tie  the  interface  quantities  to  flow  pro¬ 
perties  of  the  first  grid  point  in  the  Mellor- 
Yamada  numerical  model. 

Now  suppose  we  have  neither  measurements 
nor  a  numerical  model  to  specify  conditions  in  the 
near  surface  layer,  but  that  we  can  estimate  the 
‘far-field’  temperature,  and  salinity  and  the  stress 
acting  at  the  ice/water  interface.  Is  there  a  simple 
way  of  estimating,  say,  melt  rate?  Perhaps.  Two 
factors  make  this  problem  less  tractable  than  the 
surface  layer  approach  described  above:  first,  tur¬ 
bulent  flux  is  not  constant  across  the  entire  PBL, 
but  instead  decreases  from  its  surface  value  to  zero 


and,  second,  the  eddy  viscosity  loses  its  simple  lin¬ 
ear  ^-dependence.  Actually,  the  similarity  model  I 
used  in  McPhee  (1983)  is  quite  simple,  but  requires 
some  familiarity  with  my  notation  and  peculiar 
way  of  viewing  the  PBL.  In  a  sense,  what  follows 
is  a  further  simplification  of  that  theory,  which 
may  help  explain  its  function.  Josberger  used  an 
analogous  technique  but,  in  addition  to  our  differ¬ 
ence  about  what  constitutes  buoyancy  flux,  made 
two  assumptions  with  which  I  disagree,  as  ex¬ 
plained  below. 

It  helps  to  review  some  rudimentary  concepts 
about  steady,  horizontally  homogeneous  plane¬ 
tary  boundary  layers.  First,  almost  by  definition, 
the  PBL  extends  through  the  region  where  turbu¬ 
lent  flux  of  momentum  and  scalar  properties  oc¬ 
curs.  By  implication,  these  fluxes  must  fall  off  to 
zero  somewhere  near  the  outer  limits  of  the  PBL. 
This  is  the  source  of  one  of  my  disagreements  with 
Josberger — he  assumed  constant  flux  for  heat  and 
salt  (but  not  momentum)  through  the  entire  PBL. 
Second,  there  seem  to  be  three  length  scales  that 
affect  the  large-scale  behavior  of  turbulence  in  the 
PBL.  They  are:  1)  a  scale  proportional  to  u,/f, 
which  derives  from  a  conditional  stability  related 
to  rotation  of  the  reference  frame;  2)  the  distance 
from  the  fluid/solid  interface;  and  3)  the  Obu¬ 
khov  length,  which  is  essentially  the  vertical  dis¬ 
tance  a  fluid  parcel  moves  in  a  density-stratified 
flow  before  its  initial  kinetic  energy  is  converted  to 
gravitational  potential  energy.  In  general,  turbu¬ 
lence  at  a  given  level  in  the  PBL  responds  to 
whichever  length  scale  is  shortest  at  that  point.  Us¬ 
ing  this  paradigm,  my  studies  have  shown  that  the 
length  scale  of  the  energy-containing  eddies  in  a 

Heat  flux,  <w'T'  > 

0  < w  T' >0 


neutrally  stratified  flow  will  increase  until  it  reach¬ 
es  roughly  (l/20)*w.//.  Beyond  that  level,  the  tur¬ 
bulence  is  ‘free,’  meaning  the  eddies  are  no  longer 
affected  by  the  wall  itself.  The  neutral  boundary 
layer  seems  to  extend  to  about  8  or  10  times  the 
scale  of  the  largest  eddies,  say  0.4  to  0.5  u,/f. 

If  the  stratifying  tendency  is  strong,  as  it  is  when 
ice  melts  rapidly,  the  Obukhov  length  is  short, 
limiting  the  maximum  length  scale  of  the  eddies 
relative  to  the  neutral  case.  The  basic  idea  behind 
similarity  in  fluid  dynamical  problems  is  to  find 
scales  that  collapse  a  whole  class  of  physical  real¬ 
izations  to  one  equation,  the  solution  of  which 
then  has  a  general  nature.  My  similaiity  theory 
(McPhee,  1981)  for  the  outer  layer  combines  the 
planetary  scale  and  the  Obukhov  length  in  a  scale 
for  the  effective  vertical  influence  of  the  largest 
eddies.  When  this  scale  is  specified,  other  proper¬ 
ties  (e.g.,  the  depth  of  the  PBL,  which  is  constant 
in  the  similarity  coordinates)  fall  out  automatical¬ 
ly.  Increased  surface  buoyancy  decreases  the 
length  scale,  limits  the  size  of  the  boundary  layer, 
and  traps  the  effect  of  surface  fluxes  closer  to  the 
surface.  There  is  a  significant  feedback  here:  as 
turbulent  heat  flux  increases,  so  does  melting, 
which  inhibits  the  efficiency  of  turbulent  transfer. 
For  this  reason,  a  bulk  Stanton-number  approach 
for  calculating  heat  flux  should  be  examined  criti¬ 
cally,  if  used  in  a  regime  with  turbulence  levels  dif¬ 
ferent  from  that  used  to  derive  the  Stanton  num¬ 
ber.  An  interesting  consequence  of  the  same  effect 
is  that  by  confining  momentum  flux  to  a  shallower 
region,  increased  buoyancy  causes  faster  near- 
surface  currents  (and  ice  velocity)  for  the  same  in¬ 
terfacial  stress.  Buoyancy-caused  change  in  the 


Fig.  2.  Flux  profiles  in  boundary  layer  under  melting  ice. 


PBL  depth  is  a  second  source  of  my  disagreement 
with  Josberger:  he  sets  the  PBL  depth  to  the  limit  I 
suggested  for  the  neutral  case  (0.5  u,/f)  and  leaves 
it  there  regardless  of  surface  buoyancy  flux. 

Consider  a  melting  regime  with  no  heat  conduc¬ 
tion  through  the  ice,  as  sketched  in  Figure  2.  The 
dashed  lines  are  meant  to  indicate  that  there  is 
nothing  a  priori  known  about  the  shape  of  the 
profiles,  just  that  they  fall  from  surface  values  to 
near  zero  at  the  far  extent  of  the  PBL.  We  can 
simplify  matters  greatly  by  assuming  that  a)  the 
form  of  the  scalar  flux  profiles  is  similar,  and  that 
b)  the  eddy  diffusivities  are  the  same  for  salinity 
and  heat,  as  are  c)  the  surface  roughnesses.  Each 
of  these  assumptions  occupies  long  (sometimes  te¬ 
dious)  hours  of  argument  among  turbulence  stu¬ 
dents,  but  for  the  present  purposes  it  is  safe  to  say 
that  none  is  in  gross  error.  With  that,  we  can  im¬ 
mediately  write: 

~  MSo-Sj)] 

where  F  is  a  function  that  describes  the  shape  of 
both  flux  profiles  and  K  is  eddy  diffusivity. 

Integrating  from  the  base  of  the  PBL,  we  hnve 

UTv-TJASo-SJ]  =  [Qlw,  w(S0-Sj)]  • 


act  specification  of  conditions  at  the  interface  is 
not  crucial  to  specifying  heat  and  salt  flux.  Note 
that  this  result  follows  from  assumptions  a 
through  c  listed  above,  and  not  from  any  particu¬ 
lar  form  specified  for  the  flux  profile  shape  or  ed¬ 
dy  diffusivity. 

For  a  useful  estimate  of  melt  rate,  we  have  to 
address  the  integral  in  equation  5;  specifically,  the 
shape  of  the  flux  profile  and  the  eddy  viscosity.  It 
turns  out  that  assuming  a  linear  fall  off  to  scalar 
fluxes  is  not  a  bad  description  for  a  number  of 
modeling  and  experimental  studies.  Much  more 
important  than  the  exact  form  of  the  curve  is  the 
depth  to  which  turbulence  penetrates.  Without  go¬ 
ing  into  the  details  of  its  derivation,  1  worked  out 
a  theory  several  years  ago  (McPhee,  1981),  using 
the  three  scales  described  at  the  beginning  of  this 
section,  that  specifies  a  scale  for  the  free  turbu¬ 
lence.  As  a  first-order  approximation  to  that 
theory,  eddy  viscosity  may  be  written: 

K  —  ~**o  >  z  >  Zsi 

K  =  k$HuW,/f\  Zsi  >  z  >  Zbi 
where 
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and  with  the  freezing  line  approximation, 
-m  Sa,  and  some  rearranging. 


T0  = 


where  £N  is  a  constant,  about  0.05,  and  >7,  is  a  sta¬ 
bility  parameter  described  below.  Assuming  heat 
flux  is  linear  from  its  surface  value  to  zero  at  Zbi. 
the  equation  may  be  integrated  in  two  parts, 
through  the  outer,  free-turbulence  layer,  and 
through  a  surface  layer: 
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This  can  be  written  as  a  quadratic  in  SQ  and 
solved,  but  casting  it  in  this  form  shows  that  the 
ratio  is  dominated  by  (which  is  about  60°C),  so 
that  for  realistic  MIZ  values  of  S  and  T  it  is 
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always  within  a  few  percent  of  unity.  Similarly,  T0 
is  within  a  few  tenths  of  a  degree  of  the  freezing 
point  of  the  far  field,  -mS^.  For  practical  pur¬ 
poses,  the  difference  between  S0  and  is  over¬ 
shadowed  by  the  uncertainty  in  Sj  and  Q\,  and  ex- 
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A  straightforward  scale  analysis  shows  that  the  in¬ 
tegral  is  dominated  by  two  terms,  and  the  equa¬ 
tion  may  be  written 


A  = 


w.(VV)  1 


A, 


A 

1 


1  /«.< N  .V  ' 

wQl  k  10  I/ZOT  V'r  2Sn„. 

(7) 


A  =  A,  +  A;  is  the  inverse  of  a  heat  transfer  coeffi¬ 
cient  and  comprises  two  parts:  A,  is  the  surface 
layer  contribution  to  the  boundary-layer  integral, 
and  A;  is  the  integral  through  the  outer  (Ekman) 
layer.  Note  that  each  depends  on  ij.,  but  that  the 
surface  layer  dependence  is  logarithmic,  thus  weak¬ 
er.  A,  is  also  weakly  dependent  on  //,, 

The  buoyancy  scale  factor  »/,,  derived  in 
McPhee  (1981),  is  given  by 


/  £nw>  y ; 
7*  =  l,+  KJU 


where  L  =  fe>0wJ/(gAr  <  w'  q  '  >„)  is  the  Obukhov 
length  and  /?c  is  the  critical  flux  Richardson  num¬ 
ber,  equal  to  0.2.  r;„  is  1  for  neutral  stability  (no 
buoyancy)  and  less  than  1  for  all  stably  stratified 
flows.  The  buoyancy  flux  is  approximately 


g<w'c’>  o  _  ,c, 

-  S  gos  <  w  S  >„ 


Q  o 


where  the  thermal  expansion  term  is  neglected  at 
low  temperatures.  From  equation  2  and  with  the 
following:  «s  =  8  x  10~4  (%0)'',  g  =  9.8  m  s  ’,  and 
/  =  1.4  x  10  4  s'1,  the  expression  for  //,  is  approx¬ 
imately 


a  [l+a. 


(8) 


3  x  10“  m.  My  intuition  is  that  this  underestimates 
Cot  *n  complex  under-ice  environment,  and 
that  the  value  chosen  above,  namely  10“  m,  repre¬ 
sents  a  compromise.  Nevertheless,  1  now  think 
that  the  value  used  in  McPhee  (1983),  0.05  m,  is 
too  large  for  cOT,  though  not  necessarily  for  z0. 
For  reference,  making  c0(  30  times  smaller  in  the 
example  above  adds  about  8.5  to  A,,  the  surface- 
layer  contribution  to  the  boundary-layer  integral, 
which  is  still  considerably  smaller  than  the  buoy¬ 
ancy  contribution. 

Given  surface  stress  (y0nj),  here  is  a  method  for 
estimating  melt  rate.  First,  estimate  the  change  in 
temperature  (T0-Toc)  and  salinity  difference 
(S0-Si),  either  directly  or  from  Tx  and  S*  using 
equation  6.  Initially,  assume  »/,  =  1  (neutral  con¬ 
ditions),  and  use  equation  7  to  calculate  an  upper 
bound  for  w,  which  is  used  in  equation  8  to  obtain 
a  refined  estimate  of  >?..  The  procedure  is  repeated 
until  the  value  for  w  converges  to  within  some  tol¬ 
erance;  in  practice,  this  convergence  occurs  rap¬ 
idly.  Melt  rate  is  (e0/ej)>v.  The  following  table 
shows  the  calculation  for  u,  =  0.01  m  s  ',  (Tq-Too) 
=  2°C,  and  (S0-Sj)  =  25%„: 


Iteration 

V. 

A 

w(x  10’ m  s  ') 

1 

1 

30.5 

1.09 

2 

0.25 

53.6 

0.62 

3 

0.32 

46.2 

0.72 

4 

0.30 

47.8 

0.70 

where  a,  =  5.6  m  s'1  (%0)  '. 

The  effect  of  n.  on  heat  transfer  is  demonstrat¬ 
ed  as  follows:  let  it.  be  0.01  m  s  ',  and  assume  first 
that  there  is  no  surface  buoyancy  flux,  so  that 
=  I.  With  cqx  =  10“  m  (see  discussion  below).  A, 
is  20.5  and  A-  is  10.  If  we  now  specify  a  melt  rate 
of  50  cm  da  '  (w  =  5  x  10'*  m  s  ')  and  a  salinity 
difference  (S0-Sj)  of  25,  then  according  to  equa¬ 
tion  8,  »).  =  0.35.  The  change  in  A,  (equal  to  2  In 
[>/.]/0.4)  is  about  -5.2,  while  A2  increases  by  18.6. 
There  is  a  significant  net  increase  in  A  from  the  in¬ 
tegral  across  the  outer  layer,  and  a  corresponding 
30%  decrease  in  the  heat  transfer  coefficient. 

The  temperature  surface  roughness,  Cot-  war' 
rants  comment.  In  McPhee  (1983),  I  equated  it  to 
the  momentum  surface  roughness,  z0,  which  is  rel¬ 
atively  large  under  sea  ice.  George  Mellor  has 
pointed  out  that  the  fundamental  processes  that 
deterine  z0  include  form  drag  forces  across  rough¬ 
ness  elements,  for  which  there  is  no  analog  in  the 
heat  transfer  process.  He  suggests  using  cor  = 
v^/ku,,  where  eT  is  the  molecular  thermal  diffus- 
ivity.  For  it.  =  0.01  m  s  ',  this  puts  cOT  at  about 


Melt  rate  is  about  67  cm/day. 

This  relatively  simple  calculation  can  furnish 
additional  information,  such  as  the  depth  to  which 
the  “melt-water  lens”  penetrates,  which  is 
ky, u,/f,  or  approximately  9  m.  The  technique 
may  also  be  used  to  investigate  ice  drift  and  near¬ 
surface  currents,  but  it  is  then  probably  easier  to 
refer  to  the  original  similarity  theory.  A  reader 
who  has  followed  the  arguments  here  will  see  that 
they  are  applied  in  the  same  way  there,  except  that 
the  scale  factor  is  derived  from  more  basic  as¬ 
sumptions. 

Note  that  melt  rate  calculated  this  way  will  per¬ 
sist  only  as  long  as  heat  is  available.  With  “off¬ 
ice”  advection,  the  leading  ice  edge  will  encounter 
warm  water  continuously,  but  trailing  floes  inter¬ 
act  with  water  that  has  been  cooled  and  undergo 
less  melting.  To  treat  the  problem  in  any  detail  re¬ 
quires  consideration  of  both  time  dependence  and 
horizontal  gradients,  which  add  considerable  com¬ 
plexity. 
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Kinematics  of  Marginal  Ice:  MIZEX  83 
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The  kinematics  of  the  marginal  ice  zone,  an  im¬ 
portant  element  in  ice  dynamics  and  hence  in  the 
air/ice/water  system  of  the  area,  were  studied  in 
the  Fram  Strait,  east  of  Greenland,  as  a  part  of 
MIZEX  83  (inset.  Figure  1).  The  primary  purpose 
of  the  work  was  as  a  pilot  study  for  the  main  pro¬ 
gram  in  1984  (Johannessen  et  al.,  1983).  Some  in¬ 
teresting  phenomena  detected  directly  from  the  pi¬ 
lot  study  are  presented  in  this  article. 


Ice  margin 

At  the  ice  margin  the  frozen  sea  meets  the  open 
sea.  The  general  characteristics  of  the  transition 
are  schematically  described  first.  Two  facts  are  es¬ 
sential  to  establish  an  ice  margin.  First,  a  quantity 
(for  example,  fish  population,  air  temperature, 
floe  diameter,  to  mention  a  few)  that  has  a  con¬ 
stant  value,  Q\,  (or  changes  its  value  linearly  with 
distance)  in  the  region  far  into  the  frozen  sea.  (In 
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fig.  I.  Position  of  markers  /,  2,  3.  and  4  at  1030  and  1800  GMT,  1 1  July  1983.  I  he  path  line  <>!  marker  4  and  the 
tee  edge  loealion  a  I  1200  GMT  are  also  shown. 
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Fig.  2.  Schematic  transition  from  frozen  sea  to 
open  sea. 


the  latter  case,  its  gradient  is  used.)  Such  a  place  is 
called  the  “interior  zone”— to  be  precise,  an  inte¬ 
rior  zone  for  the  quantity — and  it  is  almost  always 
present,  since  the  extent  of  the  area  is  not  restrict¬ 
ed  and  in  the  extreme  can  be  reduced  to  a  line.  Sec¬ 
ond,  the  value  of  the  quantity  must  be  prescribed 
at  the  ice/water  boundary  (i.e.,  the  ice  edge),  Qw 
(Figure  2). 

There  may  be  a  “boundary”  zone  that  connects 
the  interior  zone  to  the  ice  edge.  How  the  quantity 
changes  beyond  the  ice  edge  into  the  water  is  not 
of  concern  here;  indeed,  it  need  not  even  be  de¬ 
fined  over  the  open  sea.  The  boundary  zone  thus 
has  a  secondary  meaning  to  construct  the  ice  mar¬ 
gin,  in  the  sense  that  a  transition  can  also  exist 
without  the  presence  of  a  boundary  zone  (e.g., 
open  water  washing  a  solid  fast  ice  field).  Never¬ 
theless,  the  boundary  zone  is  most  interesting.  All 
behaviours  that  are  associated  with  the  ice/water 
transition  can  be  expected  to  be  present  in  the 
boundary  zone  in  extreme  form.  The  final  purpose 
of  this  study  is  to  obtain  the  actual  shape  of  the 
profile  over  the  boundary  zone  for  the  quantity  in 
question. 

More  than  one  boundary  zone  can  exist,  possi¬ 
bly  with  intermediate  interior  zones  in  between.  In 
that  case  the  schema  has  only  to  be  repeated;  the 
treatment  may  become  more  complicated,  but 
nothing  essentially  new  is  added.  When  the  quan¬ 
tity  is  a  tensor,  the  schema  is  applied  on  its  com¬ 
ponents  or  on  its  invariants.  The  profile  over  the 
boundary  zone  can  change  its  shape  drastically 
with  time,  for  each  quantity  in  question  and  along 
the  ice  edge,  while  the  basic  schema  is  applicable 
throughout. 


Measurement 

Four  markers,  1  through  4,  were  set  on  ice  floes, 
and  their  relative  positions  were  noted  every  30 
minutes  during  the  time  1030  to  1800GMT,  1 1  July 
1983,  by  ordinary  survey  instruments  (Theodolite 
Wild  T2,  Distomat  Wild  DI4L).  The  absolute  po¬ 
sition  was  also  determined,  but  with  much  less  ac¬ 
curacy,  through  the  ship  position  and  the  com¬ 
pass.  The  positions  of  the  markers  at  the  first  and 
last  measurements  are  shown  in  Figure  1  with  the 
path  line  of  marker  4.  The  ice  edge  at  1200  GMT  is 
also  plotted  in  the  figure  (personal  communica¬ 
tion,  Strubing).  Since  the  ice  edge  also  moved, 
none  of  the  markers  crossed  the  ice  edge  during 
the  observation. 

The  accuracy  of  individual  distance  and  angle 
measurements  was  estimated  to  be  within  2  cm 
and  0.002  400°  respectively.  Due  to  the  rapid 
movement  (relative  to  the  time  required  to  com¬ 
plete  a  reading)  and  insufficiently  accurate  time 
determination,  the  positions  of  the  markers  are  es¬ 
timated  to  be  correct  within  several  meters  relative 
to  each  other.  The  position  of  the  markers  relative 
to  the  earth  is  assumed  to  be  reliable  within  hun¬ 
dreds  of  meters. 

The  time  interval  of  the  measurement  was  some¬ 
what  irregular,  but  the  nominal  time  of  30  minutes 
is  used  throughout  this  article,  except  in  the  cal¬ 
culation  of  velocity,  where  the  actual  time  interval 
was  used. 

The  wind  record  for  the  day  on  the  Polarstern  is 
given  below: 


0000  GMT 

280° 

from  north  at  15  knots 

0300 

290° 

at  13 

0600 

280° 

at  15 

0900 

280° 

at  12 

1200 

240° 

at  12 

1500 

250° 

at  12 

1800 

230° 

at  07 

2100 

210° 

at  12 

The  wind  was  constant  during  the  measurement 
in  terms  of  both  direction  and  speed.  Immediately 
before  measurement  started,  the  wind  direction 
rotated  counterclockwise  by  40°.  No  information 
is  available  on  the  sea  current  behaviour  for  the 
day. 

Analyses 

Preliminary  analysis  showed  that  the  classical  or 
engineering  treatment  cannot  be  used  because  of 
the  large  deformation.  More  complicated  methods 
of  analysis  (e.g.,  Truesdell  and  Toupin,  1960) 
have  to  be  employed. 

Frame.  A  moving  two-dimensional  cartesian 
frame  was  used  to  describe  the  present  configura- 
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Fig.  3.  Path  lines  of  the  markers  in  a  moving 
frame  (*).  The  movement  started  at  the  numbered 
end.  Centers  of  four  triangles  (*),  N,  E,  S,  and  W,  are 
also  shown.  (*  see  definitions  in  text.) 

tion.  The  movement  of  the  frame  was  determined 
by  two  conditions:  the  origins  of  the  frame’s  coor¬ 
dinates  were  always  at  the  center  of  the  four 
markers,  and  the  sum  of  the  directions  from  the 
origin  to  each  marker  (angles  from  an  arbitrary 
but  fixed  direction)  stayed  constant.  Such  separa¬ 
tion  of  rigid  translation  and  rotation  from  the  de¬ 
formation  is  generally  not  essential  but  it  was  con¬ 
venient  for  this  particular  study,  as  the  positions 
had  been  measured  with  differing  accuracy.  The 
positions  of  the  markers  in  this  frame  are  shown  in 
Figure  3. 

Strain.  The  shape  at  1030  GMT  was  taken  to  be 
the  reference  configuration.  A  set  of  three  mark¬ 
ers  formed  a  triangle.  There  were  four  triangles. 
Under  the  linear  assumption  within  a  triangle,  the 
deformation  gradient  was  calculated  and  that 
value  was  attached  to  the  center  of  the  triangle 
(the  centers  of  triangles  N,  E,  S,  and  W  at  1030  are 
shown  in  Figure  3).  This  procedure  is  not  com¬ 
pletely  legitimate;  two  triangles  had  a  common 
side  or  area,  where  two  different  linear  fields  were 


assumed.  Caution  is  therefore  necessary  in  the  dis¬ 
cussion.  Green’s  and  Cauchy’s  deformation  ten¬ 
sors,  rotation  tensors,  principal  axes,  and  the  val¬ 
ues  for  them  (where  applicable)  and  various  invar¬ 
iants  of  them  were  calculated.  Some  of  them  will 
be  shown  during  the  discussion. 

Strain  rate.  In  addition,  chronological  linearity 
was  assumed  during  the  preceding  30  minutes  to 
compute  the  velocity.  That  value  was  again  at¬ 
tached  to  the  ceners  of  the  triangles  and  to  the 
present  time  (neither  30  nor  15  minutes  ago). 
Various  tensors  and  invariants  were  computed  in  a 
manner  similar  to  the  strain  computation. 

Discussion 

Rigid  movement.  The  markers  moved  approxi¬ 
mately  in  the  direction  of  the  wind  with  an  average 
speed  of  700  mh'1  until  1330  (Figure  1).  Then  they 
suddenly  changed  direction  about  60°  towards  the 
right,  and  proceeded  with  a  greater  speed,  1200  m 
h*’  average,  turning  gently  more  to  the  right.  The 
wind  did  not  change  direction  (nor  speed)  through¬ 
out  the  measurement.  The  phenomenon  indicates 
one  or  both  of  two  possibilities:  the  sea  current 
changed  velocity  at  1330,  or  the  markers  were 
transferred  through  their  own  movement  relative 
to  the  zone  boundary  from  one  zone  to  another 
where  different  mechanisms  governed.  There  is  no 
means  of  checking  the  first  possibility. 

Rigid  rotation  is  plotted  in  Figure  4.  The  field 
rotated  at  first  counterclockwise  then  clockwise, 
both  almost  linearly,  with  an  average  angular 
speed  of  8°  h*1.  The  two  possibilities  given  above 
can  also  be  applied  to  explain  the  sudden  change 
in  the  rotating  direction.  The  turning  time  was, 
however,  1.5  h  earlier  in  this  case.  This  favors  the 
second  possibility:  The  geometrical  distribution  of 
zones  need  not  be  the  same  for  the  translation  and 
the  rotation.  Even  if  this  were  true,  the  sea  current 
is  not  excluded  from  the  discussion.  The  current 
may  be  an  important  element  that  specifies  the 
characteristics  of  different  zones.  The  rotation  of 
the  complete  field  to  such  an  extent  without  any 
accompanying  rotation  in  wind  is  in  itself  a  very 
interesting  phenomenon,  and  is  probably  related 
to  the  marginal  eddies. 

The  rotation  of  individual  floes  differs  from 
that  of  the  entire  field.  For  comparison,  the  ob¬ 
served  rotation  of  a  randomly  chosen  floe  on 
which  marker  4  was  located  is  plotted  in  Figure  4. 
The  curve  is  more  irregular  and  changes  value 
more  violently,  though  somewhat  systematically. 
It  is  an  indicator  of  the  local  shearing  in  the  neigh¬ 
bourhood  of  the  floe. 

The  rotation  tensor  is  an  additional  rotation  of 


Fig.  4.  Rotation  of  the  field,  positive  counterclockwise.  R— rigid  rota¬ 
tion  of  the  observation  area.  4— rotation  of  the  floe  on  which  marker  4  was  lo¬ 
cated.  W~ rotation  tensor  of  triangle  W. 


each  individual  triangle.  The  value  was  much  less 
than  the  rotations  mentioned  above,  and  nothing 
particularly  interesting  was  found.  Triangle  W  is 
shown  in  Figure  4  in  terms  of  angle,  to  indicate  the 
magnitude. 

Displacement.  Two  components  of  the  defor¬ 
mation  gradient  are  plotted  in  Figure  5.  The  area 
was  inactive  until  1200.  This  again  suggests  the 
transfer  from  one  zone  to  another  at  this  time  (or 
a  change  in  the  sea  current).  Assuming  that  the  in¬ 
active  zone  is  the  interior  zone  and  the  active  one 
is  the  boundary  zone,  the  extension  of  the  latter 
can  be  estimated  to  be  a  couple  of  kilometers  from 
the  ice  edge  (see  Figure  1).  Remembering  that  the 
ice  edge  ran  parallel  to  the  F-axis,  the  eastward 
and  southward  components  of  the  displacement 
increased  as  the  location  moved  to  the  east,  closer 
to  the  ice  edge,  after  the  observation  area  entered 
into  the  boundary  zone.  Now  the  geometrical  ar¬ 
rangement  of  the  centers  of  three  of  the  triangles 
was  in  the  order  W,  N,  E  from  west  to  the  ice  edge 
(triangle  S  is  excluded  here  for  the  reason  given  be¬ 
low).  Figure  5a  shows  the  increment  of  the  east¬ 
ward  component  decreased  at  first  from  west  to 
east  and  then,  after  1430,  increased  toward  the  ice 


edge.  This  apparent  concave-convex  transition  is 
not  essential,  but  is  associated  with  the  rotation. 

Figure  5b  shows  that  the  increase  of  the  south¬ 
ward  component  of  the  displacement  was  greater 
toward  the  ice  edge  all  the  time  after  1200.  The 
magnitude  of  the  southward  component,  hence 
also  the  difference  in  magnitude  among  the  differ¬ 
ent  triangles,  was  larger  than  that  of  the  eastward 
component. 

Two  other  components  of  the  deformation  gra¬ 
dient  are  plotted  in  Figure  6  to  show  the  change  of 
the  deformation  along  the  ice  edge.  The  behaviour 
of  triangle  S  was  quite  different  from  that  of  the 
other  three.  The  ice  motion  was  by  no  means  con¬ 
stant  parallel  to  the  ice  edge,  i.e.  one-dimensional 
treatment  and  the  concept  of  a  “band”  have  to  be 
abandoned.  Figures  6a  and  b  are  similar  to  Figures 
5a  and  b  respectively,  if  the  direction  east  (A)  is  re¬ 
placed  by  south  (-  F).  Similar  discussion  is  applic¬ 
able.  The  movement  of  the  ice  margin  is  more  ob¬ 
vious  in  the  direction  of  its  normal  for  an  “observ¬ 
er,”  as  it  is  visible  for  him,  but  the  movement  in 
the  parallel  direction  is  often  more  important. 
This  general  statement,  accepted  for  a  fixed  bound¬ 
ary  (wall)  without  question,  is  logically  supported 
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Fig.  5.  X  components  of  deformation  gradient.  N,  E,  S,  and  W  indi¬ 
cate  the  names  of  the  triangles.  X:  reference  coordinate,  x,  y:  present  coor¬ 
dinates. 


Fig.  6.  Y  components  of  deformation  gradient.  E,  S,  and  W  indicate 
the  names  of  the  triangles.  Y.-  reference  coordinate,  x,  y:  present  coordinates. 
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Fig.  7.  Direction  of  the  major  principal  axis  of  deformation,  in  de¬ 
grees  from  north.  Where  circled,  the  value  of  the  minor  axis  was  used.  N,  E, 
S,  and  W  indicate  the  names  of  the  triangles. 


in  this  particular  case  by  the  fact  that  the  ice  is 
transported  from  the  Arctic  Ocean  to  the  south, 
not  from  Greenland  to  the  east. 

Deformation.  Figure  7  shows  the  major  princi¬ 
pal  axis  of  Green’s  deformation  tensor,  in  terms 
of  angle  measured  from  north.  The  value  of  the 
minor  principal  axis  was  used  in  some  places 
where  indicated,  so  that  the  angle  does  not  change 
abruptly  by  90°.  The  angles  were  more  or  less  con¬ 
stant,  except  for  triangle  S,  but  at  different  values 
for  each  triangle,  smaller  toward  the  east  and 
south.  The  constant  directions  make  the  interpre¬ 
tation  of  the  principal  values,  which  are  plotted  in 
Figure  8,  convenient.  Exchange  of  major  and  mi¬ 
nor  values  corresponds  to  the  substitution  of  the 
angles  mentioned  above.  The  observation  area 
was  at  first  compressed  in  the  direction  given 
above,  and  then  stretched  for  a  much  greater 
amount  than  before  compression,  either  until  the 
end  of  the  measurement  (N,  E)  or  until  1500  GMT 
(W,  S),  when  the  deformation  in  the  triangle  was 
stopped. 


The  deformation  in  the  minor  direction  indi¬ 
cates  two  different  types.  Triangles  N  and  W  were 
stretched  and  either  immediately  (N)  or  after  com¬ 
pression  (W)  reached  a  constant  value  at  1200  or 
1500  GMT  respectively.  The  second  type  (E,  S) 
was  more  complicated:  small  extension,  large 
compression,  large  extension  then  compression  to 
reach  the  same  constant  value  as  the  first  group. 

There  seems  to  be  a  zone  where  the  extension 
prevailed  in  the  direction  130°  to  170°  from  north, 
the  angle  being  smaller  as  the  point  moved  to  the 
east.  Either  this  zone  was  periodically  compressed 
for  a  short  time,  or  the  observation  area  entered 
this  zone  at  1130.  Otherwise,  the  entire  zone  was 
stretched  in  that  direction  monotonically.  Tri¬ 
angles  W  and  S  left  the  zone  at  1500.  At  right 
angles  to  this  direction  the  zone  was  also  expand¬ 
ed,  but  to  a  lesser  amount;  the  northwestern  part 
was  expanded  more  at  first,  while  the  southeastern 
part  of  the  zone  was  compressed,  then  the  roles 
were  exchanged.  It  was  a  kind  of  wave  propaga¬ 
tion. 
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As  a  result,  the  areas  of  the  triangles  changed, 
as  shown  in  Figure  9.  The  increase  of  the  area 
started  almost  at  the  same  time  (1130,  only  S  at 
1330)  and  the  rates  were  also  similar  to  each  other, 
but  the  increase  was  stopped  at  different  times,  W 
at  1430,  N  at  1500,  and  E  and  S  at  1530.  The  area 
was  increased  to  double  or  more  through  the  stay 
in  the  zone  for  several  hours,  hence  the  ice  concen¬ 
tration  fell  to  less  than  half  of  its  previous  value. 

Motion.  A  similar  analysis  was  made  for  the 
strain  rate  as  for  the  strain,  but  nothing  very  inter¬ 
esting  was  found.  The  reason  for  this  disappoint¬ 


ing  finding  must  be  the  too  large  time  interval  in 
the  measurement.  For  movement  this  rapid,  the 
linearity  assumption  for  30  minutes  was  not  al¬ 
lowed. 

Principal  stretchings  are  shown  in  Figure  10. 
The  almost  linear  decrease  in  the  major  stretching 
of  triangle  E  and  the  large  repeated  areal  change 
rate  of  triangle  S  are  some  of  the  phenomena  that 
it  would  have  been  interesting  to  study  in  detail,  if 
more  information  had  been  available. 

Directly  measured  maximum  stretch  in  I  hr  was 
1 .291 ,  and  the  computed  maximum  stretching  was 
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Fig.  9.  Second  invariant  of  Green's  deformation  tensor,  C,C2 
X,!X,!  =  (area)1.  N,  E,  S,  and  W  indicate  the  names  of  the  triangles. 
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a.  Double  arrow — basic  move¬ 
ment  for  entire  ice  field;  solid 
arrow — additional  movement 
available  only  in  the  boundary 
zone;  curved  arrow — path  line 
of  the  observation  area  rela¬ 
tive  to  the  interior  zone. 


b.  i— interior  zone; 
b — boundary  zone; 
c — chaotic  zone; 
p— front  of  interior 
zone;  p ' — second 
front;  q — ice  edge. 


Fig.  II.  Mode!  of  marginal  ice. 
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0.652  h'1.  In  comparison  to  the  values  obtained  in 
AIDJEX  (e.g.,  Pritchard,  1974)  and  the  MIZEX- 
83  interior  ice  (e.g.,  Hibler  and  Lepparanta, 
1984),  the  active  motion  at  the  ice  margin  is  clearly 
indicated. 

Model 

The  model  was  an  attempt  to  reconstruct  the 
movement  of  the  ice  field,  explaining  all  the  phe¬ 
nomena  observed.  With  so  little  information  in 
hand,  great  freedom  is  left  in  constructing  a  model 
of  the  marginal  ice.  A  plausible  structure  is  pre¬ 
sented  below,  while  many  different  schemas  might 
suit  the  results  as  well. 

The  entire  ice  field  moved  approximately  in  the 
direction  of  the  wind.  It  should  not  be  understood 
from  this  that  the  coriolis  force  can  be  ignored. 
Many  factors  together,  the  wind  being  only  one  of 
them,  moved  the  field  in  a  direction  that  happened 
to  be  very  close  to  the  wind  direction.  Immediately 
at  the  ice  edge,  there  was  a  boundary  zone  where 
additional  southbound  movement  was  overlapped 
(Figure  11a).  The  additional  movement  was  vari¬ 
able  within  the  zone:  greater  speed  toward  the  ice 
edge  and  toward  the  south.  The  increase  of  the 
speed  to  the  south  was  necessary  to  maintain  con¬ 
tinuity  of  ice  transport  in  the  zone,  as  the  zone  be¬ 
came  narrow  to  the  south  (Figure  6b).  The  direc¬ 
tion  was  not  exactly  parallel  to  the  ice  edge  but 
was  slightly  inclined  to  the  east.  The  tilt  angle  was 
smaller,  because  the  point  was  located  closer  to 
the  ice  edge.  Nevertheless  the  eastward  component 
was  greater  toward  the  ice  edge  due  to  increasing 
speed  (Figure  5a). 


The  observation  area  was  in  the  interior  zone  at 
the  beginning,  and  was  transferred  over  line  p  into 
the  boundary  zone  at  1 1 30.  The  transfer  was  not 
made  continuously  but  as  a  block.  Through  such 
an  intermittent  transfer,  linep  retreated  relative  to 
the  material  movement.  Until  the  southward  addi¬ 
tional  component  increased  sufficiently  to  over¬ 
come  the  northward  basic  component  at  1330,  the 
observation  area  continued  to  move  northward. 
The  turning  of  the  path  line  (Figure  I)  thus  hap¬ 
pened  in  the  middle  of  the  boundary  zone  and  had 
no  direct  relation  to  the  crossing  of  a  zone  bound¬ 
ary. 

The  interior  zone  rotated  as  a  whole  counter¬ 
clockwise,  so  that  its  front,  line  p,  came  to  be  at 
right  angles  to  its  movement.  On  the  other  hand, 
the  boundary  zone  rotated  clockwise  so  line  p  ran 
parallel  to  the  ice  edge,  line  q  (Figure  4).  The 
change  of  the  rotation  was  encountered  immedi¬ 
ately  upon  crossing  of  line  p. 

All  the  triangles  were  stretched  in  the  direction 
of  the  boundary  zone  movement  (Figure  8). 
Through  rigid  rotation  clockwise,  triangles  W  and 
S  went  back  to  the  interior  zone  at  1 500,  when  the 
stretch  stopped  increasing.  Along  line  p  a  strong 
local  shearing  force  existed,  which  is  noticed  in  the 
rotation  of  an  individual  floe  (Figure  4),  from 
1200  (eastward  crossing)  and  from  1530  (westward 
crossing). 

The  boundary  zone  was  also  stretched  at  right 
angles  to  the  main  direction  discussed  above,  due 
to  the  increase  of  the  eastward  component  toward 
the  ice  edge.  As  the  eastward  component  was 
much  smaller  than  the  southward  one,  the  corre- 
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sponding  stretch  was  also  smaller  (Figure  8,  minor 
values).  The  behaviour  of  the  minor  values  also  in¬ 
dicates  that  the  shape  of  lines  p  and/or  q  was  not 
steady  but  a  wave  propagated  along  them. 

Where  lines  p  and  q  meet,  a  chaotic  zone  is  ex¬ 
pected  (Figure  11b).  The  irregular  behaviour  of 
triangle  S  may  be  explained  by  the  fact  that  it  lay 
closest  to  the  chaotic  zone  and  was  influenced  by 
the  turbulence  more  than  the  other  triangles.  The 
seaward  extension  of  the  chaotic  zone  may  have 
the  shape  of  an  eddy.  The  boundary  zone  in  this 
model  cannot  be  a  boundary  jet  in  the  strict  sense, 
as  it  is  not  a  continuous  stream  but  a  chain  of  re¬ 
peated  schemata.  The  extension  of  line  p  (shown 
as  p'  in  the  figure)  may  form  a  secondary  front. 

No  comment  can  be  given  from  a  single  short 
observation,  even  if  the  model  is  a  correct  one,  as 
to  whether  this  kinematic  situation  is  a  common 
occurrence  for  the  area.  Furthermore,  kinematic 
study  alone  cannot  tell  through  which  mechanism 
the  movement  was  established. 

Conclusions 

Active  marginal  ice  was  found  very  close  to  the 
ice  edge,  which  was  a  couple  of  kilometers  wide 
and  had  a  fine  structure.  Observations  have  to  be 
made  there,  and  in  detail.  The  behaviour  of  the 
marginal  ice  was  also  different  along  the  ice  edge. 
One-dimensional  analysis  and  related  measure¬ 
ments  should  be  abandoned.  If  the  strain  rate  has 
to  be  measured,  the  time  interval  of  succeeding 
readings  must  be  short. 
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It  has  become  fashionable  to  describe  sea  ice  in 
the  marginal  ice  /one  as  drifting  freely,  mainly  in 
response  to  wind  current  and  wave  effects.  Obser¬ 
vations  of  the  deformation  field  and  associated 
currents  taken  during  MIZEX  83,  however,  give 
quite  another  picture  in  which  the  internal  ice 
stress  is  not  insignificant.  In  particular  during  the 
ice  station  phase  of  the  MIZEX  83  East  Greenland 
experiment,  the  pack  ice  was  quite  compacted  and 
the  resulting  deformation  field  was  rather  damped. 
If  ice  interaction  played  no  role,  we  would  expect 
the  ice  to  continue  converging  forever  under  com¬ 
pacting  winds.  What  seems  to  happen  instead  is 
that  once  the  ice  reaches  a  certain  compact  state,  it 
ceases  converging  due  to  the  ice  interaction  and 
there  is  a  concomitant  reduction  in  the  amount  of 
shearing.  Moreover,  initial  analysis  of  MIZEX  84 
data  shows  this  more  compact  state  to  be  associ¬ 
ated  with  both  greatly  reduced  fluctuations  in  the 
deformation  signal  and  less  floe  rotation. 

During  the  1983  experiment  we  were  fortunate 
enough  to  have  simultaneous  current  measure¬ 
ment''  at  the  three  outer  vertices  of  the  deforma¬ 
tion  array  (O,  C,  and  J  in  Figure  1).  Currents  were 
measured  at  several  levels  at  10-minute  intervals, 
and  deformation  was  recorded  at  3-minute  inter¬ 
vals  with  a  better-than-l-m  accuracy  using  a  Del 
Norte  radar  transponder  system  (for  a  more  com¬ 
plete  description  of  the  deformation  data,  see  Hib- 
ler  and  Lepparanta,  1983).  Over  much  of  this  peri¬ 
od  the  winds  were  relatively  weak,  so  currents  sup¬ 
plied  the  main  motive  source  for  the  ice.  More¬ 
over,  analysis  of  the  relative  currents  between  the 


floes  shows  the  ice  typically  remains  “locked  up,” 
even  though  the  relative  currents  can  be  substan¬ 
tial  and  even  though  the  whole  region  may  be 
drifting  in  an  absolute  sense  with  the  currents. 

I _ I _ I _ I _ I _ I 

0  2  4  6  8  10  km 


Fig.  I.  Configuration  of  strain  array  and  no¬ 
menclature  for  distances  discussed.  The  triangle 
is  plotted  to  scale  for  29  June  0800  GMT;  the  ship 
(M/V  PolarbjOrnV  is  located  at  the  center  point  P 
moored  to  an  ice  floe.  Del  Norte  units  are  located  at 
solid  points.  Directly  measured  distances  (accurate  to 
1-2  m)  are  denoted  by  subscripted  b's  and  d's  and 
distances  calculated  by  triangulation  < less  accurate 
depending  on  geometry)  are  denoted  by  subscripted 
s  \. 
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Fig.  2.  Time  series  of  smoothed  relative  currents  and  relative 
ice  motion  between  stations  O  and  C. 


Such  behavior  is  quite  consistent  with  plastic 
models  (see,  for  example  Lepparanta  and  Hibler, 
1985),  but  not  with  simple  pressure-term-only  rhe¬ 
ologies  (e.g.  Rried  and  O’Brien,  1983).  Within  the 
framework  of  a  plastic  rheology,  this  situation  af¬ 
fords  enough  information  for  a  lower  bound  esti¬ 
mate  of  the  stress  in  the  ice.  Basically,  under  rela¬ 
tive  rigid  conditions  the  rigid  state  of  plasticity 
takes  the  appropriate  stress  state  necessary  to  ab¬ 
sorb  differences  in  the  body  forces  (in  this  case  the 
water  drag),  as  long  as  these  stress  states  do  not 
exceed  the  plastic  limit. 

A  typical  situation  of  this  nature  that  occurred 
is  shown  in  Figure  2  where  we  have  plotted  the  ice 
motion  and  ocean  current  differences  between 
sites  C  and  O  over  a  period  of  several  days.  These 
data  have  been  somewhat  smoothed  by  filtering 


but  that  does  not  affect  the  argument.  As  can  be 
seen  from  this  figure,  the  relative  currents  are 
quite  substantial  whereas  the  ice  is  essentially 
locked  up.  This  “locked  up’’  state  is  hard  to 
justify  by  different  drag  coefficients,  since  the 
same  relative  rigid  conditions  occur  for  both 
divergent  and  convergent  motion. 

Based  on  this  figure  the  relative  currents  have  a 
maximum  value  on  the  order  of  5  cm  s  ',  and  the 
distance  between  the  sites  is  on  the  order  of  10  km. 
If  a,  is  the  stress  at  one  site,  then,  assuming  all  var¬ 
iables  vary  linearly  over  the  deformation  array,  to 
prevent  ice  motion  must  be  at  least  a,  +  1<> 
where  -W.lv  equals  the  differential  ice  water 
stress  between  the  two  sites.  To  obtain  the  water 
stress  precisely  would  require  knowledge  of  the 
absolute  currents,  but  as  a  crude  estimate  wc  take 


18 


the  absolute  current  to  be  0  at  one  site.  Therefore, 
taking  a  quadratic  water  coefficient  of  0.0055,  the 
water  stress  differential  is 

0.0055gw(0.04);  =  0.0055  •  10'.  1.6*  10  5 
=  8.8  •  10  •’  Nt  m‘! 

.-.  Ao  =  8.8  •  10  2  •  10J 
=  10’  Nt  nr1. 

While  this  value  is  smaller  than  typical  strengths 
used  in  large-scale  plastic  models  (see  e.g.  Hibler 
and  Walsh,  1980)  on  the  order  of  10*  to  10'  Nt 
m  ',  it  is  by  no  means  trivial.  It  is  clear  from  the 
estimation  procedure  that  this  value  depends  on 
the  spacing  between  the  points.  To  support  the 
same  current  differential  over  a  larger  region,  the 
stress  must  vary  a  larger  amount.  This  is  in  fact 
one  reason  why,  on  a  larger  scale,  a  plastic  system 
will  deform  as  if  it  were  in  free  drift,  whereas 
smaller  sections  can  be  relatively  rigid. 


While  the  estimates  made  here  are  crude,  they 
do  give  an  idea  of  the  plastic  stresses  that  can  be 
expected  in  the  marginal  ice  zone.  These  compari¬ 
sons  also  indicate  that  ice  deformation  provides  a 
much  more  sensitive  test  between  ice  drift  theories 
than  ice  trajectories. 
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During  June  and  July  1984,  sea  ice  physical  pie  of  our  characterization  method  applied  to  a 

property  investigations  were  conducted  from  the  core  taken  through  undeformed  multi-year  ice.  In 

FS  Polarstern  during  MIZEX-84.  A  large  area  this  paper,  we  wish  to  report  on  specific  examples 

within  the  Fram  Strait  was  traversed  by  Polar-  of  the  crystal  structures  observed  on  MIZEX  and 

stern,  providing  an  opportunity  to  obtain  core  to  discuss  the  origin  of  these  textures  with  regard 

samples  that  very  likely  had  originated  in  different  to  the  ice  growth  and  ablation  processes.  First, 

parts  of  the  Arctic  basin.  Ice  cores  were  taken  at  however  a  general  description  of  sea  ice  growth  is 

the  locations  shown  in  Figure  1.  For  each  core  we  presented  as  background  for  later  discussion  of 

measured  temperature  and  salinity  profiles  and  the  crystal  structure. 

documented  the  crystal  structure  through  the  Sea  ice  formation  typically  begins  with  the  free 

thickness  of  the  ice  floe.  Figure  2  shows  an  exam-  growth  of  crystals  in  a  slightly  supercooled  water 


Pig.  I.  Ice  properties  sampling  locations.  Approximate  ice  edge  locations  for  19  June  and  17  July  19X4  are  also  shown 


Fif(.  2.  Typical  properties  and  structural  profile  of 
multi-year  sea  ice  composed  of  90%  congelation. 
Columnar  ice  overlain  by  granular  ice  and  a  refrozen 
melt  layer. 


column.  This  initial  growth  usually  occurs  in  the 
form  of  minute  spheres,  which  rapidly  change  to 
discs  that  typically  grow  to  diameters  of  2  to  3  mm 
(Weeks  and  Ackley,  1982).  The  discs  then  evolve 
into  hexagonal  dendritic  star-like  crystals  that  are 
much  more  effective  at  dissipating  heat  and  solute 
into  the  surrounding  liquid.  If  the  water  remains 
calm,  the  crystals  tend  to  float  with  their  basal 
planes  parallel  to  the  water  surface  with  the  result 
that  when  they  freeze  together  to  form  a  surface 
skim,  the  crystals  become  oriented  with  their 
c-axes  vertical.  This  is  rarely  the  case,  however, 
because  freezing  most  commonly  occurs  under 
turbulent  conditions  (wind-  or  wave-induced)  that 
cause  many  crystals  to  be  trapped  with  their  c-axes 
at  various  inclinations  from  the  vertical.  Generally 
this  surface  layer  of  frazil  or  granular  ice  grows  no 
thicker  than  a  few  centimeters,  but  the  action  of 
stronger  winds  and  waves  can  result  in  accumula¬ 
tions  of  frazil  up  to  I  m  thick. 

Once  a  stable  surface  layer  of  ice  is  established, 
growth  of  columnar  crystals  begins.  This  appear¬ 


ance  of  a  new  crystal  form  marks  the  transition  to 
congelation  growth,  the  direct  freezing  of  sea¬ 
water  to  the  bottom  of  the  existing  ice  sheet.  This 
transition  generally  occurs  within  5  to  10  cm  of  the 
surface,  but  in  some  instances  may  extend  to  50 
cm  deep. 

Congelation  growth  typically  produces  vertical¬ 
ly  elongated  columnar  crystals  that  can  be  tens  of 
centimeters  long  but  usually  do  not  exceed  2-3  cm 
in  cross-sectional  diameter.  The  c-axes  of  these 
crystals  themselves  also  possess  a  cellular  sub¬ 
structure  consisting  of  vertical  plates  0.5-1  mm 
thick,  separated  by  parallel  layers  of  brine  inclu¬ 
sions.  The  platelets  consist  of  essentially  pure  ice 
with  the  characteristic  salinity  arising  from  the 
brine  incorporated  between  the  plates.  Such  struc¬ 
ture  reflects  the  morphology  of  the  growing  inter¬ 
face,  which  consists  of  pure  ice  dendrites  protrud¬ 
ing  vertically  downwards  into  the  sea  water.  Some 
of  the  brine  excluded  during  freezing  is  trapped  in 
the  grooves  between  the  dendrites  as  the  growth 
interface  advances.  The  initial  salinity  of  the  ice, 
which  is  determined  by  the  amount  of  entrapped 
brine,  has  been  shown  to  be  inversely  proportional 
to  the  brine  layer  spacing  (Nakawo  and  Sinha, 

1984) .  The  layer  spacing,  in  turn,  is  related  to  the 
growth  rate  such  that  the  faster  the  growth,  the 
closer  the  spacing  and  thus  the  saltier  the  ice,  and 
vice  versa.  Once  the  initial  salinity  of  the  ice  has 
been  established,  subsequent  changes  are  deter¬ 
mined  essentially  by  the  thermal  history  of  the  ice. 
Temperature  variations  cause  changes  in  the  ge¬ 
ometry  and  distribution  of  the  brine  inclusions. 
The  largest  effect  is  the  formation  of  brine  drain¬ 
age  systems,  vertical  tubes  with  diagonal  feeder 
channels  made  up  of  coalesced  brine  inclusions. 
The  formation  of  these  drainage  systems  is  always 
associated  with  significant  warming  of  the  ice. 
Changes  in  crystalline  texture  may  also  accom¬ 
pany  desalination,  and  the  combination  of  all 
these  effects  generally  leads  to  significant  modifi¬ 
cation  of  the  mechanical  and  electromagnetic  pro¬ 
perties  of  the  ice. 

What  we  have  described  above  is  the  classical 
growth  sequence  of  sea  ice.  However,  there  can  be 
deviations  from  this  norm,  including  the  occur¬ 
rence  of  frazil  or  granular  ice  at  levels  in  the  ice 
that  are  unrelated  to  surface  turbulence  effects. 
Our  preliminary  MIZEX  findings  (Tucker  et  al., 

1985)  indicated  that  24%  of  the  total  ice  sampled 
consisted  of  granular  ice,  mainly  frazil,  but  also 
including  snow  ice,  necessarily  limited  to  the  sur¬ 
faces  of  floes.  Most  of  this  frazil  was  found  in 
cores  taken  from  previously  deformed  ice,  e  g.,  in 
ridged  ice.  Generally,  granular  ice  was  restricted 
to  frazil  occupying  the  voids  between  the  inclined 
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Fig.  3.  Horizontal  thin  section  of  fine-grained  frazil  from  a  depth  of  300  cm  in  a  multi-year  ridge.  Scale  meas- 
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Fig.  4.  Horizontal  thin  section  of  snow  ice  at  3-cm  depth  in  first-year  ice.  Scale  measures  in  centimeters  and  milli¬ 
meters. 
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blocks  in  the  ridges.  We  speculate  that  freezing 
nuclei,  consisting  of  snow  or  crushed  ice,  precipi¬ 
tated  the  formation  of  frazil  crystals  in  the  voids. 
In  constrast  to  the  small  amounts  of  frazil  ob¬ 
served  on  M1ZEX,  Gow  et  al.  (1982)  found  large 
accumulations  of  this  ice  type  in  Weddell  Sea  ice. 
Frazil  ice  made  up  60%  of  the  total  ice  sampled 
and  was  observed  in  both  first-  and  multi-year  ice 
and  occasionally  in  thicknesses  exceeding  3  m.  Ap¬ 
parently,  environmental  conditions  linked  closely 
to  the  depth  and  structural  characteristics  of  the 
oceanic  mixed  layer  are  more  amenable  to  frazil 
formation  in  the  Weddell  Sea  than  in  the  Arctic 
basin. 

We  have  used  both  horizontal  and  vertical  thin 
sections  to  examine  ice  structure.  Thin  sections 
were  prepared  from  cores  using  microtome  sec¬ 
tioning  techniques  described  in  Weeks  and  Gow 
(1978).  The  structural  plates  shown  here  were  ob¬ 
tained  by  placing  the  thin  sections  between  crossed 
polaroids  and  photographing  them  in  transmitted 
light.  This  technique,  which  utilizes  the  principle 
of  interference  colors,  clearly  reveals  individual 
crystal  boundaries,  brine  drainage  channels,  and 
the  brine  lamallae/platelet  substructure  of  individ¬ 
ual  crystals.  Color  photographs  (slides)  were  orig¬ 


inally  made  with  a  35-mm  camera  using  Koda- 
chrome®  200  ASA  color  slide  film  at  aperture  set¬ 
tings  ranging  from  F8  to  F126  at  a  1-second  shut¬ 
ter  speed.  The  slides  were  then  used  to  prepare  the 
color  internegatives  needed  to  produce  the  color 
prints  shown  here.  Black  and  white  photographs 
were  taken  with  a  10-  by  12.5-cm  press  camera  set 
up  to  provide  negatives  at  approximately  a  1:1 
scale.  The  Kodak  Plus-X-Pan®  film  used  for  this 
purpose  was  exposed  for  about  0.5  second  at  an 
aperture  setting  of  F8. 

Figures  3  and  4  illustrate  two  varieties  of  granu¬ 
lar  ice.  The  ice  in  Figure  3  consists  of  very  fine¬ 
grained  frazil  that  was  located  at  a  depth  of  300 
cm  in  a  multi-year  floe.  Inclined  columnar  ice 
both  above  and  below  the  frazil  indicate  that  it 
formed  in  the  void  between  blocks  in  a  pressure 
ridge.  Figure  4  features  a  much  coarser-grained 
type  of  ice  consisting  of  crystals  with  polyhedral 
outlines.  This  is  snow  ice,  formed  by  the  freezing 
of  water-saturated  snow.  This  is  further  con¬ 
firmed  by  its  occurrence  at  the  top  of  the  floe  and 
by  the  fact  that  its  salinity  was  less  than  0.5%o. 

Refrozen  melt  ponds  were  found  to  have  crystal 
textures  that  are  very  similar  to  those  observed  in 
ice  formed  on  freshwater  lakes.  Figure  5  shows 


Fig.  5.  Horizontal  thin  section  of  melt  pond  ice  from  7-cm  depth  in  multi-year  sea  ice.  Scale  measures  in  centi¬ 
meters  and  millimeters. 


typical  melt  pond  crystal  texture  as  observed  in  a 
horizontal  section.  Grain  cross-sectional  diam¬ 
eters  up  to  1  cm  are  commonplace  and  the  c-axes 
are  randomly  distributed  in  the  horizontal  plane. 
The  location  of  this  section  was  7  cm  below  the 
top  of  the  melt  pond  on  a  multi-year  floe.  The  sa¬ 
linity  at  this  depth  measured  less  than  0.2%,.  The 
vertical  section  of  this  ice  (not  shown  here)  dis¬ 
played  a  columnar  texture  with  vertically  oriented 
crystals.  A  much  more  dramatic  example  of  melt 
pond  structure  is  shown  in  Figure  6,  a  vertical  sec- 


Fig.  6.  Vertical  thin  section  through  a  melt  pond 
from  the  ice  surface  to  a  depth  of  15  cm  in  multi¬ 
year  ice. 


tion  of  the  top  15  cm  from  a  different  multi-year 
floe.  Here  a  5-cm-thick  layer  of  snow  ice,  made  up 
of  large,  randomly  oriented  crystals,  overlies  a 
melt  pond  exhibiting  c-axis  vertical  melt  crystal 
structure.  The  individual  crystals  here  are  very 
large,  the  result  of  slow  freezing  of  the  water.  Gow 
(1984)  points  out  that  the  type  of  crystal  structure 
occurring  in  lake  ice  covers  is  determined  by 
whether  or  not  the  lake  water  is  seeded  at  the  onset 
of  freezing.  Seeding  leads  to  the  growth  of  colum¬ 
nar,  r-axis  horizontal  crystals,  and  unseeded 
(spontaneously  nucleated)  growth  promotes  the 
formation  of  massive,  c-axis  vertical  crystals.  The 
same  is  also  true  of  melt  ponds  freezing  on  the 
tops  of  multi-year  sea  ice  floes.  Beneath  the  melt- 
pond  ice,  a  sharp  transition  to  columnar  sea  ice  is 
evident. 

As  with  granular  ice,  a  variety  of  columnar  ice 
forms  were  observed  on  M1ZEX.  Figure  7  shows 
columnar  ice  in  horizontal  section  in  which  the 
c-axes  are  not  only  horizontal  but  are  all  very 
nearly  aligned  in  the  same  direction.  (In  columnar 
sea  ice  the  c-axes  are  normal  to  the  long  dimen¬ 
sions  of  the  crystals.)  This  particular  section  was 
prepared  from  multi-year  ice  at  a  depth  of  203  cm. 
In  sharp  contrast,  Figure  8  shows  a  section  from 
near  the  bottom  of  a  236-cm-thick  first-year  floe. 
In  addition  to  exhibiting  a  variety  of  crystal  sizes 
and  shapes,  very  little  crystal  alignment  is  observed 
here.  Weeks  and  Gow  (1978,  1980),  on  the  basis  of 
measured  and  inferred  current  direction,  believe 
that  the  alignment  of  crystal  c-axes  is  controlled 
by  the  direction  of  the  current  at  the  growing  in¬ 
terface.  This  control  of  current  directions  on  the  c- 
axis  alignments  in  sea  ice  has  now  been  verified  in 
a  series  of  laboratory  experiments  by  Langhorne 
(1983).  A  likely  explanation  is  that  a  current  mov¬ 
ing  across  (parallel  to  the  c-axes),  rather  than  par¬ 
allel  to  the  protruding  dendrites  (normal  to  the 
c-axes)  presumably  provides  more  effective  heat 
transfer  away  from  the  dendrites,  resulting  in  a 
growth  advantage.  Strong  horizontal  alignments 
would  indicate  that  the  ice  had  remained  in  a  fixed 
position  relative  to  the  steady  current  direction.  A 
weak  alignment  or  lack  of  alignment  implies  that 
the  ice,  if  in  fixed  position,  must  have  grown  in  an 
area  of  highly  variable  current  directions  or,  in  the 
case  of  drifting  pack  ice,  must  have  been  rotating 
as  it  drifted. 

Figure  9  shows  much  finer  grained  columnar  ice 
than  either  of  the  previous  two  figures.  We  do  not 
fully  understand  all  the  processes  controlling  the 
crystal  size  of  columnar  sea  ice,  but  we  speculate 
that  the  freezing  rate  of  the  sea  water  must  be  an 
important  factor.  Note  also  in  this  fine-grained  ice 
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Fig.  9.  Horizontal  thin  section  of  fine-grained  columnar  ice  from  a  depth  of  225  cm  in  multi-year  ice.  Scale 
measures  in  centimeters  and  millimeters. 


ures  in  centimeters  and  millimeters. 
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the  absence  of  any  significant  alignment  of  the 
crystals.  This  section  of  ice  was  obtained  from  a 
depth  of  225  cm  in  a  multi-year  ridge  fragment. 

One  of  the  more  interesting  structures,  which  to 
our  knowledge  has  not  been  previously  reported, 
is  shown  in  Figure  10.  Ice  in  this  thin  section,  from 
a  depth  of  15  cm  in  a  multi-year  floe,  shows  con¬ 
siderable  modification  of  its  structure  due  to  ele¬ 
vated  temperatures  during  the  previous  summer. 
The  more  obvious  structural  modifications  in¬ 
clude  nearly  complete  obliteration  of  the  original 
ice  platelet/brine  lamallae  substructure  and  sub¬ 
stantial  smoothing  of  the  grain  boundaries.  Al¬ 
though  there  appear  to  be  no  significant  changes 
in  the  overall  crystal  lattice  structure,  as  evidenced 
by  the  retention  of  an  original  c-axis  alignment, 
this  retexturing  always  accompanies  other  changes 
in  the  physical  properties  of  the  ice.  For  example, 
it  1)  is  invariable  associated  with  ice  of  very  low 
salinity  (less  than  0.5%o),  and  2)  is  restricted  to  the 
top  meter  of  multi-year  ice.  The  ice  is  much  less 
opaque  (often  semi-transparent)  in  appearance 
and  is  much  more  resistant  to  drilling  and  sawing 
than  ice  from  the  same  depth  in  a  first-year  floe. 


Although  the  sequence  of  horizontal  thin  sec¬ 
tions  in  the  multi-year  profile  in  Figure  2  clearly 
indicates  some  diversity  of  crystal  structure,  evi¬ 
dence  for  retexturing  of  the  ice  is  best  observed  in 
magnified  sections,  such  as  those  shown  in  Figure 
11.  These  sections  were  obtained  from  depths  of 
20,  93,  and  192  cm  respectively,  in  a  294-em-thick 
floe  that  was  very  similar  structurally  to  that 
shown  in  Figure  2.  Vertical  sectioning  of  the  floe 
verified  that  it  consisted  entirely  of  congelation 
(columnar)  ice. 

In  ice  from  192  cm  (Figure  1 1 c),  crystals  are  still 
interpenetrated  structurally  and  still  retain  the 
well-defined  platelet/brine  lamalla  substructure  so 
characteristic  of  congelation  ice.  Flowever,  the  sa¬ 
linity—  2. 5%o — clearly  indicates  some  loss  of  the 
originally  entrapped  brine  even  at  this  level  in  the 
floe.  Since  this  floe  was  294  cm  thick,  it  is  very 
likely  that  ice  now  at  192  cm  was  formed  two  win¬ 
ters  ago  but  was  located  deep  enough  in  the  floe  to 
protect  it  from  the  dramatic  effects  of  surface 
warming  during  the  intervening  summer  (1983). 
However,  this  situation  is  not  true  of  ice  at  93  cm 
(Figure  1  lb)  where  only  vestiges  of  the  original  ice 


platelet /brine  layer  substructure  of  crystals  are  re¬ 
tained  (though  original  t-axis  alignment  appears 
to  be  preserved),  and  the  salinity  has  been  reduced 
to  l.7‘7oo,  indicative  of  substantial  desalination. 
The  situation  is  even  more  extreme  in  the  case  of 
ice  from  20  cm  (Figure  1  la),  which  bears  little  re¬ 
semblance  whatsoever  to  congelation  sea  ice. 
Crystal  boundaries  have  developed  sutured  out¬ 
lines,  and  all  trace  of  substructure  has  disap¬ 
peared.  Total  loss  of  brine  pocket  structure  in  this 
ice  is  also  reflected  in  the  salinity,  which  measured 
less  than  0.2%o. 

Retexturing  of  crystals  in  the  upper  levels  of 
multi-year  floes  and  the  associated  loss  of  brine 
(and  brine  pockets)  also  produces  very  substantial 
changes  in  the  mechanical  and  electrical  properties 
of  the  ice.  In  the  case  of  the  latter  such  changes 
should  result  in  remote  sensing  signatures  being 
significantly  different  from  those  of  first-year  ice. 

In  this  brief  report  we  have  described  the  vari¬ 
ous  forms  of  crystal  structure  that  characterize  the 
sea  ice  in  Fram  Strait.  We  have  found  that  frazil  is 
a  minor  component  of  the  floes,  and  that  retextur¬ 
ing  of  the  crystal  structure  is  a  prevalent  and  diag¬ 
nostic  property  of  the  congelation  ice  in  the  top 
meter  or  so  of  multi-year  floes.  Based  on  the  struc¬ 
tural  properties  described  here,  in  conjunction 
with  other  physical  properties  of  the  ice,  we  have 
been  able  to  distinguish  between  first-  and  multi¬ 
year  ice  and  to  determine  the  ages  of  multi-year 
floes.  In  addition,  we  have  obtained  a  much  better 
understanding  of  the  processes  that  control  the 
growth  characteristics  and  metamorphic  proper¬ 
ties  of  the  ice  entering  Fram  Strait  from  the  Arctic 
basin. 

The  authors  very  much  appreciate  the  continu¬ 
ing  support  provided  by  fhe  Arctic  Program,  Of¬ 
fice  of  Naval  Research,  which  allows  us  to  pursue 


these  studies.  The  logistical  support  provided  by 
the  scientific  staff  and  crew  of  the  FS  Polarslern  is 
gratefully  acknowledged.  Additional  funds  were 
provided  to  the  MIZEX  Bulletin  to  cover  costs  as¬ 
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INTRODUCTION 

The  principal  parameters  used  in  modeling  the 
microwave  response  of  a  surface  are  the  dielectric 
constant  of  the  surface  material  and  the  surface 
and  subsurface  roughness.  The  dielectric  constant 
depends  on  the  physical  properties  of  the  material: 
temperature,  density,  and  composition.  It  is  great¬ 
ly  affected  by  the  presence  of  liquid  water  in  the 
material  because  of  the  highly  absorptive  nature 
of  water  in  the  microwave  portion  of  the  spec¬ 
trum.  In  the  case  of  snow  cover,  snow  wetness  can 
be  a  critical  factor  in  the  value  of  the  dielectric 
constant,  and  therefore  in  its  microwave  response. 

The  snow  cover  on  arctic  sea  ice  is  of  interest  for 
both  its  own  signature  and  its  effect  on  the  signa¬ 
ture  of  the  underlying  ice.  Surface  albedo  is  a  key 
parameter  in  thermodynamic  models  of  sea  ice 
growth  and  can  differ  significantly  for  dry  snow 
and  melting  snow  conditions  (Shine  and  Hender- 
son-Sellers,  1985).  The  state  of  the  snow  also  de¬ 
termines  whether  it  is  transparent  or  opaque,  or 
totally  masking  to  radiation  coming  from  the  sea 
ice.  Snowpack  with  high  liquid  water  content  can, 
depending  on  the  operating  frequency  of  the  sen¬ 
sor,  absorb  sufficient  radiation  that  it  appears 
similar  to  open  water,  or  limit  penetration  such 
that  surface  scattering  dominates  and  the  sea  ice 
has  a  uniform  signature  regardless  of  ice  type  (On- 
stott  and  Moore,  1984).  Both  effects  introduce  er¬ 
rors  into  estimates  of  ice  concentration  and  ice 


thickness  (i.e.  ice  type)  distribution  derived  from 
remote  sensing  instruments. 

During  MIZEX  84  several  groups  working  from 
different  research  vessels  made  measurements  of 
the  dielectric  constant,  snow  wetness  and  physical 
properties  of  the  snow  cover  on  ice  floes.  We  re¬ 
port  here  results  obtained  by  the  group  operating 
from  the  FS  Polarstern.  Nineteen  sites  were  visited 
during  the  period  18  June-14  July  1984  for  which 
both  dielectric  constant  and  snow  wetness  meas¬ 
urements  were  made.  These  sites  were  distributed 
across  the  region  traversed  by  Polarstern,  gener¬ 
ally  78°10  'N  to  80°40  'N  and  1 1  °W  to  9°E,  as  seen 
in  Figure  1 .  Exact  site  locations  are  given  in  Burns 
and  Larson  (1985).  Snow  cover  conditions  at  these 
sites  ranged  from  dry  and  cold  to  melting,  with 
snow  depths  between  0  and  85  cm.  Ice  cores  were 
obtained  at  most  sites  by  the  Cold  Regions  Re¬ 
search  and  Engineering  Laboratory  (CRREL) 
team,  so  that  ice  type  was  positively  determined. 
Dielectric  constant  measurements  of  two  of  the 
core  samples  are  also  reported  here. 

Although  some  sites,  especially  the  ones  visited 
early  in  the  experiment,  do  not  have  coincident  re¬ 
mote  sensing  data,  the  aim  was  to  coordinate 
measurements  with  shipscat  and  heloscat  opera¬ 
tions  and  with  aircraft  SAR,  SLAR  or  PM1  data 
collection.  Nine  of  the  site  measurements  were  co¬ 
ordinated  with  shipscat  measurements,  five  with 
heloscat  flights,  two  with  aircraft  PMI  overflights, 
and  seven  with  aircraft  SAR  or  SLAR  overflights. 
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Fig.  1.  Site  positions  for  dielectric  constant  measurements  made  from 
Polarstern  during  MIZEX  84. 


EQUIPMENT  AND  OPERATION 

Equipment  used  in  making  surface  measure¬ 
ments  included  ERIM  dielectric  constant  instru¬ 
mentation  at  100  MHz,  L-band  and  X-band  (Lar¬ 
son,  1986),  CRREL  alcohol  calorimetry  apparatus 
to  measure  snow  wetness  (Fisk,  1983),  a  snow  den¬ 
sity  kit  (sample  tubes  and  balance),  stem  ther¬ 
mometers,  and  collapsible  rulers.  The  100-MHz 
Q-meter  proved  the  most  useful  of  the  dielectric 
instrumentation  because  it  was  the  most  portable, 
could  be  used  to  profile  the  significant  snow  cover 
found  throughout  the  experimental  area,  and  pro¬ 
vided  values  for  both  the  real  part  of  the  dielectric 
constant  and  the  loss  tangent  most  directly. 

At  sites  accessed  directly  from  the  ship,  a  full  set 
of  dielectric  constant,  snow  wetness,  and  snow 
characterization  measurements  were  carried  out, 
as  it  was  possible  to  transport  all  necessary  equip¬ 
ment  to  the  ice.  For  these  sites  snow  samples  were 
taken  back  to  the  ship  and  the  liquid  water  frac¬ 
tion  measured  in  a  coldroom  under  controlled 
conditions.  Usually,  but  not  always,  ship- 
mounted  scatterometer  measurements  and  ice 
cores  were  also  obtained  for  these  sites. 

When  accessing  the  ice  with  helicopter,  both 
time  and  space  were  very  limited  so  that  only 
100-MHz  dielectric  measurements  were  carried 
out.  On  one  helicopter  mission  snow  wetness  was 


measured  on  the  floe,  but  due  to  the  lengthiness  of 
the  measurement  this  was  not  repeated,  and  on  the 
next  mission  samples  were  brought  back.  The  in¬ 
evitable  change  of  state  caused  by  transporting 
samples  will  of  course  introduce  additional  uncer¬ 
tainties  in  the  measurement  results. 

The  dielectric  constant  results  reported  here  are 
restricted  to  those  obtained  at  100  MHz,  as  these 
represent  the  most  complete  record  both  spatially 
and  temporally.  A  full  tabulation  of  the  surface 
measurements  is  given  in  Bums  and  Larson 
(1985). 


Snow  cover 

The  principal  reason  for  making  this  set  of  sur¬ 
face  measurements  is  to  gain  a  better  understand¬ 
ing  of  how  the  complex  dielectric  constant,  and 
therefore  remotely  sensed  microwave  signatures, 
can  be  used  to  obtain  information  on  snow  cover 
conditions.  It  is  of  interest  then  how  the  measured 
surface  parameters  varied  throughout  the  experi¬ 
ment  and  how  the  dielectric  constant  is  related  to 
them. 

The  temporal  record  of  observed  snow  wetness 
between  21  June  and  16  July  is  shown  in  Figure  2 
superimposed  on  the  surface  air  temperature 
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Fig.  2.  Variation  of  the  liquid  water  content  of  the  snow  cover  contrast¬ 
ed  with  air  temperature  variation  observed  during  the  experiment. 
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DEPTH  OF  SNOW  COVER  (cm) 

Fig.  3.  Frequency  of  observed  mean  snow  depth  as  a  function  of  ice  type. 


record  from  Polarstern.  L  iquid  water  values  lor  sparse  for  this  trend  to  be  quantified.  It  should  be 

subsurface  layers  show  much  less  extreme  varia-  noted  that  liquid  water  variability  seen  in  Figure  2 

tion  (I  to  l3ro  by  mass)  than  do  the  surface  values  also  includes  spatial  variability,  as  the  ship  did  not 

(0  to  21  ro  by  mass),  indicating  the  insulating  cl-  remain  at  one  floe  for  more  than  a  few  hours,  and 

feet  of  the  snow  .  This  may  aUo  be  indicated  by  the  it  docs  not  show  the  high  liquid  water  content  of 

observed  lag  in  the  liquid  water  plots  relative  to  snow  cover  observed  near  the  snow/icc  interface, 

the  temperature  plot,  which  is  most  pronounced  which  increased  throughout  the  experiment  as  the 
for  the  subsurface,  although  sample  spacing  is  too  melt  season  progressed. 
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Fig.  4.  Observed  variations  in  snow  density,  grain  size  and  air  temperature 
during  the  experiment  period. 


As  seen  from  Figure  2,  the  majority  of  liquid 
water  mass  fraction  values  for  subsurface  meas¬ 
urements  fall  in  the  5  to  10%  range.  The  extreme 
values  outside  of  this  range  correspond  to  first- 
year  floe  measurements.  Although  only  a  few  data 
points  are  available,  different  behaviors  for  first- 
year  and  multiyear  ice  may  be  expected  because  of 
differences  in  amount  of  snow  cover.  On  the  aver¬ 
age,  multiyear  floes  had  twice  as  much  snow  cover 
as  first-year  floes.  The  distribution  of  snow  depths 
observed  on  both  types  of  floes  is  shown  in  Figure 
3;  contributions  to  the  0-  to  10-cm  bin  for  multi¬ 
year  floes  come  mainly  from  those  floes  visited 
late  in  the  experiment  when  melting  was  signifi¬ 
cant. 

The  transition  to  melt  conditions  also  produced 
a  metamorphosis  in  snow  cover  from  the  less 
dense,  small  crystal  snow  to  dense,  large  crystal 
(corn  snow)  firn.  This  was  generally  observed  dur¬ 
ing  the  experiment  and  is  to  some  extent  quanti¬ 
fied  by  grain  size  and  density  measurements  plot¬ 
ted  as  a  function  of  time  in  Figure  4.  Density  val¬ 
ues  plotted  are  for  surface,  or  subsurface  in  cases 
of  a  surface  crust.  Again  these  plots  reflect  spatial 
as  well  as  temporal  variability,  although  over  the 
area  sampled  the  incident  radiation  is  probably 
uniform. 

Given  the  variation  observed  for  liquid  water, 
temperature,  density  and  grain  size  (Figures  2  and 
4),  significant  variation  in  the  dielectric  constant, 
which  depends  on  all  these  factors,  is  also  expected. 
The  temporal  record  of  the  dielectric  constant  at 


100  MHz  is  shown  in  Figure  5  for  both  the  snow 
surface  and  snow  at  depths  of  2-15  cm.  Values 
range  from  1.5  for  frozen  crust  (Ta u  =  -1.5°C)  to 
4.5  for  firn  with  5%  liquid  water  content  by  vol¬ 
ume.  The  surface  measurements  show  a  close  cor¬ 
relation  with  surface  temperature  fluctuations  (see 
Figure  6),  whereas  the  subsurface  measurements 
indicate  more  stable  dielectric  properties  at  depth. 
The  spatial  variability  expected  is  indicated  in  Fig¬ 
ure  7,  which  shows  contours  of  dielectric  constant 
based  on  a  grid  of  measurements  made  on  a  site 
with  fairly  advanced  melt.  These  results  have  two 
major  implications  for  remotely  sensed  microwave 
signatures:  (1)  the  microwave  signatures  will  show 
significant  temporal  variability  if  the  frequency  is 
such  that  only  the  surface  is  sensed;  and  (2)  given 
the  typical  sensor  footprint  size,  the  response 
reflects  a  mean  dielectric  constant  with  possibly 
large  variance. 

To  understand  the  microwave  signatures  in 
terms  of  snow  cover  conditions  (or  in  other  words 
the  effect  of  snow  cover  on  the  signatures)  re¬ 
quires  examination  of  the  relationship  between  the 
complex  dielectric  constant  and  snow  cover 
parameters.  During  the  summer  season,  the  dielec¬ 
tric  properties  of  the  snow  cover  are  primarily  de¬ 
termined  by  its  liquid  water  content.  The  correla¬ 
tion  between  the  electrical  properties  and  liquid 
water  content  is  demonstrated  in  Figure  8.  Here 
the  real  part  of  the  dielectric  constant,  f ',  and  the 
loss  tangent,  tan5,  are  plotted  against  volume  li¬ 
quid  water  content  per  unit  volume,  Wv.  A  linear 
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Fig.  5.  Variations  in  dielectric  constant  values  for  surface  and  subsur¬ 
face  snow  contrasted  with  air  temperature  variation  for  the  experiment 
period. 
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Fig.  6.  An  approximate  linear  dependence  observed  between  surface 
dielectric  constant  measurements  and  air  temperature. 
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Fig.  7.  2-D  contour  plot  shows  spatial  variation  of  1 00-M Hz  dielectric 
constant  values  over  an  area  representing  approximately  20  SAR  reso¬ 
lution  cells.  Values  of  10  represent  melt-pond  areas  where  no  measurements 
were  obtained. 
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Fig.  8.  Measured  values  of  both  the  real  part  of  the  dielectric  constant 
and  the  loss  tangent  increase  with  increasing  snow  liquid  water  content. 
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a.  Observed  relationship  between  measurement  values 
(correlation  coefficient  ~  0.6).  The  effect  of  liquid  water 
content  on  the  measured  density  has  been  removed. 
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b.  Relationship  between  dielectric  constant  t  snow  den¬ 
sity  eDS,  and  liquid  water  content  Wv  derived  from  the 
data  plotted  for  Wv  ranging  from  0  to  6%  by  volume. 

Fig.  9.  Dielectric  constant  (real  part )  versus  snow  density. 


regression  of  against  t 1  results  in  a  correlation 
coefficient  of  0.8.  The  regression  coefficient  for 
vs  tan6  is  0.6. 

Snow  temperature,  grain  size,  and  density  varia¬ 
tions  can  also  affect  dielectric  properties,  although 
their  importance  is  greatly  reduced  in  the  summer 
season.  During  the  measurement  period  the  snow 
temperature  remained  near  the  melting  point, 
ranging  between  -1°C  and  0°C.  Such  small  varia¬ 
tion  was  insufficient  to  produce  detectable  varia¬ 
tions  in  the  dielectric  properties. 


The  measured  values  of  the  real  part  of  the  di¬ 
electric  constant  do  show  a  weak  correlation  with 
snow  density  (Figure  9a)  but  no  correlation  with 
grain  size.  Doing  a  multiple  linear  regression  of  li¬ 
quid  water  and  snow  density  against  t  ’  resulted  in 
the  relationship 

f '  =  I  +  2.2eDS  +  2%W^  (1) 

where  gDS  is  an  equivalent  dry  snow  density  in 
which  the  liquid  water  is  replaced  by  air: 


r.tv  sfT 


a.  Observed  relationship  between  measurement  values 
(correlation  coefficient  =  0.4). 


b.  Relationship  between  loss  tangent  LT,  grain  size  d, 
and  liquid  water  content  Wv  derived  from  the  data  plot¬ 
ted  for  W\  ranging  from  0  to  6%  by  volume. 

Fig.  10.  Loss  tangent  versus  grain  size. 
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though  weakly  so,  as  indicated  in  Figure  10a.  A 
multiple  linear  regression  of  liquid  water  content 
and  grain  size,  d,  produced  the  relationship 


This  relationship  between  t '  and  dry  snow  density 
for  a  range  of  liquid  water  content  values  is  plot¬ 
ted  in  Figure  9b. 

Measurements  of  loss  tangent  appear  to  be  cor¬ 
related  not  with  density  but  with  grain  size,  al¬ 


tan5(  x  10’)  =  3.1  +  1.91FV+  \A5d.  (3) 

Following  this  relationship  the  loss  tangent  is  plot¬ 
ted  against  grain  size  in  Figure  10b  for  several  val¬ 
ues  of  liquid  water  content. 
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Fig.  11.  Comparison  of  salinity  (a)  and  100-MHz  dielectric  con¬ 
stant  measurements  (b  and  c)  as  a  function  of  depth  into  the  ice 
for  multiyear  floe  core  11-1  shows  dose  correlation  between  loss 
tangent  and  salinity. 


Ice  cores 

Under  certain  conditions,  such  as  relatively  dry 
and/or  thin  snow  cover  and  the  use  of  relatively 
long  wavelengths,  penetration  through  the  snow 
and  into  the  ice  can  be  significant.  The  dielectric 
properties  of  the  sea  ice  then  become  important. 
Depth  profiles  of  the  real  part  of  the  dielectric 
constant  and  loss  tangent  were  obtained  for  ice 
cores  from  two  multiyear  floes,  one  core  from  a 
relatively  uniform  region  of  congelation  ice  and 
one  from  a  ridge.  These  are  shown  in  Figures  1 1 
and  12  along  with  salinity  profiles  determined 
from  adjacent  cores  by  the  CRREL  team. 


Measurements  were  made  using  the  100-MHz 
Q-meter  device  in  the  ship’s  coldrooms  approxi¬ 
mately  10-12  hours  after  the  cores  were  extracted 
from  the  floes.  Samples  were  cut  from  labeled 
core  segments  that  were  being  stored  for  structural 
analysis  at  -10°C.  Actual  measurements  were  car¬ 
ried  out  in  an  adjoining  coldroom  at  0°C,  which 
may  have  induced  some  surface  melting  in  the 
samples  as  will  be  discussed  later. 

For  both  cores  the  loss  tangent  measurements 
appear  to  be  positively  correlated  with  salinity, 
showing  a  general  increase  with  increasing  salinity 
as  expected  from  previous  measurements  (e.g. 
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Fig.  12.  In  contrast  to  Figure  II,  depth  profiles  for  multiyear 
ridge  core  13-1  show  little  correlation  between  salinity  and  loss 
tangent. 


Vant  et  al.,  1978).  The  loss  tangent  profile  for  core  structural  variations  found  in  this  multiyear  ridge 

11-1  closely  follows  the  form  of  the  salinity  pro-  core.  Because  of  these  uncertainties  and  the  time 

file,  whereas  that  of  core  13-1  shows  significantly  lag  between  the  extraction  of  the  cores  and  these 

more  deviation.  Part  of  this  discrepancy  may  be  measurements,  a  relationship  between  loss  tangent 

attributable  to  mislabeling  of  core  segments.  A  and  salinity  cannot  be  reliably  derived  from  these 

separate  salinity  measurement  made  on  the  ice  data. 

sample  used  for  the  dielectric  constant  measure-  The  measurements  of  the  real  part  of  the  dielec- 

ment  at  "360  cm”  produced  a  salinity  value  of  trie  constant  show  much  less  correlation  with  the 

0.5°/oo,  corresponding  to  the  values  at  270-320  cm  salinity  profiles  than  do  the  loss  tangent  measure- 

in  the  salinity  profile.  The  dip  in  loss  tangent  at  ments.  The  real  part  is  more  sensitive  to  porosity 

"360  cm”  in  the  loss  tangent  profile  may  therefore  than  to  salinity  so  that  some  of  the  variations  ob- 

actually  be  due  to  the  dip  at  270-320  cm  in  the  sa-  served  in  the  dielectric  constant  profiles  probably 

Unity  profile.  Other  deviations  may  be  due  to  the  reflect  structural  changes,  especially  for  core  13-1, 
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which  is  composed  of  several  different  segments 
of  frazil  and  congelation  ice  and  mixtures  of  the 
two  (A.  Gow,  private  communication).  However, 
no  clear  correlation  between  dielectric  constant 
and  structure  is  observed  when  the  two  profiles  are 
compared.  The  extremely  high  values  (>  5.5)  can¬ 
not  be  attributed  to  structural  fluctuations  and 
probably  reflect  the  presence  of  liquid  water  either 
on  the  walls  of  the  100-MHz  cavity  or  in  the  sam¬ 
ple.  This  unavoidably  results  from  the  necessity  of 
making  measurements  in  the  0°C  coldroom.  These 
values  are  therefore  not  representative  of  the  in 
situ  dielectric  properties  of  the  sea  ice  itself  and 
should  be  discounted. 


CONCLUSIONS  AND 
RECOMMENDATIONS 

Analysis  results  from  measurements  of  snow 
cover  properties  made  from  the  Polarstern  during 
MIZEX  84  have  been  presented  here  with  the  aim 
of  providing  relationships  between  dielectric  prop¬ 
erties  and  surface  and  environmental  conditions 
that  can  be  used  in  the  interpretation  of  remotely 
sensed  microwave  data  taken  during  the  experi¬ 
ment.  The  major  results  obtained  from  this  analy¬ 
sis  are: 

1 .  Documentation  of  the  temporal  variation  in 
dielectric  properties  of  snow  cover.  Results  show  a 
linear  increase  in  dielectric  constant  with  increased 
air  temperature,  which  probably  reflects  the 
change  in  snow  liquid  water  content.  These  meas¬ 
urements  need  to  be  compared  with  values  of  sur¬ 
face  radiant  heat  flux  also. 

2.  Determination  of  the  relationship  between  li¬ 
quid  water  content  and  dielectric  properties  for 
snow  cover  in  the  MIZ.  It  was  found  that  the  real 
part  of  the  dielectric  constant  depended  on  liquid 
water  content  and,  to  a  lesser  extent,  snow  densi¬ 
ty.  The  loss  tangent  (imaginary  part/real  part) 
also  depended  on  liquid  water  content  and  was 
weakly  correlated  to  grain  size.  Additional  meas¬ 
urements  are  needed  over  a  wider  range  of  densi¬ 
ties  and  snow  temperature. 

3.  Comparison  of  the  depth  profiles  of  dielectric 
constant  and  salinity  for  multiyear  ice  cores.  This 
showed  that  a  high  correlation  does  exist  between 
loss  tangent  and  salinity,  but  that  other  factors 


(ice  structure)  control  the  behavior  of  the  real  part 
of  the  dielectric  constant. 

In  carrying  out  this  analysis,  several  potential 
improvements  to  the  measurement  routine  have 
become  evident.  One  parameter  that  was  not 
measured  consistently  was  the  air  temperature  at 
the  surface  where  measurements  were  being  made. 
Observations  through  several  diurnal  cycles  of  the 
same  snowpack  are  also  needed  to  monitor  snow 
metamorphosis.  These  data,  combined  with  radi¬ 
ant  heat  flux  measurements,  would  provide  furth¬ 
er  insight  into  the  dependence  of  electrical  surface 
properties  and  meteorological  parameters.  With 
regard  to  sampling  techniques,  liquid  water  meas¬ 
urements  need  to  be  made  in  situ  from  the  same 
sample  level  as  dielectric  constant  measurements, 
and  dielectric  constant  measurements  on  ice  cores 
need  to  be  made  immediately  after  they  are  ex¬ 
tracted,  with  simultaneous  temperature  and  densi¬ 
ty  measurements. 
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1.  INTRODUCTION 

The  East  Greenland  Current  (EGC)  is  a  primary 
oceanographic  circulation  feature  in  the  Green¬ 
land  Sea  (Figure  1).  This  current  flows  southward 
through  Fram  Strait  and  is  initially  a  broad,  strait¬ 
wide  flow  with  a  major  component  entering  south- 
westward  from  the  region  north  of  Svalbard.  Far¬ 
ther  south  this  flow  occupies  only  that  portion  of 
the  Fram  Strait-Greenland  Sea  region  located 
west  of  the  East  Greenland  Polar  Front  (EGPF), 


and  contains  a  pronounced  jet-like  high-speed 
core  just  west  of  this  front.  The  EGPF  separates 
cold,  low-salinity  water  exiting  the  Arctic  Basin 
from  warmer,  more  saline  water  in  the  central 
Greenland  Sea.  In  central  Fram  Strait  part  of  the 
current  becomes  incorporated  into  a  cyclonic  gyre 
which  overlies  the  Molloy  Deep.  By  the  time  it  has 
come  as  far  south  as  about  79°N,  the  jet-like  core 
of  the  current  has  impinged  upon  the  steep  East 
Greenland  continental  slope,  characterized  at  that 
point  by  north-south-trending  isobaths  which  con- 
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Fig.  1.  Geographical  location  of  the  1984-85  MIZEX  current 
moorings  (I  and  2),  schematic  representation  of  the  major  regional 
oceanic  circulation  features  farrows),  location  of  the  East  Green¬ 
land  Polar  Front  (dashed  line)  and  major  bathymetric  features. 
Depths  in  the  Molloy  Deep  exceed  5500  m. 
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strain  the  flow  to  a  north-south  orientation.  The 
current  coincides  approximately  with  the  outer 
marginal  ice  zone  (MIZ)  from  the  vicinity  of  Sval¬ 
bard  southward  to  the  southern  tip  of  Greenland, 
and  it  is  suspected  that  the  current  exerts  a  strong 
influence  over  physical  processes  in  the  MIZ. 

Though  the  EGC  has  long  been  recognized  as  a 
distinct  feature,  it  has  been  only  within  the  past 
decade  or  so  that  sufficient  field  work  has  been 
carried  out  to  allow  a  preliminary  quantitative  de¬ 
scription.  As  part  of  the  ongoing  effort  to  describe 
and  understand  this  current,  and  in  an  attempt  to 
identity  its  interactions  with  the  MIZ,  two  current 
moorings  were  deployed  within  the  EGC  over  a  13- 
month  period  in  1984-85.  The  moorings  were 
located  near  the  core  of  the  current  at  about  the 
1000-m  isobath  along  the  continental  slope,  at 
about  78°30'  and  78°45'N  (see  Figure  1).  These 
locations  lay  west  of  the  ice  edge  and  were  ice- 
covered  throughout  most  of  the  observation  peri¬ 
od.  Each  mooring  contained  two  current  meters 
which  sampled  the  upper  and  lower  layers  of  the 
EGC,  resulting  in  four  records  of  current  speed 
and  direction,  each  about  13  months  long.  This 
report  describes  the  observation  program  and  pre¬ 
sents  and  summarily  discusses  the  most  significant 
features  seen  in  the  resulting  current  data. 

The  data  discussed  in  this  report  were  obtained 
in  conjunction  with,  and  as  part  of,  the  summer 
1985  Marginal  Ice  Zone  Experiment  (MIZEX). 


2.  THE  OBSERVATION  PROGRAM 

The  two  current  moorings  were  deployed  nomi¬ 
nally  along  the  1000-m  isobath  of  the  East  Green¬ 
land  continental  slope  in  mid-June  1984  (Figure  1, 
Table  1).  Variations  in  actual  bottom  depths  of 
about  1000  m  at  the  times  of  deployment  were  due 
to  ship  drift  which  occurred  during  the  deploy¬ 
ment  operations.  Mooring  configuration  was  of 
the  “standard”  oceanographic  taut-wire  type, 


with  36-in.  spherical  steel  floats  positioned  at 
about  90-m  depth  for  flotation.  This  depth  is  well 
below  that  at  which  wave-induced  rotor  pumping 
is  believed  to  have  a  significant  biasing  effect  on 
the  recorded  current  data.  Acoustic  releases  were 
positioned  below  the  deeper  current  meters  on 
each  mooring.  The  deployments  took  place  during 
the  early  part  of  the  summer  1984  MIZEX  field 
program,  and  were  carried  out  from  the  German 
research  icebreaker  Polarstern.  The  moorings 
were  recovered  in  mid- July  1985  from  the  same 
vessel. 

The  instruments  used  were  Aanderaa  RCM-4 
recording  current  meters.  These  units  measure 
speed  to  a  factory-advertised  accuracy  of  about 
±  1  cm/s  (at  the  speeds  which  were  observed),  and 
direction  to  about  ±5°.  The  sampling  intervals 
were  60  minutes,  with  speeds  integrated  over  this 
period  and  direction  measured  once  instantane¬ 
ously  during  each  interval.  The  raw  data  were  first 
inspected  visually  for  bad  data  points.  This  inspec¬ 
tion  revealed  that  the  first  4  days’  upper  level  rec¬ 
ord  at  mooring  2  was  irretrievably  contaminated 
by  “spiking,”  and  this  short  segment  of  the  record 
was  deleted.  Some  minor  spiking  during  the  final 
month’s  upper  level  record  from  mooring  1  was 
edited  out  of  the  data,  leaving  a  satisfactorily 
clean  record.  As  a  final  check  on  data  quality,  har¬ 
monic  tidal  analyses  were  carried  out  on  sequen¬ 
tial  30-day  record  segments.  Invariant  amplitudes, 
phases  and  major  axis  orientation  of  the  tidal  cur¬ 
rent  constituents  indicated  that  the  speed,  time 
base  and  directional  data  are  reliable  and  that  ac¬ 
curacy  falls  within  the  above  specifications. 

To  provide  a  check  on  whether  the  moorings 
“heeled  over”  during  periods  of  strong  currents, 
which  would  have  pulled  the  current  meters  into 
deeper  water  than  under  static  conditions  and  re¬ 
sulted  in  recording  of  current  data  from  depths 
greater  than  planned,  the  uppermost  meters  were 
equipped  with  pressure  gauges.  Some  such  deep¬ 
ening  was  observed  on  the  few  occasions  when 


Table  1.  Particulars  of  the  June  1984-July  1985  MIZEX  current  moorings. 


Mooring 

ID 

Latitude 

(North) 

Longitude 

(West) 

Date 

deployed 

(Z) 

Date 

recovered 

IZ) 

Bottom 

depth 

(m) 

Instrument 

depth 

(m) 

1 

78°43.988 ' 

4°5 1.131  ' 

15  June  84 

16  July  85 

994 

94 

394 

2 

78°29.I60' 

4  ”33.288  ’ 

15  June  84 

18  July  85 

1020 

120 

420 

42 
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strong  (about  30  cm/s)  currents  were  present  at 
both  depths,  but  these  depth  increases  amounted 
to  only  a  few  meters  and  the  incidents  were  not 
frequent  or  severe  enough  to  affect  the  longer- 
term  statistical  validity  of  the  data. 


3.  THE  OBSERVED  CURRENTS 

The  four  current  records  obtained  were  each 
about  13  months  long.  This  length  is  adequate  to 
detect  seasonal  variations  in  current  and,  in  addi¬ 
tion,  mesoscale  variability  and  higher-frequency 
tidal  and  near-inertial  motions  which  can  also  be 
characterized.  Because  these  very  generalized 
classes  of  motion  reflect  different  dynamic  proc¬ 
esses,  the  following  section  is  subdivided  based 
upon  the  time  scales  of  the  observed  currents  into 
low-frequency  (seasonal),  mesoscale,  tidal,  and 
near-inertial. 

3.1  Low-frequency  currents 

Each  of  the  four  current  records  was  first  re¬ 
solved  into  its  u  (cross-isobath,  east  positive)  and  v 
(along-isobath,  north-positive)  components;  u  and 
v  were  taken  to  be  in  the  same  direction  at  both 
moorings  so  as  to  simplify  comparisons  between 
the  records  and  also  because  the  isobaths  at  each 
mooring  trended  within  a  few  degrees  (i.e.  within 
the  directional  accuracy  of  the  current  meters)  of 
true  north.  The  resulting  hourly  u  and  v  compo¬ 
nents  were  then  subjected  to  a  35-hour  low-pass 
filter  to  remove  tidal  and  inertial  constituents.  The 
low-pass  filtered  records  were  then  passed  through 
a  30-day  “boxcar”  filter  which  effectively  com¬ 
putes  30-day  mean  values  with  a  15-day  overlap 
between  consecutive  30-day  intervals.  A  95 Vo  er¬ 
ror  estimate  was  computed  for  each  interval, 
where  error  =  (2  standard  deviations)/(square  root 
[record  length/2  x  integral  time  scale])  and  the  in¬ 
tegral  time  scale  is  defined  as  the  integral  from 
zero  to  infinity  under  the  autocorrelation  curve 
(Lagerloef  et  al.,  1981).  The  results  of  these  com¬ 
putations  are  presented  in  Figure  2,  and  the 
discussion  below  refers  to  this  figure. 

The  30-day  mean  v  components  (Figure  2a) 
were  almost  always  southward,  the  only  exception 
being  a  brief  period  of  northward  flow,  not  signif¬ 
icant  at  the  95  Vo  level,  which  occurred  at  mooring 
1  in  June  1985.  The  southward  30-day  mean 
speeds  were  greatest  in  summer  and  for  the  upper 
records,  and  were  smallest  in  winter  for  both  the 
upper  and  lower  records.  Maximum  southward 
summer  speeds  of  order  10  cm/s  occurred  in  the 
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a.  Along-shelf  (v)  components. 
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b.  Cross-shelf  (\i)  components. 

Fig.  2.  Very  low  pass  (30-day)  filtered  current  com¬ 
ponents.  Vertical  bars  indicate  error  estimates  computed 
as  described  in  the  text.  Last  month  of  series  (8  June-8  July 
1985)  is  omitted  for  convenience  in  plotting. 


upper  record  at  mooring  1,  and  speeds  of  more 
than  15  cm/s  were  observed  in  the  upper  record  at 
mooring  2.  Summer  lower  level  speeds  were  of 
order  5  and  10  cm/s  at  moorings  1  and  2  respec¬ 
tively,  so  that  the  speed  differences  between  upper 
and  lower  levels  were  about  5  cm/s  at  both  moor¬ 
ings.  Speeds  were  lower  in  winter  than  in  summer, 
being  less  than  5  cm/s  at  mooring  1  and  of  order  5 
cm/s  at  mooring  2.  The  differences  between  upper 
and  lower  level  speeds  were  generally  less  than  2 
cm/s  in  winter  and  were  not  significant  at  the  95% 
level.  The  large  error  bars  during  winter  were  due 
to  occurrence  of  energetic  mesoscale  disturbances, 
noted  in  Section  3.2,  which  showed  strong  separa¬ 
tion  between  upper  and  lower  layers  in  contrast  to 
the  weak  vertical  shear  in  the  30-day  means. 

The  30-day  mean  u  components  (Figure  2b) 
were  smaller  than  v  at  both  moorings  and  at  both 
levels;  u  was  always  less  than  about  3  cm/s  and 
was  often  not  significantly  (at  the  95%  level)  non¬ 
zero.  There  was  no  apparent  preferred  long-term 
direction  in  u\  however,  there  was  a  tendency  for 
weak  (1-2  cm/s)  eastward  flow  in  summer  1984  at 
both  depths  at  mooring  1  and  for  westward  flow 
in  summer  1985  at  both  moorings.  In  winter,  u  for 
both  levels  at  the  two  moorings  was  less  than  1 
cm/s  and  was  often  negligibly  small. 

The  above  describes  a  mean  (30-day  averaged) 
flow  which  was  consistently  along-isobath  and 
southward  at  speeds  of  5-15  cm/s  (upper  level) 
and  2-10  cm/s  (lower  level),  and  which  had  much 
smaller  (less  than  3  cm/s)  cross-isobath  speeds 
with  no  apparent  preferred  direction.  The  higher 
speeds  occurred  in  summer  1984.  Vertical  speed 
differences  between  the  two  levels  were  greatest  at 
both  mooring  locations  in  summer  1984,  with  up¬ 
per  level  speeds  typically  5  cm/s  greater  than  those 
in  the  lower  level.  In  winter,  vertical  speed  differ¬ 
ences  were  typically  2  cm/s  or  less. 

There  are  few  available  observations  with  which 
to  compare  these  mean  currents.  Paquette  et  al. 
(1985)  provide  a  useful  summary  of  past  work  on 
the  EGPF  and  present  new  data  from  October- 
November  1981.  They  have  constructed  dynamic 
topographies  which  show  a  southward  net  baro- 
clinic  flow,  both  at  the  surface  and  at  150  db,  at 
the  mooring  locations.  Their  data  were,  however, 
insufficient  to  define  the  net  flow  orientation. 
They  cite  frontal  jet  speeds  approaching  100  cm/s, 
far  highe'  than  observed  at  any  time  in  our  data, 
and  it  seems  apparent  that  moorings  1  and  2  were 
not  situated  correctly  to  sample  the  highest  speed 
front-associated  currents  (“jet”).  The  100-m 
mean  currents  were  10-15  cm/s  in  autumn  1984, 


and  the  400-m  currents  were  still  between  5  and  10 
cm/s.  Paquette  et  al.  used  the  500-db  level  as  a  ref¬ 
erence  in  computing  dynamic  heights,  and  our  ob¬ 
servations  show  that  there  was  still  considerable 
motion  at  400  m  in  autumn  1984.  A  more  quanti¬ 
tative  representation  of  the  upper  layer  dynamic 
topography  would  therefore  require  extrapolation 
of  a  deeper  reference  level  onto  the  continental 
shelf  from  deeper  water.  Paquette  et  al.  observe 
that  surface  current  speeds  are  typically  3  to  10 
times  higher  than  at  150  db,  consistent  with  baro- 
clinic  shear  computed  from  historical  data.  Ap¬ 
plied  to  our  current  observations,  the  overlying 
surface  current  speeds  could  have  been  as  high  as 
about  150  cm/s,  perhaps  an  unrealistically  high 
value  but  not  incompatible  with  the  abovemen- 
tioned  100-cm  7s  currents  associated  with  the  fron¬ 
tal  jet. 

The  presence  of  greater  mean  current  speeds  at 
mooring  2  than  at  mooring  1  can  be  qualitatively 
accounted  for  by  referring  to  dynamic  topogra¬ 
phies  presented  by  Paquette  et  al.,  and  also  by  re¬ 
ferring  to  Figure  1 .  Atthe  approximate  locations  of 
the  moorings,  the  EGC  is  impinging  from  the 
northeast  onto  the  Greenland  continental  slope,  so 
that  the  moorings  were  located  in  a  zone  of  con¬ 
vergence  between  the  slope-constrained  southward 
flow  and  the  main  core  of  the  EGC.  Higher  speeds 
at  mooring  2  probably  reflect  a  greater  transport, 
due  to  this  convergence,  at  that  mooring  than  at 
mooring  1. 

3.2  Mesoscale  fluctuating  currents 

Mesoscale  is  taken  here  to  include  those  fluctua¬ 
tions  which  occur  primarily  within  the  2-  to  12-day 
period  band  and  exclude  shorter  period  tidal  and 
inertial  motions.  Mesoscale  features  have  been  ob¬ 
served  throughout  the  world  ocean  to  be  typical  of 
fronts,  meanders,  eddies  and  wind-driven  motions. 

The  current  records  were  resolved  into  u  and  v 
components  and  run  throuigh  a  35-hour  low-pass 
filter  as  described  in  Section  3.1.  The  qualitative 
nature  of  the  mesoscale  fluctuations  can  be  seen  in 
the  resulting  time-series  of  u  and  v  components 
(Figure  3).  For  the  early  part  of  the  records,  from 
their  start  until  about  late  October,  the  greatest 
fluctuations  were  typically  about  15  cm/s  in  am¬ 
plitude.  A  notable  exception  was  the  large  (about 
30  cm/s)  single  oscillation  which  occurred  in  late 
October,  primarily  at  mooring  I,  though  it  was 
present  also  in  the  lower  level  at  mooring  2.  From 
late  October  until  early  February  all  records  were 
"quieter”  than  earlier.  Through  the  mid-winter 
February-April  period  the  currents  were  charac- 
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Fig.  3.  Low  pass  (35-hour)  filtered  u  and  v  current  components. 


terized  by  consistently  large  (20-30  cm/s),  persis¬ 
tent  fluctuations  having  shorter  time  scales  than 
those  observed  through  summer  1984.  Finally,  by 
May-June  1985  the  fluctuations  had  returned  to  a 
pattern  similar  to  that  observed  in  summer  1984. 

The  above  observations  were  quantified  by 
computing  energy-preserving  spectra  for  the  u  and 


v  components  from  all  four  time  series  (Figure  4). 
These  spectra  were  computed  for  sequential  90- 
day-Iong  records,  with  a  30-day  overlap  at  each 
end.  Each  spectrum  was  band-averaged  to  give  es¬ 
timates  with  10  degrees  of  freedom.  This  method 
provides  a  sufficiently  long  record  for  each  spec¬ 
tral  computation  to  give  confidence  in  the  results. 
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and  also  breaks  the  entire  record  into  enough  seg¬ 
ments  to  identify  the  long-term  differences  which 
were  noted  visually  from  the  time-series  plots. 
These  spectra  yield  the  following  information. 

1.  The  significant  mesoscale  energy  was  con¬ 
tained  in  the  period  band  between  about  2 
and  12  days. 

2.  The  upper  level  observations  at  both  moor¬ 
ings  showed  a  tendency  for  energy  to  be  con¬ 
centrated  at  longer  periods  (8-12  days)  in 
summer  and  shorter  periods  (2-4  days)  in 
winter. 

3.  No  significant  changes  in  spectral  energy  or 
shape  occurred  during  June-November 
1984,  except  for  some  increase  in  longer  peri¬ 
od  energy  in  the  upper  layer  at  mooring  1  in 
August-November,  in  qualitative  agreement 
with  our  above  observation  that  the  time- 
series  appeared  steady  through  this  period. 

4.  October  1984-January  1985  was  character¬ 
ized  by  low  energy  levels,  in  contrast  to  the 
preceding  and  following  periods. 

5.  December-March  was  characterized  by  a 
dramatic  increase  in  upper  layer  energy 
levels.  This  affected  primarily  the  u  compo¬ 
nent  at  mooring  1,  but  affected  u  and  v  near¬ 
ly  equally  at  mooring  2.  This  energy  was 
dominated  by  a  very  sharp  peak  at  about  a 
4-day  period,  which  reflects  presence  of  the 
strong  oscillations  which  were  noted  in  the 
time-series  plots  for  February-March  1985. 
This  energy  level  persisted  into  the  late  win¬ 
ter  (February-May)  period  in  both  the  upper 
and  lower  layer  v  components  at  mooring  2 
but  not  at  mooring  1 .  The  significant  vertical 
difference  in  mesoscale  energy  during  this 
period  is  in  striking  contrast  to  the  30-day 
mean  flow  (Figure  2),  which  implied  mini¬ 
mum  vertical  shear  during  the  same  period. 

6.  By  spring-summer  (April-July)  1985  the 
spectral  signatures  appeared  to  be  reverting 
to  those  observed  during  the  previous  sum¬ 
mer. 

Finally,  empirical  orthogonal  functions  (EOF’s) 
were  computed  for  four  different  periods,  each  90 
days  in  length  and  chosen  based  upon  inspection 
of  both  the  time-series  and  the  spectra  to  encom¬ 
pass  periods  characterized  by  quasi-steady  behav¬ 
ior.  These  periods  were: 

l.  Ninety  days  starting  15  June  1984,  falling  in 
the  interval  of  strongest  vertical  mean  shear 
(as  inferred  from  the  differences  between  up¬ 
per  and  lower  level  currents)  and  longer  peri¬ 
od  mesoscale  oscillations. 


11.  Ninety  days  starting  15  October  1984,  falling 
in  the  interval  of  intermediate  vertical  mean 
shear  and  low  mesoscale  variance. 

III.  Ninety  days  starting  15  January,  encompass¬ 
ing  the  period  of  weakest  mean  currents, 
smallest  vertical  mean  shear  and  large,  rela¬ 
tively  short  period  mesoscale  oscillations. 

IV.  Ninety  days  starting  15  April,  falling  in  an 
interval  of  intermediate  current  speeds,  ver¬ 
tical  mean  shear  and  mesoscale  variance 
energy. 

Each  mooring  was  analyzed  separately.  With 
two  sample  depths  each,  the  records  were  decom¬ 
posed  into  u  and  v  components  which  generated 
four  input  time  series  for  the  EOF  analyses.  These 
analyses  produced,  in  turn,  four  EOF  modes  de¬ 
scribing  the  vertical  empirical  structure  at  each  lo¬ 
cation  for  each  period.  The  resulting  eigenvectors 
representing  the  u  and  v  components  for  each  rec¬ 
ord  are  recombined  into  a  vector  representation 
(Figure  5).  (We  note  that  use  of  complex  EOF  an¬ 
alyses  would  have  produced  two  rather  than  four 
modes  and  would  only  have  shown  relative  as 
compared  to  absolute  orientation  of  the  eigenvec¬ 
tors.) 

In  all  cases  more  than  80%  of  the  variance  was 
contained  in  the  first  two  modes;  therefore,  only 
these  two  modes  are  discussed  here.  Forming  dy¬ 
namical  interpretations  of  empirically  derived 
modes  requires  that  decoupled  dynamic  modes 
dominated  the  system  throughout  the  period 
under  analysis.  For  example,  definition  of  baro- 
tropic  and  baroclinic  modes  requires  that  the 
former  have  nearly  equal  eigenvectors  at  each 
depth.  This  situation  is  clearly  not  prevalent  in  the 
results  depicted  in  Figure  5.  Consequently,  the  fol¬ 
lowing  discussions  emphasize  the  statistical  prop¬ 
erties  which  can  be  described  from  EOF  analyses. 

Mode  1,  the  dominant  mode,  contained  be¬ 
tween  46%  and  65%  of  the  total  mesoscale  vari¬ 
ance  energy.  This  mode  rotated,  at  both  moorings 
and  depths,  from  a  primarily  cross-shelf  direction 
in  June-September  to  primarily  along-shelf  by  the 
following  spring.  The  strong  winter  cross-shelf 
orientation  at  mooring  1  (also  evident  as  a  cross¬ 
shelf  mode  2  at  mooring  2)  was  an  exception,  and 
was  due  to  the  strong  February-March  east-west 
fluctuations  which  were  noted  above.  The  second 
mode  contained  19-41%  of  the  total  variance,  so 
that  for  any  specified  period  more  than  80%  of 
the  variance  was  contained  in  modes  1  and  2.  The 
second  mode  behaved  in  inverse  fashion  to  the 
first,  i.e.  it  rotated  from  primarily  along-shelf  in 
June-September  to  primarily  cross-shelf  by  the 
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Fig.  5.  EOFs  computed  for  four  separate  periods.  Upper  plots  are  for 
first  mode  EOFs,  while  lower  plots  are  for  second  mode.  Dashed  lines  show 
local  bathymetric  trend.  Solid  dots  indicate  mooring  locations  (uppermost  is 
/  and  to  wer  is  2  for  each  EOF  plot),  and  nearby  numbers  indicate  percentage 
variance  contained  in  that  mode.  Smaller  numbers  near  "heads”  of  vectors 
indicate  depths. 


following  spring.  In  June-September  46%  of  the 
variance  was  cross-shelf  (mode  1)  while  36-41% 
was  along-shelf  (mode  2);  the  difference  between 
these  two  modes  is  probably  not  significant  for 
this  period,  which  we  would  expect  because  u  and 
v  fluctuations  were  similar  during  this  period  (Fig¬ 
ure  3).  In  October-January  the  modes  were  both 
inclined  to  the  isobaths.  This  period,  however, 
was  one  of  quite  low  variance  (see  the  Octo¬ 
ber-January  spectra  in  Figure  4)  so  that  the  signif¬ 
icance  of  the  modes  for  this  period  is  uncertain.  In 
January-April  the  upper  layer  at  mooring  1  shows 
the  influence  of  the  vigorous  cross-shelf  fluctua¬ 
tions,  an  influence  which  was  evident  in  the  weak¬ 
er  cross-shelf  mode  2  at  mooring  2.  By  April- 
July  the  majority  of  the  variance  was  still  along- 
shelf,  65%  at  mooring  1  (mode  1)  and  59%  at 
mooring  2  (mode  2). 

The  EOFs  also  provide  information  on  vertical 
differences  in  the  variance.  In  June-September  the 
lower  level  EOFs  were  typically  about  half  those  in 
the  upper  level.  In  Octooer-January  the  upper  and 
lower  level  values  were  more  nearly  equal.  In  Jan¬ 
uary-April  the  picture  was  more  confused,  with 
nearly  all  of  the  variance  at  mooring  1  occurring  in 
the  upper  layer  whereas  mooring  2  had  significant 
mode  1  energy  in  the  lower  layer.  By  April-June 
the  ratio  between  upper  and  lower  layer  variance 
appeared  to  be  approaching  that  which  had  char¬ 
acterized  the  previous  summer,  though  there  was 
proportionally  more  lower  level  variance  than  had 
been  present. 

To  summarize,  the  mesoscale  variance  energy  in 
summer  1984  fell  primarily  in  the  8-  to  12-day  per¬ 
iod  band,  was  about  twice  as  intense  in  the  upper 
as  in  the  lower  level,  and  was  about  equally  direct¬ 
ed  across-  and  along-shelf.  Autumn  1984  saw  the 
lowest  observed  variances,  with  distribution 
across  the  frequency  band  and  poor  directional 
definition.  In  winter  1985  variance  was  high,  was 
confined  primarily  to  the  2-  to  6-day  period  band, 
was  primarily  in  the  upper  level,  and  was  both 
cross-  and  along-shelf.  In  spring  1985  the  variance 
had  intermediate  values,  was  broad  band  and  was 
primarily  along-shelf. 

The  EGC/EGPF  system  has  been  characterized 
as  having  a  rich  mesoscale  eddy  field  (Johannessen 
et  al.,  1983).  Time-series  current  observations 
which  coincided  roughly  with  the  summer  1984 
MIZEX  field  experiment  indicate,  in  fact,  propa¬ 
gation  of  eddies  or  eddy-like  features  past  moor¬ 
ings  1  and  2  (Figure  6).  The  most  obvious  such 
feature,  occurring  on  about  28  June-2  July,  was 
readily  identifiable  as  an  anticyclonic  feature  hav- 
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Fig.  6.  Sample  segment,  overlapping  in  time  much 
of  the  summer  1984  MIZEX  field  experiment,  of 
low-pass  ( 35-hour )  filtered  vector  current  plots 
showing  signatures  of  eddies. 


ing  associated  speeds  of  10-15  cm/s  in  the  upper 
level  at  mooring  1.  This  feature  was  less  apparent 
at  mooring  2,  but  this  may  have  been  due  in  part 
to  masking  by  the  stronger  mean  southward  flow 
at  the  latter  location  than  at  mooring  1 .  A  less  pro¬ 
nounced  anticyclonic  eddy-like  feature  was  appar¬ 
ent  at  mooring  2  on  about  23-24  June,  but  was 
much  less  apparent  at  mooring  1.  Other  fluctua¬ 
tions  which  occurred  throughout  the  time-series 
segment  can  be  interpreted  variously,  if  not  unam¬ 
biguously,  as  combinations  of  anticyclonic  and  cy¬ 
clonic  eddies.  This  pattern  persisted  throughout 
the  13-month-long  current  records,  though  eddy 
activity  was  minimal  during  the  autumn  period  of 
low  spectral  energy.  In  general,  the  current  obser¬ 
vations  substantiate  prior  claims  that  the  region  is 
rich  in  eddy  activity. 

Finally,  the  highly  energetic  east-west  oscilla¬ 
tions  which  occurred,  with  about  a  4-day  period, 
during  February  1985  at  mooring  1  deserve  com¬ 
ment  because  of  their  high  energy  and  uniqueness 
within  these  records  (no  similar  features  were  ob¬ 
served  at  any  other  time  during  the  mooring  peri¬ 
od).  Current  amplitudes  approaching  30  cm/s 
were  associated  with  these  fluctuations  (Figure  3), 
and  the  resulting  variance  dominated  both  the 
spectra  (Figure  4a)  and  the  mode  1  EOFs  (Figure 
5).  They  were  evident  also  at  mooring  2,  where 
much  of  the  variance  had  been  shifted  to  lower 


frequencies.  They  were  highly  baroclinic,  being 
confined  to  the  upper  layer.  Their  origin  is  un¬ 
known. 
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3.3  Tidal  currents 

The  raw  (i.e.  unfiltered)  current  records  were 
analyzed  for  tidal  current  constituents.  These 
analyses  were  carried  out  on  a  12-month-long  seg¬ 
ment  for  each  record  using  the  method  of  Fore¬ 
man  (1978).  The  tidal  currents  were  dominated  by 
the  semidiurnal  M2  and  S2  constituents.  Of  the 
smaller  diurnal  constituents,  only  K1  was  signifi¬ 
cant.  These  three  constituents  are  graphically 
shown  in  Figure  7.  It  is  apparent  from  these  plots 
that  the  semidiurnal  tidal  currents  were  greater  in 
the  upper  than  in  the  lower  layers.  The  diurnal 
constituents  were,  conversely,  slightly  larger  in  the 
lower  than  in  the  upper  layer.  The  results  from  the 
sequential  monthly  tidal  analyses  which  were  car¬ 
ried  out  to  check  the  current  data  (see  Section  2) 
indicated  that  these  variations  with  depth  were 
nearly  constant  throughout  the  observation 
period. 

A  short  time-series  plot  of  the  non-filtered  cur¬ 
rents  from  mooring  1  indicates  qualitatively  the 
relative  significance  of  tidal  currents  with  respect 
to  the  mean  and  fluctuating  flows  (Figure  8). 
While  the  individual  constituents,  of  which  M2  is 
the  largest  with  a  magnitude  of  about  5  cm/s,  are 
of  the  same  order  or  less  than  the  mean  flow,  the 
interactions  between  constituents  can  create  large 
enough  tidal  period  fluctuations  to  be  significant 
when  viewed  within  the  context  of  short-term  field 
experiments. 

3.4  Near-inertial  currents 

The  unfiltered  current  time-series  plots  showed 
well-defined  semidiurnal  current  fluctuations  (Fig¬ 
ure  8).  As  noted  above,  much  of  this  semidiurnal 
variability  was  due  to  the  M2  and  S2  tidal  consti¬ 
tuents.  However,  the  observed  day-to-day  fluctua¬ 
tions  in  diurnal  current  speeds  were  too  great  to  be 
explainable  in  terms  of  the  tides.  For  example, 
sharp  increases  in  the  semidiurnal  u  and  v  currents 
occurred  on  about  29  November  and  3  December. 
It  is  hypothesized  that  these  fluctuations  were  due 
to  near-inertial  oscillations,  which  at  the  latitudes 
of  the  moorings  should  have  a  limiting  (i.e.  iner¬ 
tial)  period  of  about  12.24  hours. 

To  test  for  the  presence  of  near-inertial  mo¬ 
tions,  the  current  records  were  detided  and  de¬ 
meaned  by  subtracting  the  predicted  tide  and  long¬ 
term  mean  from  the  unfiltered  series.  A  sample 
segment  of  the  unfiltered  data,  the  corresponding 
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Fig.  7.  Plots  of  tidal  current  ellipses  (major  and 
minor  axes)  for  the  three  major  observed  consti¬ 
tuents  M2,  S2andKl.  Upper  plots  are  for  mooring  l, 
and  lower  are  for  mooring  2.  Solid  plots  are  for  nominal 
100-m  depth,  and  dashed  plots  are  for  nominal  400-m 
depth  (see  Table  l  for  actual  depths). 
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Fig.  8.  Time  plots  of  unfiltered  u  and  v  current 
components  showing  raw  data  (upper),  artificially 
generated  tidal  current  series  (middle),  and  detid¬ 
ed  series  (lower). 
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a.  Upper  layer;  the  near-inertial  bands  have  been 
shaded  for  comparison  purposes. 
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b.  Lower  layer. 

Fig.  9.  Spectra  showing  variance  in  the  vicinity  of 
the  semidiurnal  band  before  (left)  and  after  (right) 
detiding  the  time-series  by  subtracting  the  tidal 
series  from  the  raw  data. 


predicted  tidal  current  and  the  residual  currents  is 
shown  in  Figure  8.  Prominent  semidiurnal  current 
oscillations  remain  after  subtraction  of  the  tidal 
currents.  Spectral  analyses  were  then  carried  out 
to  determine  the  nature  of  these  fluctuations.  The 
spectra  were  generated  by  ensemble  averaging  two 
4096-point  FFT’s  with  prewhitening  and  a  Han¬ 
ning  window  with  a  10%  cosine  overlap.  Two  ex¬ 
amples  of  the  spectra  for  the  unfiltered  original 
data  and  for  the  residual  currents  are  presented  in 
Figure  9.  The  spectra  have  been  plotted  in  vari¬ 
ance-preserving  format  (linear  y-axis)  with  a  linear 
frequency  scale  to  display  detail  in  the  vicinity  of 
the  semidiurnal  tidal  band.  These  spectra  demon¬ 
strate  that  most  of  the  tidal  variance  has  been  ef¬ 
fectively  removed;  this  is  especially  evident  for  the 
deeper  of  the  two  records,  where  no  significant 
semidiurnal  variance  remains  in  the  residuals.  The 
upper  level  residual  indicates  a  significant  (at  the 
90%  level)  variance  peak  at  the  true  inertial  fre¬ 
quency,  and  also  two  significant  superinertial 
peaks.  Nothing  can  be  said  concerning  the  relative 
magnitudes  of  the  peaks  because  of  the  broad  con¬ 
fidence  limits;  however,  it  is  apparent  that  signifi¬ 
cant  near-inertial  variance  energy  was  present  in 
the  upper  layer  currents  (Figure  9a)  but  not  in  the 
lower  layer  (Figure  9b). 

The  presence  of  near-inertial  motions  having 
been  established,  certain  features  shown  in  the  re¬ 
sidual  time-series  can  be  explained  (Figure  8).  The 
approximate  equality  of  the  semidiurnal  u  and  v 
components,  coupled  with  an  observed  90-degree 
phase  lead  of  v  ahead  of  u,  are  consistent  with  an- 
ticyclonic,  circular  particle  orbits  which  typify 
near-inertial  motions.  The  large  semidiurnal  cur¬ 
rents  which  were  present  on  about  29  November 
and  3  December  appear  to  be  inertial  “ringing”  in 
response  to  some  strong  forcing— probably  a 
storm  event.  During  these  events  the  near-inertial 
motions  clearly  dominated  the  tidal  currents.  Such 
near-inertial  “ringing”  has  been  observed  at  sites 
throughout  the  world  ocean  (see  for  instance  the 
current  time-series  shown  by  Anderson  et  al., 
1983),  and  is  generally  believed  to  be  generated  by 
surface  wind  events  (see  for  example  d’Asaro, 
1985).  The  more  or  less  continual  lower  amplitude 
near-inertial  “background”  motions  are  of  uncer¬ 
tain  origin.  They  may  be  wind-generated  and 
propagate  into  the  region  from  elsewhere,  or  may 
be  generated  by  such  theoretically  feasible  means 
as  interactions  between  surface  or  internal  waves 
(Hasselmann,  1970).  The  frequency  shifting  to  su¬ 
perinertial  bands  is  likewise  of  uncertain  origin 
and  probably  due  to  a  combination  of  factors  in¬ 
cluding  mean  current  shear  and  lateral  density  gra- 


dients,  as,  for  instance,  reported  along  the  Oregon 
coast  by  Kundu  (1976). 

Virtually  no  near-inertial  motions  were  present 
in  the  lower  layer  at  either  mooring;  these  motions 
were  confined  to  the  upper  layers.  Concentration 
of  these  motions  in  the  upper  layers  is  consistent 
with  upper  layer  generation  by  winds  or  wave  in¬ 
teractions.  Within  this  context,  it  is  also  worth 
noting  that  the  observations  reported  here  were 
from  a  nominal  depth  of  100  m.  Near-inertial  mo¬ 
tions  in  shallower  layers  or  at  the  surface  were 
most  certainly  present  and  were  probably  consid¬ 
erably  more  energetic  than  those  described  here, 
particularly  if  generated  by  surface  forcing  such  as 
storm  winds. 

If  there  are  uncertainties  regarding  the  modes  of 
generation  of  these  near-inertial  motions,  the  un¬ 
certainties  concerning  mechanisms  for  decay  are 
still  greater.  One  possible  theoretical  mechanism 
is,  however,  of  special  interest.  Bell  (1978)  has  hy¬ 
pothesized  that  near-inertial  motions  can  decay 
through  transfer  of  energy  to  the  internal  wave 
field.  Since  internal  waves  are  of  special  signifi¬ 
cance  in  the  marginal  ice  zones  for  a  number  of 
reasons  (see  for  example  Muench  et  al.,  1983),  fur¬ 
ther  studies  should  perhaps  focus  in  part  on  this 
possible  mechanism. 


SUMMARY 

Our  preliminary  analysis  of  the  overwinter 
1984-85  current  observations  from  the  East 
Greenland  Current/Polar  Front  system  has  led  to 
the  following  conclusions: 

1.  Net  flow  was  southward,  along-isobath,  for 
nearly  the  entire  13-month  observation  peri¬ 
od.  Cross-isobath  flow  was  usually  negligi¬ 
ble. 

2.  Southward  current  speeds  (30-day  averaged) 
were  greatest,  and  vertical  speed  differences 
(interpreted  loosely  as  shear)  were  also  great¬ 
est,  during  summer-autumn  and  were  least 
in  mid-winter.  Summer-autumn  upper  layer 
( -  100  m)  speeds  were  10-15  cm/s,  and  low¬ 
er  layer  ( -400  m)  speeds  were  5-10  cm/s.  In 
winter  the  respective  speeds  were  5-10  and 
less  than  5  cm/s. 

3.  Energetic  current  fluctuations  occurred  in 
the  2-  to  12-day  period  band.  In  summer 
these  fluctuations  were  present  in  both  upper 
and  lower  layers  and  were  predominantly 
cross-isobath,  whereas  in  winter  they  were 
confined,  with  some  exceptions,  to  the  upper 


layer,  with  episodes  of  both  along-  and 
cross-isobath  dominance.  In  summer- 
autumn  the  fluctuations  tended  toward  long¬ 
er  (8-  to  12-day)  periods,  whereas  in  winter 
they  were  of  shorter  (4-day)  periods. 

4.  Mesoscale  eddies  were  identifiable  features 
at  the  moorings  throughout  the  observation 
period,  and  were  responsible  for  much  of  the 
observed  variance.  They  were  predominantly 
anticyclonic,  and  were  distributed  similarly 
to  the  variance  as  to  season  and  depth. 

5.  The  tidal  currents  were  dominated  by  the  M2 
(lunar  semidiurnal)  constituent,  which  had  a 
magnitude  of  about  5  cm/s.  The  S2  and  K1 
constituents  were  also  significant.  The  semi¬ 
diurnal  (M2  and  S2)  constituents  were  signif¬ 
icantly  baroclinic  and  decreased  with  in¬ 
creasing  depth.  The  diurnal  (Kl)  constituent 
increased  slightly  with  increasing  depth. 

6.  Near-inertial  current  oscillations  were  pres¬ 
ent  throughout  the  upper  (-100  m)  layer 
records.  Typical  speeds  were  5  cm/s,  but  os¬ 
cillation  amplitudes  as  great  as  about  15 
cm/s  were  observed  on  occasion.  No  signifi¬ 
cant  near-inertial  motions  were  present  in 
the  lower  ( -400  m)  layer. 
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The  aim  of  this  study  is  to  investigate,  by  means 
of  a  coupled  ice/air  model,  the  generation  of 
strong  along-ice-edge  winds  in  the  MIZ,  the  for¬ 
mation  of  the  ice  edge  jet,  and  the  instability  cri¬ 
teria  of  ice  drift.  Ice-drift  observations  in  the 
Greenland  Sea  from  28  April  to  3  September  1978 
(Figure  1)  show  two  different  types  of  ice  motion. 
Oscillation  occurs  mostly  in  spring  (April,  May, 
and  June)  and  nonoscillation  takes  place  in  sum¬ 
mer  (July  and  August).  Ice  melts  and  becomes 
thinner  from  winter  to  summer.  Although  ocean 
eddies  probably  produce  oscillations  in  drift  track, 
another  relation  between  the  pattern  of  ice  motion 
and  ice  thickness  may  still  exist.  One  purpose  of 
the  present  work  is  to  find  this  relation. 


A  possible  mechanism  for  the  air/ice  interaction 
is  presented  in  Figure  2.  The  low-level  air  flow 
generated  by  differential  surface  heating  is  water- 
ward  (icebreeze)  and  equatorward  along  the  ice 
edge  in  both  eastern  Greenland  and  the  eastern 
Antarctic  Peninsula  and  influences  the  ice  through 
surface  air  stress.  However,  movement  of  the  ice 
edge  in  the  MIZ  toward  warm  water  changes  the 
thermal  conditions  near  the  surface  and  produces 
an  air  surface  temperature  gradient  across  the  ice 
edge. 

The  ice/air  interaction  model  depicted  in  this 
article  is  intended  to  simulate  the  main  physical 
processes  and  to  determine  some  instability  cri¬ 
teria  for  the  prediction  of  MIZ  migration. 


Fig.  I .  Ice-drift  observations  in  the  Greenland  Sea. 
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where  8 ^  is  the  basic  air  potential  temperature  at 
the  surface  and  N  is  the  Brunt-Vaisaila  frequency. 

Waterward  migration  of  the  MIZ  increases  the 
surface  temperature  gradient,  and  iceward  migra¬ 
tion  of  the  MIZ  decreases  the  surface  temperature 
gradient.  Surface  temperature  perturbation  is  par¬ 
ameterized  as 

9'.  =  -DT<,[\  +f(0]cos™,  (2) 


Fig.  2.  Ice/air/water  coupled  system. 


Thermally  forced  boundary  layer  air  flow 

The  surface  isotherms,  which  are  nearly  parallel 
to  the  eastern  coast  of  Greenland  over  the  Green¬ 
land  Sea  (Prik,  1959),  indicate  that  air  tempera¬ 
ture  monotonically  increases  waterward  near  the 
MIZ.  Such  differental  surface  heating  generates  a 
local  air  flow  near  the  MIZ.  In  this  section  we 
utilize  a  planetary  boundary  layer  model  treated 
by  Kuo  (1973)  and  Chu  (1985)  to  simulate  a  ther¬ 
mally  forced  boundary  layer  air  flow.  The  coor¬ 
dinate  system  is  chosen  as  moving  with  the  edge  of 
the  MIZ.  The  *-axis  is  in  the  cross-edge  direction, 
and  the  .v-axis  parallels  the  ice  edge,  as  shown  in 
Figure  3.  The  jr-coordinate’s  unit  length  is  twice 
the  MIZ  width,  L  (200  km),  and  that  of  the  verti¬ 
cal  coordinate  is  12  =  (iVO)'1,  where  v  is  the  verti¬ 
cal  eddy  viscosity  and  0  is  the  angular  velocity  of 
the  earth’s  rotation.  The  line  x  =  Vi  is  located  at 
the  ice  edge.  The  MIZ  covers  the  zone  (0  <  x  < 
0.5,  y)  where  the  .y-axis  (*  =  0)  is  near  the  boun¬ 
dary  between  the  MIZ  and  the  interior  ice  pack.  It 
is  considered  that  spatial  variations  in  the  MIZ  are 
much  larger  perpendicular  to  the  ice  edge  than 
parallel  to  it,  so  derivatives  with  respect  to  y  are 
assumed  to  be  zero.  Since  acceleration  of  the  ice 
edge  is  small  compared  to  that  of  the  air  flow,  we 
may  ignore  the  inertial  force  due  to  the  use  of  a  co¬ 
ordinate  system  (moving  relative  to  the  earth). 

The  potential  temperature  of  air  is  divided  into 
two  parts:  a  basic  state  0B.(z)  and  perturbation  8',. 
The  basic  state  is  given  by 


where  DT0  is  the  mean  surface  temperature  differ¬ 
ence  across  the  MIZ  and  f(/)  is  the  nondimension- 
al  displacement  of  the  ice  edge.  The  coordinates 
and  atmospheric  variables  are  nondimensionalized 
by  setting 

(x.,z.,t.)  =  (xL.zb.tT), 

s'  *  0.V0BO  =  (DTo/Swts  (3) 

( U.,v.,w .)  =  U(u,  v,  wb/L), 

P.  =  QAgbDTJ8^P, 

where 

U  =  gbDT„/(2LW)m)  (4) 

is  the  scale  of  the  icebreeze  wind  and  T  is  the  time 
scale  for  the  change  of  surface  heating  due  to  the 
movement  of  the  ice  edge.  If  we  assume  that  the 
local  air  flow  satisfies  the  modified  Boussinesq  ap¬ 
proximation  (Kuo,  1973),  then  the  vorticity  equa¬ 
tion,  the  momentum  equation  (both  in  they' direc¬ 
tion),  and  the  heat  equation  for  air  disturbances 


y 


Fig.  3.  The  model  MIZ  and  the  coordinate  system. 


generated  by  differential  heating  near  the  M1Z  are 
written  by  (Chu,  1985) 


(12) 


Lim  ( !  1 ,  \dyp/dz\ ,  i  v’l ,  si)  <  oo. 


£Tf  V;  =  /„  dv/dz  -  ds/dx, 

EV2  v  =  -fa  d\J,/dz, 

I 

£T2s  =  R,d\J,/dx, 

and 

u  =  -  dyp/dz,  w  =  dyp/dx, 

! 

i  VJ  =  82  vii/(L2v)d1/dx1  +  d2/dz2, 

i 

!  rj  =  8VL2  dvdx2  +  d2/dz2  =  d2/dz2. 


(5) 

The  boundary  conditions  at  z  =  0  are: 

(6) 

yp  =  0, 

d<P/dz  =  Kd2yP/8z2, 

(13) 

(7) 

v  =  Kdv/dz, 

S  =  -[1  +  f(01  COS7TJC, 

where  K  is  a  measure  of  the  effective  depth  of  the 
constant  stress-sublayer  (Kuo,  1971).  Substituting 
Equation  1 1  into  Equation  10  we  obtain  the  fol¬ 
lowing  sixth-order  ordinary  differential  equation 
for  the  Fourier  coefficient  yp: 

(8) 

dty/df-Aydty/d?  +  4  (fi  +  y2) 


|  where /0  =  sinX,  X  is  the  latitude,  and 

1 

E  =  v/(2Q82)  =  Vi ,  /?j  =  SJ/W(4L!fiJ)  (9) 


cPl/dz^WRii  =  0 


where 


(14) 


are  the  Ekman  number  and  global  Richardson 
number,  respectively.  Eliminating  i>  and  s  from 
Equations  5  through  7,  we  find  that  the  stream- 
function  satisfies  the  following  partial  differential 
equation: 

( w  v;  +/o2)  d2i/dz2  +  r,  d^/dx2  =  o.  (io) 

We  solve  Equation  10  for  streamfunction  \P,  and 
obtain  the  solutions  of  v  and  s  from  Equations  6 
and  7  after  substituting  \P.  The  local  air  flow  is 
thermally  forced  by  differential  surface  heating  as 
indicated  in  Equation  2;  therefore,  the  stream- 
function  should  be  written  as 

i(x,  z,t)  =  i  (z,  t)  sin  nx.  (11) 

The  boundary  conditions  in  the  vertical  direction 
are  derived  as  follows.  The  dependent  variables 
should  remain  finite  as  z  —  oo,  i.e.. 


y  =  xJ62  vh/(2 L2v).  (15) 

The  general  solution  of  Equation  14  has  the  fol¬ 
lowing  form: 

6 

Uz,t)  =  [1 1-|-(01 V  5jeVz,  (16) 


where  tjie  eigenvalues  Xj  (j  =  1,2 . 6)  are  the 

roots  of  the  succeeding  sixth-order  algebraic  equa¬ 
tion: 

X‘-4yX4  +  4  (fa1  +  y2)\‘-4x2R,  =  0.  (17) 

Table  1  lists  all  the  eigenvalues  under  typical  strat¬ 
ification  at  three  distinct  latitudes  (65°,  70°,  and 
75  °).  Notice  that  fi  appears  in  Equation  17,  so  the 
eigenvalues  in  Table  1  represent  both  northern  and 


Table  1.  The  eigenvalues  of  the  boundary  layer  air  flow  model  at 
three  different  latitudes. 


_ Latitude _ 

65°  70°  75° 


(-0.32371,  0) 
(-1.02800,  0.87796) 
(-1.02800,  -0.87796) 
(  0.32371,  0) 

(  1.02800,  0.87796) 
(  1.02800,  -0.87796) 


(-0.31223,  0) 
(-1.04508,  0.89586) 
(-1.04508,  -0.89586) 
(  0.31223,  0) 

(  1.04508,  0.89586) 
(  1.04508,  -0.89586) 


(-0.30378,  0) 
(-1.05823,  0.90973) 
(-1.05823,  -0.90973) 
(  0.30378,  0) 

(  1.05823,  0.90973) 
(  1.05823,-0.90973) 


southern  polar  regions.  According  to  the  upper 
boundary  conditions  listed  in  Equation  12  we 
must  set  up  coefficients  that  correspond  to  those 
eigenvalues  with  positive  real  parts  to  zero.  Conse¬ 
quently  the  general  solution  Equation  16,  satisfy¬ 
ing  the  top  boundary  conditions,  may  be  written 
by: 


i(z,t)  =  |l  +  f(/)]^  a^f,  (18) 

7=1 

where  the  eigenvalues  X;  all  have  the  negative  real 
parts.  Substituting  Equation  18  into  11  we  obtain 
the  streamfunction 


3 


^  =  ° 


7=1 

3 


7=1 

3 


-/o^Xjd-KXj)  FfX J,  -y)5j 
7=1 

+  6(1  +  Ks/T',)  =  0 


(24) 


(25) 


(26) 


2 


v  (A,  C,  I)  =  y  fljll  +  v  (/)l  e'/  v/nir.v.  (19) 


i  i 


xR,  ^  F(Xj,  7)5j  +  fljfo  b/r  =  -1.  (27) 


7=1 


Integrating  the  momentum  elation  (6)  and  heat 
equation  (7)  with  respect  to  z  after  substituting 
Equation  19,  we  find  that  v  and  s  are  given  by 


v(jt,  z,  i ) 


7=1 

F{\,,y)  ex/][l  +  f(f)l  sinxjr, 


(20) 


sfx.  z,  t) 


7  =  1 


/^Xj.ijeV'HI  +  {‘(01  costt.v. 


(21) 


where 

F\Ky)  m  l/(XV2->).  (22) 

Substituting  solutions  19-21  into  the  vorticity 
equation  (5)  we  find  that 

c  =  vly/oS/w.  (23) 


Substituting  solutions  19-21  into  the  surface 
boundary  conditions  listed  in  Equation  13,  we  ob¬ 
tain  the  following  four  algebraic  equations  for  a) 
and  b. 


The  four  constants  d,,  3j,  a,,  and  b can  be  readily 
obtained  by  solving  the  four  linear  nonhomo- 
geneous  algebraic  equations. 

Figures  4-6  reveal  the  solutions  for  zero  ice  mo¬ 
tion  at  X  =  75°,  DT„  =  6°C,  and  N  =  10  J  s  '. 
Figure  4  shows  the  distribution  of  air  temperature 


Fig.  4.  The  distribution  of  air  temperature 
(°K)  in  the  cross-ice  section. 
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Fig.  5.  The  distribution  of  thermally  forced 
along-ice-edge  winds  (m  s  ')  in  the  cross-ice 
section. 


Fig.  6.  Thermally  forced  icebreeze  circulation  in  the 
cross- ice  section  with  u  in  m  s  '  and  w  in  cm  s  '. 


i 

i 


in  the  cross-ice  section.  The  horizontal  tempera¬ 
ture  gradient  is  very  deep  near  the  ice  edge  and  de¬ 
creases  with  altitude.  Above  z  =  3  (z.  =  800  m) 
there  is  almost  no  horizontal  temperature  gra¬ 
dient.  Figure  5  shows  the  distribution  of  the  wind 
along  ice-edge  v  in  the  cross-ice  section.  The  long¬ 
shore  wind  is  equatorward  at  low  levels  (below  z 
=  2)  and  poleward  at  high  levels.  This  strong 
equatorward  wind  along  the  ice  edge  near  the  MIZ 
(around  7  ms"')  is  consistent  with  the  observations 
along  the  eastern  Antarctic  Peninsula  (Schwerdt- 
feger,  1979).  Figure  6  shows  the  cross-ice  circula¬ 
tion.  The  maximum  value  of  the  ice  breeze,  which 
is  around  3  ms'1,  is  located  at  the  ice  edge. 

Ice  drift  model 

Suppose  that  the  MIZ  is  considered  to  be  pack 
ice  with  linear  viscous  rheology  and  constant  shear 
and  bulk  viscosities.  The  MIZ  has  a  mean  thick¬ 
ness  of  H0.  The  linearized  momentum  equation 
and  continuity  equation  are  written  by: 

itu/Ht ,  =  (i  +  >t)  il'u/Hx •’ 


+/Vj  +  (rax  +  T„x)/(Q,Ho)-ghhx/dx,  (28) 


dv/dt .  =  r/  dJVj/dxJ-/Mj 

+  (Tay  F  7wy)^(  P  JF) 

(29) 

dh/dt .  =  -(Ho/L)  du/dx. 

(30) 

where  /(  =  2ft  sin  X)  =  Coriolis  parameter 
Pi  =  ice  density 
Wj, Vj  =  Vi(  ice  velocity 
7ax»Tay  =  7a»  wind  Stress 
7wx>T*y  =  Tw.  water  stress 

hx  =  ice  thickness  perturbation 
generated  by  forcing  terms 
x  =  nondimensional  horizontal 
coordinate 

and  rj  and  I  are  defined  by 

S  ■  {/(Pi//o L1),  ij  =  r,/(e,HoL!),  (31) 

where  v  and  £  are  the  shear  and  bulk  viscosities  of 
ice.  The  air  and  water  stresses  are  given  by 

ra  =  Ci  V„  (32) 

r»  =  C;(VW-V  i),  (33) 
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where  Va  is  the  surface  air  velocity,  V*  the  surface 
water  velocity,  and  Ca  and  C*  the  dimensional  air 
and  water  drag  coefficients.  Ca  is  further  given  by 

Ca  =  CaQaU[(d^/dZ)!  +  V'o]  ;=0.x  =  0.5  .  (34> 

where  C o  and  v0  are  obtained  by  Equations  19  and 
20.  In  contrast  to  the  treatments  by  Hibler  (1984), 
Lepparanta  and  Hibler  (1984),  and  Lepparanta 
(1984),  the  surface  air  and  water  stresses  ra  and  rw 
are  directly  computed  by  the  bojndary  layer  flow 
model  presented  in  the  previous  section.  Substitut¬ 
ing  the  expressions  for  air  and  water  stresses  (32) 
and  (33)  into  Equations  28  and  29  and  assuming 
that  V*  =  0,  we  obtain 

du/dt,  =  (I  4  v)  d'-u/dx^Qu,  4  fv { 

Qua0[  1  4  f(/)]  sin  xx-g/L  dh/dx,  (35) 

dv/dl,  =  ijd'Vj/A^-QVj/uj 

+  Qva0(  1  +  . C(/)]  sin  «,  (36) 

where 

q  S  Ca/(ci«o),  C  =  Cy(Q,H0), 

3 

wao  s  -gbDTa/(2LSU>w) 

y  =  i 

(37) 
3 

va0  -  g6DTa/(2LQem)[b-fo  V  fljXjftXj.7)]. 

fa 

Equations  30,  35,  and  36  are  the  basic  equations 
for  ice  motion.  The  solutions  of  these  two  equa¬ 
tions  should  have  the  same  Fourier  component  as 
the  forcing  terms,  i.e., 

u,(x,t.)  =  u,(t.)  sin**, 

(38) 

Vj(jr,r.)  =  v,(t,)  sin  kx. 

By  definition,  f(f)  is  the  nondimensional  ice-edge 
displacement  in  the  x-direction.  It  is  related  to  ice 
velocity  by 


LdS/dt.  =  uy/ij.)  =  «,(/.).  (39) 

Eliminating  two  of  the  three  dependent  variables, 
w„  vif  and  h„  from  Equations  30,  35,  and  36  we  get 


the  following  third-order  equation  for  the  coeffi¬ 
cients: 

[D5  +  («,  +  <*;)£>•’  +  (x'gHo/ L 

4  a, nr 2  4 P)D  4  a}(itlgHa/ L!-Quaa/ L) 
-/Q va0/ /.](«;,  v j)  =  0,  (40) 


where 

D  =  d/di„ 

a ,  =  7TJ(?  4  7))  4  q, 

a,  =  *Vq 


(41) 


The  general  solutions  of  Equation  40  have  the  fol¬ 
lowing  form: 


J 

«,(/.)  =  ^  d]  exp(/i/.), 


j- 1 
3 


(42) 


i\(U  =  es  exp^j/.)- 
j  - 1 


The  eigenvalues  /*,,  n,,  n,  are  the  roots  of  the  fol¬ 
lowing  equation: 

H1  +  (oi  +  at2)fi2 

4  ( Tr2gH0/L2-Qu30/L  4  a, a:  4 f)n 
+  aiWgH„/L!-Qua0/L)-fQva0/L  =0.  (43) 

Stability  and  oscillation  criteria 

In  the  present  research  the  bulk  ice  viscosity  £  is 
taken  to  be  twice  the  shear  ice  viscosity  »/.  The 
standard  values  of  the  parameters  are  given  in 
Table  2.  We  solve  the  characteristic  equation  (43) 
with  different  values  of  H0,  y,  and  DT0.  H0  varies 
from  0.5  to  10.5  m,  rj  from  0  to  2  x  10"  kg  s'1,  and 
DTo  from  1  °C  to  21  °C.  The  instability  criterion  of 
the  MIZ  migration  in  the  ice/air-coupled  model 
can  be  written  by 


Re(/t) 


<  0  decaying 

=  0  neutral  £  =  (44) 

>  0  growing 
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Table  2.  The  standard  model  parameters. 


m 


i’|V| 

I#1 


L  =  200  km 
«  =  0.729  x  10  *  s  ' 
N  =  10  1  s' 
q  =  1.2x10’ 
t>,  =  910  kg  m’ 


i’  =  5  m1  s  '1 
g  =  9.81  m1  s'! 
flao  =  270  °K 
ea  =  1.29  kg  m’ 

Uj  <  1ms'1 


i'h  =  10’  m;  s-‘ 

K  =  0.9 

C'w  =  0.55  kg/m!  s"‘ 
ew  =  10’  kg  m’ 

X  =  (65  °,70o,75°) 


f 


»’,Vl 


where  y  is  the  root  of  the  cubic  equation  (43).  The 
oscillation  criterion  of  the  MIZ  migration  is  given 
by 


=  0  nonoscillatory 
*  0  oscillatory 


>  (45) 


Dependence  of  eigenvalues  on  parameters 

We  compute  all  the  roots  of  Equation  43.  Here 
u,  is  real  in  the  entire  space,  and  y2  and  y>  are  real 


for  some  values  of  (H„,  ij,  DTo)  and  are  complex 
conjugates  for  others. 

Figure  7  shows  the  surface  of  y,  =0  in  three- 
dimensional  space  {Ho,  v,  DT0).  This  surface  di¬ 
vides  the  space  into  two  parts  corresponding  to 
growing  and  decaying  modes.  The  growing  mode 
appears  when  DT0  exceeds  some  critical  value 
around  2  °C,  and  is  located  in  the  small  mean  ice- 
thickness  Ho  region.  Ice  motion  corresponding  to 
the  eigenvalue  yt  is  nonoscillatory. 

Figure  8  indicates  the  surface  of  Re(/i2)  =  0  [or 


‘/vt’L 


<  o 


ft*  °' 


Non-osclllatory 
Crowing  Modes 


DT0(°C) 


Fig.  7.  Separation  of  nonoscillatory  decaying  and  growing  modes  referring  to  real 
eigenvalue  y ,. 
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/■/,?.  <8.  Separation  of  decaying  and  growing  modes  referring  to  eigenvalues  io  and  70. 


Re(#«.)  =  0)  in  the  parameter  space  (H0,  17,  DT„). 
This  surface  separates  the  whole  space  into  grow¬ 
ing  and  decaying  parts.  The  growing  mode  is  pres¬ 
ent  when  DTo  exceeds  some  critical  value  that  is  a 
function  of  Ho  and  >7  (in  case  of  H«  -  4.5  m  and  >7 
=  1 0“  kg  s  ',  DTo  =  6°C),  and  is  concentrated  in 
the  large  mean  ice-thickness  Ho  region. 

Figure  9  reveals  the  segregation  of  nonoscilla- 
tory  and  oscillatory  modes  relating  to  70  and  70. 
Comparing  Figure  9  with  Figure  8,  we  find  that 
the  decaying  mode  of  70  and  70  is  nonoscillatory, 
whereas  the  growing  mode  of  71.  and  70  is  oscilla¬ 
tory. 

Whether  ice  motion  grows  or  decays  substan¬ 
tially  depends  on  the  parameter  DTo.  When  DTo  is 
small,  the  motive  force  is  so  small  that  it  cannot 
overcome  the  dissipative  effect  and  does  not  make 
ice  motion  unstable.  If  DTo  becomes  large  enough 
to  overcome  the  dissipative  forces  originated  by 
ice  stress  and  friction,  ice  motion  becomes  larger. 


Whether  ice  motion  is  oscillatory  or  nonoscilla¬ 
tory  largely  depends  on  parameter  H„.  In  cases  of 
small  Ho,  the  restoring  force  is  so  small  that  it  can¬ 
not  balance  the  motive  force  and  cannot  produce 
oscillations.  However,  when  parameter  H0  exceeds 
some  critical  value  (which  varies  with  DTo  and  >7) 
and  becomes  large  enough  to  balance  the  motive 
force  brought  on  by  the  ice  breeze  generated  by 
surface  temperature  difference,  the  ice  motion  be¬ 
comes  oscillatory. 

For  the  nonoscillatory  growing  mode,  the  time 
during  which  ice  doubles  its  speed  is  computed  by 

T  =  In  2/71,  (70  >  0).  (46) 

If  Ho,  »).  and  £  are  assigned  as  commonly  used  val- 


10“  kgs 


£  =  2x10“  kgs 


DT0(°C) 


Ira(Vo)  >  0 


I»(p,)  <  0 


Oscillatory 


INI  On  I 


Nonos  dilatory 


Fig.  9-  Separation  of  nonosci/latorv  and  oscillatory  states  referring  to  eigenvalues  /t: 
and  /ii. 

Z5| - - 

The  doubling  time  treated  as  a  function  of  DT0  is 
shown  in  Figure  10,  which  displays  a  decrease  of 
T i  with  an  increase  of  DT„.  T,  changes  from  7.3  hr  20  -  \ 

to  0.73  hr  as  DT„  varies  from  6°C  to  16  °C.  It  takes  \ 

only  12.4  hr  for  the  ice  velocity  to  increase  from  6  1 

cm  s'1,  which  is  a  steady-state  estimation,  to  20  cm  „  is  -  \ 

s  '  when  DT„  =  6°C.  |  l 

The  oscillatory  growing  mode  is  ice  motion  due  l 

to  n;  (or  g}),  whose  real  part  is  positive.  10  -  \ 

Re(//;)  =  Re(n>)  >0,  V 

the  growth  rate  Re(/i2  3)  increases  with  an  increase  \. 

of  DTa  (increase  of  forcing  term),  with  an  enlarge-  _ 

ment  of  Ha  (increase  of  restoring  force),  and  with  0 - 1 - 1 - 1 - 1 - * - 1 - 

a  reduction  of  r?  (decrease  of  dissipation).  With 

DT„  =  6°C,  rj  =  10’  kg  s  ',  and  H„  =  4.5  m,  the  „ 

growth  rate  is  0.22  x  10'4  s  '  and  the  corresponding  0 

doubling  time  is  8.75  hr.  The  period  of  growing  Fig.  10.  Dependency  of  doubling  time  on  DT 

oscillatory  mode,  defined  by  for  the  nonoscillatory  growing  mode. 


T  (hr) 


Fig.  II.  Dependency  of  doubling  time  and  period 
on  DTo  for  oscillatory  growing  mode. 


Tp  =  271-/ 1  Im(/i2.j)  I  (47) 

is  9.7  hr. 

Figure  1 1  shows  the  growth  rate  and  period  as 
functions  of  DT '*  for  Ho  =  4.5  m  and  r\  =  10’  kg 
s'1.  The  growth  rate  decreases  monotonically  with 
an  increase  of  DT0.  However,  the  period  decreases 
slightly  when  DT0  varies  from  6°C  to  10  °C,  and 
then  increases  with  DT0  rapidly.  When  DT0  = 
19  °C,  the  period  is  nearly  one  day. 

Conclusion 

This  ice/air-coupled  model  is  intended  only  to 
depict  the  mesoscale  processes  of  ice/aii  interac¬ 
tion  in  the  M1Z.  The  synoptic  scale  pressure  gra¬ 
dient  associated  with  the  semi-permanent  Iceland¬ 
ic  low  may  in  addition  produce  equatorward  winds 
in  the  Greenland  Sea,  and  large-scale  ocean  cur¬ 
rents  near  the  East  Greenland  Sea  or  the  Antarctic 
Peninsula  may  drive  ice  drift.  These  processes  are, 
however,  beyond  the  scope  of  the  current  re¬ 
search.  Nevertheless,  when  the  ice  to  open  water 
temperature  gradient  is  strong,  the  mesoscale 
feedback  mechanism  discussed  here  may  become 
as  strong,  or  stronger,  than  the  synoptic  scale  and 
oceanic  forcings. 

The  equatorward  and  waterward  surface  winds 
near  the  MIZ  in  the  Greenland  Sea  and  the  east 
Antarctic  Peninsula  are  thermally  generated  by 
differential  surface  heating  over  ice  and  water. 
The  surface  wind  has  its  maximum  speed  along  the 


ice  edge.  The  component  of  surface  wind  parallel 
to  the  ice  edge  can  reach  7  ms'1  when  the  tempera¬ 
ture  gradient  in  the  cross-ice  direction  is  chosen  as 
0.03 °K  km  '.  This  agrees  with  the  observational 
results  at  the  east  Antarctic  Peninsula  reported  by 
Schwerdtfeger  (1979). 

The  migration  of  the  MIZ  is  forced  by  local 
winds  generated  by  differential  surface  heating 
over  ice  and  water  near  the  ice  edge.  The  ice  veloc¬ 
ity  V,  has  sinusoidal  forms  (38),  showing  an  ice 
edge  jet.  The  icv  velocity  has  its  maximum  at  the 
ice  edge  and  decreases  iceward.  This  implies  that 
the  ice  edge  jet  is  produced  indirectly  by  the  ther¬ 
mal  effect  of  the  ice/water  contrasts. 

The  ice  motion  has  two  bifurcations.  It  first  bi¬ 
furcates  into  decaying  and  growing  modes,  de¬ 
pending  on  th  mean  surface  temperature  differ¬ 
ence  over  ice  and  water,  DT0.  When  DT0  is  small, 
the  decaying  mode  predominates.  However,  when 
DTo  exceeds  some  critical  value  (i.e.,  6°C  as  H„  = 
4.5  m  and  t)  =  10’  kg  s'1,  and  2°C  as  Ho  =  1.5  m 
and  j]  =  10’  kg  s'1),  the  growing  mode  appears. 
The  growing  mode  further  bifurcates  into  nonos- 
cillatory  and  oscillatory  states,  depending  on  ice 
properties.  In  cases  of  small  ice  thickness,  which 
causes  an  inadequate  restoring  force  and  generally 
occurs  in  summer,  the  ice  motion  is  nonoscilla- 
tory.  In  cases  of  large  ice  thickness,  which  gener¬ 
ates  an  adequate  restoring  force  and  is  usually 
present  in  winter,  the  ice  motion  is  oscillatory. 

This  research  provides  a  possible  way  to  predict 
the  MIZ  migration.  The  procedure  is  as  follows: 

a)  Observe  surface  temperature  difference  DT„ 
and  determine  the  ice  properties,  i.e.,  H0  and  t). 

b)  Use  Figures  8  and  9  to  determine  the  pattern 
of  motion,  i.e.,  growing  or  decaying,  oscillatory 
or  nonoscillatory,  etc. 

c)  Employ  Figures  10  and  11  to  determine  the 
growth  rate  and  period  (in  case  of  the  oscillatory 
state). 

d)  Observe  initial  ice  velocity. 

e)  Make  a  prediction  of  MIZ  migration  accord¬ 
ing  to  the  initial  ice  velocity,  growth  rate,  and  peri¬ 
od  (for  the  oscillatory  state). 

The  author  is  grateful  to  Prof.  H.L.  Kuo  and 
Prof.  D.R.  MacAyeal  of  the  University  of  Chi¬ 
cago  for  invaluable  discussion  and  comments.  The 
reviewers’  comments  are  also  highly  appreciated. 

This  research  was  supported  by  grant  ATM 
83-14206  from  the  National  Science  Foundation. 


REFERENCES 

Brown,  R.A.,  Planetary  boundary  layer  modeling 
for  AIDJEX.  In  Sea  Ice  Processes  and  Models. 


63 


University  of  Washington  Press,  387-401,  1980. 
Chu,  P.C.,  A  boundary  layer  theory  of  coastal 
desert  formation  and  seabreeze  circulation.  Ph.D. 
thesis.  The  University  of  Chicago,  194  pp.,  1985. 
Glen,  J.W.,  Thoughts  on  a  viscous  model  for  sea 
ice.  A/DJEX  Bull.,  2,  18-27.,  1970. 

Hibler,  W.D.,  III,  Ice  dynamics.  CRREL  Mono¬ 
graph  84-3,  52  pp.,  1984. 

Kuo,  H.L.,  Axisymmetric  flows  in  the  boundary 
layer  of  a  maintained  vortex.  J.  Atmos.  Sci.,  28, 
20-41,  1971. 

Leavitt,  E.,  Surface-based  air  stress  measurements 
made  during  AIDJEX.  In  Sea  Ice  Processes  and 
Models.  University  of  Washington  Press, 
419-429,  1980. 

l.epparanta,  M.  and  YV.I).  Hibler  111,  On  the  role 
of  ice  interaction  in  marginal  ice  zone  dynamics. 
MIZEX  Bulletin  III,  USA  Cold  Regions  Research 
and  Engineering  Laboratory,  Special  Report  84-7, 
pp.  23-29,  1984. 

Macklin,  S.A.,  Wind  drag  coefficient  over  first- 
year  sea  ice  in  the  Bering  Sea.  J.  Geophvs.  Res., 
88,  2845-2952,  1983. 

McPhee,  M.G.,  The  effect  of  the  oceanic  boun¬ 
dary  layer  on  the  mean  drift  of  pack  ice:  Applica¬ 


tion  of  a  simple  model.  J.  Phvs.  Oceanogr. ,  9, 
388-400,  1979. 

Overland,  J.E.,  H.O.  Mofjeld  and  C.H.  Pease, 

Wind-driven  ice  drift  in  a  shadow  sea.  MIZEX 
Bulletin  III,  USA  Cold  Regions  Research  and  En¬ 
gineering  Laboratory,  Special  Report  84-7,  pp. 
49-53,  1984. 

Pease,  C.H.,  S.A.  Salo,  and  J.E.  Overland,  Drag 
measurements  for  first-year  sea  ice  over  a  shallow 
sea.  J.  Geophvs.  Res.,  88,  2853-2862,  1983. 

Prik,  Z.M.,  Mean  position  of  surface  pressure  and 
temperature  distribution  in  the  Arctic.  Tr.  Ark- 
ticheskogo  Nauchn.-lssled.  Inst.,  217,  5-24  (in 
Russian),  1959. 

Schwerdlfeger,  W„  Meteorological  aspects  of  the 
drift  of  ice  front  the  Weddell  Sea  toward  the  mid¬ 
latitude  westerlies.  J.  Geophvs.  Res.,  84,  6321- 
6328,  1979. 

Smith,  S.D.,  E.G.  Banke,  and  O.M.  Johannessen, 

Wind  stress  and  turbulence  over  ice  in  the  Gulf  of 
St.  Lawrence.  J.  Geophvs.  Res.,  75,  2803-2812, 
1970. 

Vinje,  T.E.,  The  drift  pattern  of  sea  ice  in  the  Arc¬ 
tic  with  particular  reference  to  the  Atlantic  ap¬ 
proach.  In  The  Arctic  Ocean,  John  Wiley  and 
Sons,  New  York,  83-96,  1982. 


The  Planetary  Boundary  Layer  in  the  Marginal  Ice  Zone 
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The  object  of  the  Marginal  Ice  Zone  Experi¬ 
ments  (MIZEX)  meteorology  analysis  program  is 
to  obtain  information  for  all  scales  of  atmospheric 
motion  that  are  directly  affected  by  the  marginal 
ice  /one  (M1Z)  or  vice  versa.  The  relatively  abrupt 
transition  from  pack  ice  to  ocean  surface  provides 
an  excellent  location  for  investigating  the  air  mod¬ 
ification  in  the  Planetary  Boundary  Layer  (PBL) 
How  and  the  effects  on  the  surface.  The  most  sig¬ 
nificant  parameter  of  interaction  between  air  and 
pack  ice  or  sea  surface  is  the  aerodynamic  drag 
parameter,  Q}  =  t/qU'.  This  is  best  obtained 
from  measurements  of  stress,  r  and  U,  at  the  sur¬ 
face.  However,  the  measurement  of  stress  is  diffi¬ 
cult  and  requires  fairly  complicated  measure¬ 
ments.  Consequently,  stress  data  over  the  pack  ice 
or  the  ocean  has  rarely  been  obtained  and  drag  co¬ 
efficients  established  at  smooth  land  sites  were 
used.  Within  the  past  decade,  experiments  have  re¬ 
vealed  new  trends  in  the  variation  of  the  drag  coef¬ 
ficient  over  these  surfaces. 

One  of  the  goals  of  air/sea  interaction  model¬ 
lers  has  been  to  get  an  upper  boundary  condition 
for  the  oceanic  flow,  the  stress.  Traditionally,  this 
upper  boundary  condition,  the  surface  stress,  was 
used  as  a  free  parameter  in  dynamic  models.  The 
oceanic  problem  was  simply  solved  or  modelled 
with  respect  to  the  input  forcing  represented  in  the 
drag  coefficient,  and  this  constant  could  be  looked 
up  somewhere.  In  the  case  of  a  large  experiment, 
all  that  was  needed  was  to  find  a  wind  measure¬ 
ment  in  the  region,  square  it,  multiply  by  constant 
density,  and  obtain  an  estimate  of  the  stress  for 
the  entire  area.  This  had  the  advantage  that  if 
things  didn’t  work  out  right  in  the  model,  one 
could  always  change  the  Cp  as  a  free  parameter. 
This  was  probably  good  enough  in  some  cases, 
and  it  may  have  been  justified  by  the  difficulty 
that  atmospheric  scientists  seemed  to  be  having  in 
determining  the  value  of  Cp  over  (he  ocean.  How¬ 
ever,  there  are  many  current  oceanic  models  that 


deal  with  small  scales  and  are  accurate  enough  to 
respond  to  variations  in  Cp  that  naturally  occur. 

In  field  experiments,  the  stress  is  usually  ob¬ 
tained  at  a  point  using  sophisticated  turbulence 
measurements,  theory,  and  analysis.  The  eddy 
correlation  method  requires  high-frequency  turbu¬ 
lence  measurements  to  determine  the  time  average 
of  the  velocity  fluctuations  correlation,  u'w'  = 
t/q.  This  uses  sophisticated  instruments  that  are 
difficult  to  maintain  in  the  harsh  environment  of 
the  MIZ.  Nevertheless,  enough  measurements 
have  been  made  over  pack  ice  to  give  an  approxi¬ 
mation  of  the  stress  based  on  a  single  mean  wind 
measurements  using  the  drag  coefficient  method. 
The  problem  is  then  to  extrapolate  these  data  to 
the  mesoscale  and  synoptic  scale.  This  is  done  with 
a  mixture  of  theory  and  data.  The  theory  defines 
similarity  parameters  that  can  relate  point  values 
of  wind,  stress,  and  heat  flux  to  large-scale  mean 
flow  parameters,  such  as  horizontal  pressure  and 
temperature  fields.  Unfortunately,  most  of  the 
similarity  parameters  are  not  universal.  They  vary 
with  PBL  stratification  and,  to  a  lesser  degree, 
with  surface  roughness  and  windspeed.  All  of 
these  variations  become  particularly  significant  in 
the  MIZ,  since  these  factors  change  dramatically 
in  this  region. 

The  model 

In  applications,  initially  the  drag  coefficient  is 
often  assumed  to  be  constant.  The  surface  stress 
then  depends  only  on  the  surface  wind  speed,  den¬ 
sity  being  essentially  constant  within  the  boundary 
layer.  Moreover,  since  surface  wind  measurements 
are  generally  not  widely  available,  it  is  often  neces¬ 
sary  to  obtain  the  surface  wind  from  the  geostroph- 
ic  windfield.  The  latter  can  be  approximated  from 
a  surface  pressure  field  by  assuming  vertical  co¬ 
herence  of  the  pressure  gradient.  This  condition  is 
a  good  approximation  as  long  as  significant  baro- 
clinicity  (the  horizontal  temperature  gradient) 


doesn’t  exist.  If  it  does,  then  a  value  for  the  hori¬ 
zontal  temperature  field  is  needed  to  calculate  this 
influence  (the  thermal  wind).  FBL  similarity  theory 
provides  a  direct  link  between  the  synoptic  flow, 
G,  and  the  surface  stress.  The  surface  friction  vel¬ 
ocity,  u*  -  \Jt/q,  can  be  used  to  define  the  geo- 
strophic  drag  parameter,  u*/G.  The  angle  of  turn¬ 
ing  between  the  geostrophic  flow  and  the  surface 
wind  or  stress,  a,  can  also  be  calculated.  All  that  is 
required  is  a  set  of  good  measurements  of  u*,  G, 
and/or  fj,0  to  evaluate  the  similarity  parameters 
A  '  and  B  in 


(sina)  =  -B 

(1) 

kG  , 

~u»-  (cos a)  =  -A 

where  B  is  a  constant  and  A'  *  A  +  In  (ku*/fz0) 
is  a  function  of  surface  roughness,  as  represented 
in  the  roughness  parameter  z0.  k  is  von  Karman’s 
constant  -  0.4,  and /is  the  Coriolis  parameter. 


The  data 

To  illustrate  the  history  of  oceanic  drag  coeffi¬ 
cients,  Figure  1  is  a  montage  of  a  collection  of 
anonymous  published  results  of  Co  vs.  surface 
windspeed  and  stratification  parameter,  z/L, 
where  L  is  the  Obukhov  stratification  parameter, 
which  must  be  designated  along  with  CD.  Scatter  is 
a  problem  in  the  data.  Cp  is  not  constant;  it  varies 
with  windspeed,  stratification,  and  surface  rough¬ 
ness.  These  effects  have  been  slowly  identified  on¬ 
ly  from  many  stress  measurements.  There  is  varia¬ 
tion  in  Cq  with  magnitude  of  U  due  to  surface 
roughness  (wave  spectrum)  changes.  Stratification 
measurements  were  generally  not  available,  since 
they  require  measurements  of  vertical  temperature 
gradients,  or  at  least  air/surface  temperature  dif¬ 
ferences.  Early  graphs  of  Cp  vs.  roughness  conse¬ 
quently  were  even  more  scattered.  Most  plots  of 
Cjj  vs.  any  single  variable— stratification,  rough¬ 
ness,  or  windspeed — will  generally  show  large 
scatter  due  to  the  difficulty  in  eliminating  the 
other  effects.  Since  this  ancillary  data  is  not  gener¬ 
ally  available,  most  MIZ  and  polar  measurements 
have  intrinsic  large  error. 


0  10  20  30 

Ur  (m/s) 


fig.  /.  The  results  of  experimental  measurements  of  the  10-m  (Iran  coefficient  over  the  ocean  and 
the  roughness  parameter  z„  v.v.  mean  wind  speed  or  stratification  parameter. 
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Fig.  2.  A  representative  summary  of  drag  param¬ 
eters  measured  over  pack  ice. 

Figure  2  shows  representative  values  of  Cq 
taken  over  the  ice  during  the  past  few  decades.  A 
more  comprehensive  compilation  by  Overland 
(1985)  shows  the  same  trend.  The  drag  seems  to  be 
increasing  with  time.  This  increase  is  actually  due 
to  two  factors:  better  control  in  measurements, 
and  data  taken  in  increasingly  rough  ice  areas.  Up 
to  1975,  measurements  were  usually  taken  on  very 
large,  very  smooth  ice  floes.  During  AIDJEX 
(Arctic  Ice  Dynamics  Joint  Experiment),  several 
different  measurements  of  stress  were  taken.  The 
traditional  methods:  the  profile  and  eddy  correla¬ 
tion  with  sonic  anemometer  data,  all  were  taken 
on  a  very  large,  smooth  ice  floe,  and  yielded  val¬ 
ues  of  1-2  x  I0"3.  There  were  indications  of  higher 
values.  Banke  and  Smith  (1973)  took  measure¬ 
ments  in  rougher  areas  and  obtained  values  up  to 
2.4  x  10  \  Aircraft  data  taken  during  AIDJEX  in¬ 
cluded  large  horizontal  averages  and  yielded  ex¬ 
trapolated  values  of  CD  =  2.5  x  10'3  (Katz,  1980). 
Another  method,  which  incorporated  the  entire 
boundary  layer  wind  profile,  was  the  momentum 
integral  method  from  Brown  (1974).  This  was 
used  with  acoustic  sounder  data  by  Carsey  (1980) 
to  obtain  Cps  up  to  3.1  xIO"3.  In  recent  years. 
Pease  (personal  communication)  (using  a  momen¬ 
tum  integral),  Macklin  (1983)  (a  profile  tower), 
Walter  et  al.  (1984)  (aircraft),  and  Andreas  et  al. 
(1984)  have  worked  in  the  MIZ,  obtaining  large 
Cq  for  this  rough  region. 

During  the  1970’s,  several  large-scale  experi¬ 
ments  were  conducted  where  Cp  was  measured  ac¬ 
curately  in  high  winds  over  the  ocean.  Figure  3  is  a 
representative  data  set  from  one  of  these  experi¬ 
ments,  the  Joint  Air-Sea  Interaction  Experiment 
(JASIN),  calculated  by  Large  and  Pond  showing  a 
two-regime  behavior  of  Co  with  windspeed.  The 
drag  coefficient  increases  with  windspeed  at  high 
winds — the  most  important  case.  The  next  illustra¬ 
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Fig.  3.  Drag  parameters  calculated  in  the  JASIN 
experiment  on  the  weathership  Meteor.  The  solid 
line  is  an  average  found  at  the  Bedford  Institute  tower. 


tion,  Figure  4,  shows  a  similar  two-regime  behav¬ 
ior  of  heat  flux  with  respect  to  stratification. 
There  is  a  dramatic  change  across  neutral  stratifi¬ 
cation.  Since  momentum  is  transported  by  the 
same  turbulent  eddies,  one  might  expect  a  similar 
behavior  in  stress  variation  with  stratification. 
Heat  and  momentum  flux  (stress)  depend  on  the 
turbulence  distribution  in  the  PBL,  which  depends 
on  surface  roughness  (mechanical  mixing),  strati¬ 
fication  (via  convective  instabilities,  buoyancy  ef¬ 
fects),  and  windspeed  (via  dynamic  instabilities, 
secondary  flows). 

Using  the  observations  to  determine  the  similar¬ 
ity  parameters  resulted  in  large  scatter  in  the 
similarity  functions.  The  /4’s  and  B’s  must  be  de¬ 
termined  by  measuring  all  other  parameters  in 
Equation  1:  u*.  G,  k,  and  z0.  Since  they  were 
some  (unknown)  function  of  stratification,  a  pa¬ 
rameter  must  be  determined  from  measurements 


Fig.  4.  Same  as  Figure  3  for  the  heat  flux  param¬ 
eter  vs.  the  stratification  parameter  z/L. 
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Fig.  5.  Similarity  functions  vs.  stratification  par¬ 
ameter  n  =  lu*/fL.  Taken  from  Brown  ( 1982 ). 


Fig.  7.  Air  photo  of  ice  ridges  near  the  AlDJFX 
camp,  1976.  This  ice  ridge  is  about  2  m  high.  10  km 
long. 


Fig.  6.  Ice  feathers  on  pack  ice.  The  width  at  the 
nearest  point  is  about  /  meter. 


of  heat  flux.  Figure  5  shows  typical  data-produced 
A  and  B  vs.  stratification.  It  would  appear  that  we 
are  not  much  better  off  with  this  method.  How¬ 
ever,  the  scatter  in  this  formulation  may  often 
result  from  improper  representation  of  the  turbu¬ 
lence.  When  secondary  flows  in  the  Ekman  layer 
are  included  in  the  turbulence  analysis,  the  (modi¬ 
fied)  similarity  laws  offer  good  parameterization 
of  the  stress. 

Figures  6  through  8  show  some  variation  in  the 
surface  roughness  of  the  ice.  Figure  6  is  a  form  of 
ice  feathers  (frosting);  it  is  probably  not  common 
enough  to  worry  about,  but  it  is  illustrative  of  how 
varied  the  ice  surface  can  be.  Figure  7  shows  the 
pack  ice  near  AIDJEX.  This  is  a  relatively  smooth 
floe,  but  it  contains  ridges  1-2  m  high  separated 
by  a  few  hundred  meters  to  kilometers.  Figure  8 
shows  the  profile  tower  at  AIDJEX,  26  m  high, 
used  to  measure  stress.  Even  this  tall  (for  the  Arc¬ 
tic)  tower  reaches  only  about  5%  of  the  PBL,  and 
reflects  the  surface  roughness  characteristics  only 
in  the  vicinity  of  the  tower— a  relatively  smooth 
floe.  Figure  9  shows  stress  data  obtained  from  this 
tower,  plotted  versus  a  stratification  parameter, 
n  =  ku*/fL.  L  is  the  Obukhov  parameter  =  ku*'/ 
gH0\  H0  is  the  heat  flux.  The  lines  are  from  the 


Fig.  8.  AIDJEX  meteorological  tower.  It  was  26  in 
tall  anti  was  frequently  climbed.  The  smooth  ice  of  the 
large  (10-km)  drift  is  broken  only  by  the  black  drums  at 
the  fuel  depot. 
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Fig.  9.  Geostrophic  drag  parameter  values  ob¬ 
tained  from  the  AIDJFX  tower  by  E.  Leavitt. 
The  lines  are  from  the  theory  of  Brown  (1981/  for 
various  geostrophic  windspeeds  Ci  and  roughness 
parameters  /„. 


theory  of  Brown  (1981).  The  scatter  is  typical. 
There  is  a  clear  trend  toward  smaller  u*/G  for 
stable  stratification.  Since  stress  is  proportional  to 
u*\  the  change  in  stress  (dynes/cm2)  is  significant. 

The  general  picture  over  the  ocean  is  shown  in 
Figure  10.  Cp  vs.  surface  wind  magnitude  at  dif¬ 
ferent  air/sea  temperature  differences  are  shown 
for  many  experiments.  The  heavy  line  is  the  com¬ 
posite  parametric  curve  used  in  the  model  for  neu¬ 
tral  stratification.  The  heavy  dashed  lines  are  the 
predicted  curves  with  stratification  effects.  Again, 
the  data  did  not  have  good  stratification  informa¬ 
tion,  and  this  could  account  for  some  of  the  dif¬ 
ferences,  e.g.  the  Davidson  data  were  in  somewhat 
stable  stratification  (personal  communication). 

The  PBL  model 

Figure  1 1  illustrates  one  of  the  characteristics  of 
the  MIZ  PBL,  the  omnipresence  of  cloud  streets 
when  wind  >  7.0  m  s  ’.  In  the  past  decade  labora¬ 
tory  experimenters  have  used  new  flow  visualiza¬ 
tion  techniques  to  discover  a  fundamental  charac¬ 
teristic  of  turbulent  flow — the  fact  that  often 
large,  organized  eddies  or  secondary  flows  can  ex¬ 
ist  within  the  fully  developed  turbulent  flow  (e.g., 
Meyer,  1981).  This  happens  most  often  when  the 
system  is  in  a  rotating  frame  of  reference.  The 
clouds  provide  excellent  flow  visualization  in  the 
atmosphere,  and  the  rotating  frame  of  reference  is 
built  in.  Furthermore,  there  is  an  analytic  solution 
for  the  mean  flow:  the  Ekman  spiral.  There  is  also 
an  analytic  solution  for  the  mean  flow  containing 
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Fig.  10.  The  drag  parameter  vs.  windspeed  sum¬ 
marized  for  many  experiments.  The  model  function 
is  shown  Jor  neutral,  stable,  and  unstable  stratification. 
Taken  from  Brown  and  I  iu  (1982/. 


Fig.  11.  Satellite  photograph  of  northerly  flow  off 
the  pack  ice  over  the  ocean  in  the  Bering  Sea,  Feb¬ 
ruary  1981.  The  cloud  rows  occur  in  the  convergence 
and  upward  motion  of  the  underlying  secondary  flow 
roll  circulation  in  the  PBL.  The  wavelength  is  about  5 
km.  The  stratification  is  near-neutral  over  the  pack  ice 
and  convectively  unstable  over  the  ocean.  The  flow  be¬ 
comes  cellular  at  low  Reynolds  numbers  over  the  ocean. 


the  secondary  flow  large  eddies:  the  modified  Ek- 
man  spiral.  Figure  12  is  a  sketch  of  the  predicted 
flow  pattern  in  the  PBL.  Condensation  occurs  in 
the  updraft  regions,  where  the  convective  activity 
and  the  moisture  flux  are  enhanced.  The  lateral 
wavelength  of  these  vortical  eddies  is  about  2-6 
km,  depending  on  stratification  in  the  layer.  This 


means  that  there  is  a  significant  variation  in  the 
vertical  wind  profile  if  taken  at  any  point  in  the 
secondary  vortices,  slich  as  A,  B,  or  C  in  the  fig¬ 
ure.  The  expected  vertical  wind  profiles  are  shown 
in  Figure  13.  These  are  the  “mean”  winds,  which 
would  be  determined  from  a  radiosonde  launched 
into  the  PBL’s  vortical  eddies,  or  rolls.  Each  pro¬ 
file  is  different  from  the  mean  profile,  which 
would  result  from  averaging  across  many  rolls. 
One  expects  similar  variation  in  the  phase  of  heat, 
momentum,  and  any  passive  agent  such  as  COj. 
The  onset  of  secondary  flow  and  enhanced  flux  re¬ 
gions  as  the  wind  gradients  increase  may  account 
for  the  jump  in  momentum  and  heat  flux  indicat¬ 
ed  by  Figures  3  and  4.  While  it  is  impractical  to 
launch  enough  radiosondes  to  obtain  an  average, 
an  aircraft  can  traverse  the  rolls  and  readily  obtain 
average  flows. 

The  influence  of  the  boundary  reaches  quite 
deep  into  the  atmospheric  flow,  from  100s  of 
meters  to  several  kilometers,  depending  on  wind- 
speed,  surface  roughness,  and  stratification.  The 
scatter  in  point  measurements  with  radiosonde 
balloons  prevented  any  meaningful  understanding 
of  the  upper  PBL  flow.  Consequently,  most  data 
are  from  the  surface  layer.  The  logarithmic  veloci¬ 
ty  profile  and  its  variation  with  stratification  for 
this  layer  over  land  have  been  thoroughly  param¬ 
eterized. 

Our  boundary  layer  model  incorporates  both 
surface  and  modified  Ekman  boundary  layer  solu¬ 
tions,  patching  them  together  at  some  height 
around  15%  of  the  height  of  the  boundary  layer. 
Since  the  rolls,  and  consequently  the  mean  flow 
dynamics,  vary  with  convective  energy  input,  both 
layers  vary  with  stratification.  Figure  14  is  a 
sketch  showing  the  predicted  flows  in  each  layer. 
The  continuous  velocities  from  the  two-layer  mod¬ 
el  include  the  contributions  from  surface  layer 
stratification,  PBL  stratification,  surface  rough¬ 
ness,  and  thermal  wind,  as  given  by  Brown  and 
Liu  (1982).  Note  that  the  secondary  flows  are  a  re¬ 
sult  of  dynamic  instabilities  in  the  basic  Ekman 
flow,  and  are  always  present  above  a  threshold 
wind  velocity  of  about  7  m  s"1  neutral — less  when 
convection  is  important,  and  not  at  all  when  sig¬ 
nificant  stable  stratification  exists  (local  Ri  > 
0.25). 

When  the  equations  for  patching  the  modified 
Ekman  solution  to  the  log  solution  are  written,  the 
similarity  characteristics  of  the  PBL  emerge.  The 
“classic”  similarity  relations  for  PBL  parameter¬ 
ization  of  the  neutral,  stable,  barotrophic  layer  are 
given  in  (1).  When  stratification  effects  are  includ¬ 
ed,  the  similarity  equations  become  (Brown,  1981) 


Typicol  Secondary  Flow 

in  the  Planetary  Boundary  Layer 

(Modified  Ehman  Layer) 


Fig.  12.  Sketch  of  the  coherent  helical  flow  vortices  contained  in  the  boundary  layer. 


Fig.  13.  Typical  wind  hodographs  to  be  expected 
from  radiosondes  taken  at  stations  A,  B,  or  C  in 
Figure  12.  The  freestream  flow  is  Ug,'  the  heights  are 
with  respect  to  the  Ekman  depth  f>. 


(sina  +  /3)  ®  -B 
kG  , 

(cosa  +  y)  =  -A 

where  B  =  (l-X)^and/!  '  =  ln|-2£u*/[/i0(l-^)]l 
+  i p-B.  0  and  -y  are  corrections  that  are  non-zero 
whenever  secondary  flows  or  baroclinicity  are  im¬ 
portant.  4/  is  the  stratification  correction  to  the 
surface  layer  (Paulson,  1970).  In  the  general  case 
with  stratification  the  predicted  drag  law  can  be 
written,  u*/G  ( A ,  B,  0,  *y),  where  A,  B,  0,  and  > 
are  similarity  functions.  This  relation  is  shown  in 
Figure  15.  These  functions  depend  on  stratifica¬ 
tion  and  the  single  similarity  parameter,  X,  which 
is  approximately  a  constant,  0.15.  Using  this  value 
and  stratification  information  one  can  analytically 
calculate  A,  B,  0,  and  >.  This  was  done  by  Brown 
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Fig.  14.  Sketch  of  the  velocity  profile  through  the 
boundary  layer  for  the  two-layer  model.  Velocity 
profiles  from  the  model  are  shown. 

(1982)  and  the  results  agree  with  traditional  curve 
Fits  to  several  data  sets.  The  result  is  that  the  sur¬ 
face  wind  and  stress  vectors  and  (less  accurately) 
heat  flux  are  determined  from  the  surface  pressure 
gradients,  which  determine  G  for  the  neutral, 
barotrophic  case.  One  needs  the  surface/air  tem¬ 
perature  difference  for  the  important  stratifica¬ 
tion  effect  on  both  layers.  The  horizontal  air  tem¬ 
perature  gradient  allows  calculation  of  the  thermal 
wind  effect.  The  neutral  roughness  parameter  and 
the  single  similarity  parameter  are  known  input 
values.  The  relative  humidity  adds  a  small  buoy¬ 
ancy  effect,  which  has  been  included.  The  model- 
predicted  1/,,/G  and  turning  angle  o  are  shown  in 
Figure  16. 


Fig.  15.  Summary  of  the  drag  relation  for  the  simi¬ 
larity  formulation  (from  Brown  1982).  Effects  includ¬ 
ed  are:  C — surface  layer  stratification:  (U,.  V,) — secondary 
circulation;  (U, ,  V,) — thermal  wind;  — variable  rough¬ 
ness. 
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Fig.  16.  Results  of  two-layer  model  for  windspeed  ratio,  U,o/G,  and  angle  of 
turning  <*  for  various  neutral  drag  parameters  (i.e.  various  roughnesses )  vs. 
stratification. 


The  MIZ 

The  basic  driving  force  for  the  PBL  flow  and 
the  fluxes  is  the  freestream  flow,  the  geostrophic 
or  gradient  wind  G.  This  must  be  determined  from 
the  balance  between  Coriolis  force  and  pressure 
gradient  force.  The  first  is  determined  by  the  lati¬ 
tude,  but  for  the  second  we  need  a  good  pressure 
field.  Data  on  surface  pressure  fields  are  routinely 
produced  by  the  national  weather  services.  How¬ 
ever,  particularly  in  the  high  latitudes,  the  data  are 


often  sparse  and  the  pressure  field  is  inaccurate. 
Figure  17  shows  such  a  predicted  field  for  AIDJEX 
compared  to  the  pressure  field  determined  by  the 
pressure  measurements  at  the  manned  stations  and 
a  ring  of  buoys  at  about  400  km  diameter.  While 
this  clearly  indicates  that  we  have  trouble  getting 
adequate  pressure  input  data,  other  correction 
factors  are  often  minimal  in  the  polar  regions. 
Over  the  ice,  roughness,  stratification,  and  the 
thermal  wind  (horizontal  temperature  gradients) 
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Fig.  18.  The  variation  in  the  surface  wind  for  a  constant  geos  tro¬ 
phic  wind  field  over  the  M1Z.  The  variation  in  drag  parameter, 
air/sea/ice  temperature  difference,  and  relative  humidity  produces 
variation  in  angle  of  turning,  speed,  ratio,  and  fluxes. 


are  often  fairly  uniform.  Thorndyke  and  Colony 
(1982)  have  shown  excellent  correlation  between 
simple  geostrophic  flow  from  surface  pressure 
Fields  and  ice  motion. 

However,  at  the  MIZ,  both  roughness  and  strat¬ 
ification  change  relatively  abruptly.  This  means 
careful  evaluation  of  the  PBL  How  field  assump¬ 
tions,  such  as  flow  equilibrium,  is  necessary.  The 
MIZ  produces  excellent  conditions  for  studies  of 
PBL  modification  on  the  mesoscale.  We  have  run 
our  /VxD  equilibrium  model  for  various  condi¬ 
tions  across  an  ice  edge.  The  variable  inputs  are 
the  surface  roughness,  directly  related  to  the  Cp 
for  neutral  stratification,  the  stratification  repre¬ 
sented  in  the  air/surface  temperature  difference, 
and  the  surface  relative  humidity.  The  outputs  are 
the  angle  of  turning,  the  wind  ratio,  t/l0/C,  and 
the  latent  and  sensible  heat  fluxes.  Figure  18 
shows  representative  values  across  an  ice  edge  for 
off-ice  winds.  The  turning  and  acceleration  of  the 
wind  are  evident.  A  large-scale  windfield  was  cal¬ 
culated  for  the  North  Sea  Experiment  (NORSEX) 
period,  shown  in  Figure  19.  Here,  only  the  change 
from  moderately  rough  ice  to  moderately  convec¬ 


tive  marine  conditions  is  calculated  on  a  120-km 
grid.  A  hypothetical  ice  edge  with  a  wavelength  of 
about  200  km  was  inserted  into  a  10-m  s'1  uniform 
geostrophic  flow  field.  The  magnitude  of  the  vari¬ 
ation  in  the  flow  at  the  MIZ  depends  on  the  width 
of  the  MIZ  and  the  rapidity  of  adjustment  to  the 
new  surface  conditions.  However,  it  is  instructive 
to  look  at  the  predicted  divergence  fields  for  this 
particular  case  (Figure  20).  The  imposed  ice  edge 
wavelength  is  evident.  The  interaction  with  the 
flow  field,  which  is  only  slightly  curved,  produces 
indications  of  large  divergence  values  in  the  MIZ 
vicinity. 

Finally,  to  get  down  to  the  scale  of  the  MIZ,  I 
will  speculate  a  bit  on  the  effects  of  non-equilibri¬ 
um.  Figure  21  shows  a  cross  section  with  transi¬ 
tion  zones  as  suggested  by  the  Bering  Sea  MIZEX. 
The  stratification  and  roughness  changes  from 
ocean  to  ice  are  as  shown.  The  Con  corresponds  to 
the  z0  roughness  parameter.  The  model  output  for 
the  equilibrium  flow  over  the  pack  ice,  the  MIZ, 
and  the  ocean  has  been  extrapolated  across  the 
edge  zones  by  imposing  a  linear  variation  in  Cq  to 
show  the  expected  variation  in  turning  angle  and 
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Fig.  19.  The  mesoscale  wind  field  calculated  during  NORSEX.  The  ice  edge  has  a  repre¬ 
sentative  wave  variation. 


t/,0/G,  depending  on  the  direction  of  flow.  The  Once  the  smoother,  convective  conditions  over  the 

bottom  graph  shows  the  resulting  effect  on  surface  water  are  encountered,  the  flow  speeds  increase 

stress.  For  instance,  on-ice  flow  has  higher  than  towards  the  ratio  t/.o/G  =  0.75.  However,  there  is 

equilibrium  surface  winds  when  it  first  hits  the  ice.  a  region  where  the  stress  is  lower,  reflecting  the 

This  means  that  U, » is  higher  than  is  represented  in  lower-than-equilibrium  windspeeds.  The  ramifica- 

the  Cp  where  the  wind  was  carefully  in  equilibri-  tions  for  ice  band  creation  and  movement  will  be 

um  for  parameterizing  measurements.  This  higher  left  to  the  legions  of  band  afficionados. 

momentum  will  produce  a  compact  edge.  When  When  the  airflow  is  parallel  to  the  ice  edge, 

the  flow  is  off-ice,  the  gradual  buildup  to  the  there  will  be  differential  drag  and  wind  turning, 

rough  ice  may  allow  equilibrium  values  of  surface  The  result  is  shear  and  convergence  or  divergence 

winds  with  a  ratio  of  t/,o/G  =  0.64  to  develop.  at  the  PBL  along  the  MIZ. 


Fig.  21.  A  summary  of  the  modified  PBL  flow  across  the  ice  edge. 
The  model  (equilibrium  flow )  results  show  constant  values  in  each  of  the 
four  zones  (thin  lines).  The  actual  variations  due  to  linearly  represented 
inertial  effects  are  shown  in  solid  lines. 


Conclusion 

It  is  evident  that  the  polar  regions  have  certain 
advantages  to  the  theoretical  PBL  modeller.  The 
ability  to  model  stress,  and  perhaps  heat  flux,  on 
the  mesoscale  using  similarity  and  geostrophic 
forcing  has  greatly  increased.  The  drag  coefficient 
is  available  to  an  accuracy  sufficient  to  reflect 
changes  in  PBL  stratification  and  ice  or  oceanic 
roughness.  The  conditions  at  the  MIZ  present  a 
relatively  sharp  and  uncomplicated  transition  re¬ 
gion  for  the  study  of  PBL  adjustment  times  and 
processes. 


There  are  suggestions  that  the  divergence  (or 
convergence  for  on-ice  flow)  is  significant  along 
the  MIZ.  There  is  probably  an  ice/sea  breeze  cir¬ 
culation  component.  This  resulting  wind  stress  can 
cause  an  oceanic  circulation  at  the  MIZ.  For  the 
along-ice  flows,  the  surface  winds  are  larger  over 
the  ocean  than  over  the  ice  for  the  same  geo¬ 
strophic  flow,  yet  the  surface  stress  is  larger  over 
the  ice.  There  are  evidently  several  possible  sec¬ 
ondary  circulations  forced  by  this  transition  zone. 

From  a  practical  standpoint,  to  get  an  accurate 
evaluation  of  these  variations,  it  will  be  essential 
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to  get  an  accurate  determination  of  the  forcing 
flow.  This  will  require  good  pressure  measure¬ 
ments  on  as  small  a  scale  as  instrument  accuracy 
will  allow.  At  80°  latitude,  a  1-m  s"‘  accuracy  in  G 
corresponds  to  0.7  mb  accuracy  in  instruments 
separated  by  400  km.  On  the  MIZ  scale  of  40  km, 
0.07  mb  accuracies  are  needed.  The  available 
NWS  pressure  maps  probably  yield  no  better  than 
2  or  3  mb  accuracy  on  this  scale  and  show  no  vari¬ 
ation  with  wavelengths  less  than  about  600  mb. 
They  will  have  to  be  supplemented  to  provide 
mesoscale  gradient  wind  information.  Alternate 
prospects  for  accurate  wind  input  data  are  aircraft 
survey  flights  in  the  PBL  and/or  satellite  ocean 
scatterometer-determined  winds  together  with 
buoy  winds  over  the  ice. 

When  the  most  recent  drag  coefficient  data, 
which  include  variations  due  to  stratification,  hu¬ 
midity,  and  roughness,  are  used  in  conjunction 
with  a  similarly  responsive  PBL  modei  for  the 
wind,  details  of  the  MIZ  stress  and  flux  field  can 
be  calculated.  The  limiting  factor  to  these  rela¬ 
tively  simple  bulk  formula  methods  for  the  meso¬ 
scale  to  synoptic  fields  is  the  accuracy  of  the  input 
wind  fields. 
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SUMMARY 

A  seabreeze  mechanism  which  may  contribute 
to  the  formation  of  leads  or  polynyas  in  the  mar¬ 
ginal  ice  zone  is  investigated  using  a  steady-state 
coupled  air-ice-ocean  model.  The  model  contains 
three  parts:  a  thermally  forced  boundary  layer  air 
flow,  a  mechanically  forced  ice  drift,  and  a  re¬ 
duced  gravity  ocean  model.  The  three  components 
are  linked  through  the  surface  temperature  gra¬ 
dient  and  various  stresses,  i.e.  air  stress  (on  ice), 
ice  stress  (on  water),  and  water  stress  (on  ice). 

In  the  marginal  ice  zone,  a  thermally  generated 
surface  wind  blowing  from  ice  to  water  (“ice- 
breeze”)  with  some  deflection  due  to  the  earth’s 
rotation  is  generally  divergent.  Such  divergent 
local  wind  fields  can  make  the  ice  flow  either  di¬ 
vergent  or  convergent,  depending  on  the  proper¬ 
ties  of  the  ice  and  water  beneath  them.  For  thin  ice 
and  a  thick  surface  water  layer,  the  turning  angle 
(the  angle  between  the  directions  of  the  surface 
wind  and  ice  flow)  is  small  and  the  ice  flow  is  di¬ 
vergent.  In  contrast,  for  thick  ice  and  a  thin  sur¬ 
face  water  layer,  the  turning  angle  is  large  and  the 
ice  flow  is  convergent.  The  air-ice-ocean  model 
confirms  this  result  and  gives  criteria  for  the  for¬ 
mation  of  leads  or  polynyas. 


INTRODUCTION 

In  this  study,  a  steady-state  coupled  air-ice- 
ocean  model  is  employed  to  investigate  the  condi¬ 
tions  for  ice  divergence  in  the  marginal  ice  zone 
(MIZ)  which  may  break  the  ice  and  cause  the  for¬ 
mation  of  leads  or  polynyas.  Although  upwelling 
of  warm,  deep  water  and  convective  processes 
probably  produce  polynyas  by  supplying  heat 
from  the  ocean,  another  mesoscale  seabreeze 
mechanism  which  causes  ice-breaking  may. exist. 


A  possible  mechanism  for  the  mesoscale  air- 
ice-ocean  interaction  in  the  MIZ  is  depicted  in 
Figure  1.  The  low-level  air  flow  generated  by  a 
surface  temperature  gradient  moves  from  ice  to 
water  (“icebreeze”),  with  some  deflection  to  the 
right  (or  lefi)  in  the  non'n  (or  south)  polar  region. 
This  icebreeze  should  have  a  maximum  speed  near 
the  ice  edge  and  decrease  both  iceward  and  water- 
ward,  i.e.  it  is  divergent  over  the  ice  and  conver¬ 
gent  over  the  water.  The  ice  flow  driven  by  the  sur¬ 
face  wind  field  is  to  the  right  (or  left)  of  the  wind 
in  the  north  (or  south)  polar  region.  The  turning 
angle,  defined  as  the  angle  between  the  directions 
of  ice  flow  and  surface  air  flow,  varies  from  0°  to 
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Fig.  I.  Ice-air-ocean  coupled  system. 
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90°,  depending  on  the  properties  of  the  ice  and 
water. 

The  steady-state  air-ice-ocean  interaction 
model  depicted  in  subsequent  sections  is  used  to 
determine  the  main  criteria  for  ice  divergence/ 
convergence,  and  therefore  for  lead  or  polynya 
formation. 


where  6^  is  the  basic  air  potential  temperature  at 
the  surface  and  N  is  the  Brunt-Vaisala  frequency 
of  the  atmosphere. 

The  mathematical  manipulation  is  almost  the 
same  as  in  the  second  section  of  Chu  (1986),  ex¬ 
cept  that  the  steady  surface  temperature  gradient 
is  used  here: 


THERMALLY  FORCED 
BOUNDARY  LAYER  AIR  FLOW 

We  utilize  the  same  planetary  boundary  layer 
model  treated  by  Kuo  (1973)  and  Chu  (1985,  1986) 
to  simulate  a  thermally  forced  boundary  layer  air 
flow.  The  coordinate  system  is  chosen  such  that 
the  x-axis  is  in  the  cross-ice-edge  direction  and  the 
y-axis  parallels  the  ice  edge,  as  shown  in  Figure  2. 
The  x-coordinate’s  unit  length  is  twice  the  M1Z 
width  L  (i.e.  200  km),  and  that  of  i'ne  vertical  co¬ 
ordinate  is  6  =  (i'/ft)° 5,  where  i>  is  the  vertical  ed¬ 
dy  viscosity  and  1)  the  angular  velocity  of  the 
earth’s  rotation.  The  line  x  =  Vi  is  located  at  the 
ice  edge.  The  MIZ  covers  the  zone  (0  <  x  <  0.5,  y) 
where  the  y-axis  (x  =  0)  is  near  the  boundary  be¬ 
tween  the  MIZ  and  the  interior  ice  pack.  It  is  con¬ 
sidered  that  spatial  variations  in  the  MIZ  are  much 
larger  perpendicular  to  the  ice  edge  than  parallel 
to  it,  and  hence  derivatives  with  respect  to  y  are  as¬ 
sumed  to  be  zero  (Lepparanta  and  Hibler,  1986). 

The  potential  temperature  of  air  is  divided  into 
two  parts:  a  basic  state  0B.(z)anda  perturbation  0.'. 
The  basic  state  is  given  by 

0b*(z)  =  Obo  +  (N^)z  (1) 


Y 


Fig.  2.  The  model  MIZ  and  the  coordinate  system. 


#.'(*.0)  =  -DT0  cosirx  (2) 

where  2 DT0  is  the  mean  surface  temperature  dif¬ 
ference  across  the  MIZ.  The  coordinates  and  at¬ 
mospheric  variables  are  nondimensionalized  by 
setting 

(x.,z.)  =  (xL,  zb) 
s’  =  0.7080  =  (DTo/flHo)* 
l«;  >  v.  ,  w,-  j  =  oi u  •,  v •  w-  u/L.i 

pf  =  (ghDT0/em)PM 
where 

U  =  gbDTo/(2LQ0w)  (4) 

is  the  scale  of  the  icebreeze. 

The  surface  wind  driven  by  the  horizontal  tem¬ 
perature  gradient  is  computed  by 

n(a,(x)  =  u,a)  simrx  (5) 

vla)(x) 5  v,a)  sin  7rx  (6) 

where 

3 

U(a)  m  U  £  fljXj 

7=1 

3 

vfa)  m  U[b-fo  £  fljXj  P(Xj,  r)] 

7=1 

/o  =  sin<Z> 

r  =  xJ62i'(,/(2/.2i>) 

F(\,  r)  =  l/(X2/2-r) 

and  <t>  is  the  latitude  and  vh  the  horizontal  eddy  vis¬ 
cosity.  Eigenvalues  X>  {J  =  1,2,3)  are  the  roots 
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to 


with  negative  real  parts  of  the  following  sextic 
equation: 


X‘-4/V  +  Mfo 1  +  r5)X'-4xJ/?i  =  0  (8) 

where  the  Richardson  number  R{  is  defined  by 

Ri  =  6WV(4L'n!).  (9) 

The  coefficients  a,,  au  a,  and  b  in  eq  5  and  6  are 
determined  by  (Chu,  1986) 

L  «j  =  0 

j=  i 

3 

E  Xj(i-/fXj)flj  =  o 

;=i 

3 

-/o  ^Xjd-ATXj)  F(V)aj  +  />(l  =  0 

y  =  i 

3 

xRi  £  fXXj,r)oj  +  V2r /ob/ it  =  -1 

;=i 

where  K  is  a  measure  of  the  effective  depth  of  the 
constant  stress  sublayer  (Kuo,  1973): 

K  =  i»/(Ca6)  (N) 

where  Ca  is  a  dimensional  (m  s'1)  air  drag  coeffi¬ 
cient. 


FREE  ICE  DRIFT  MODEL 

The  linearized  steady-state  momentum  equa¬ 
tions  for  a  free  drift  ice  model  are 

/v(0  +  [T<«.0  +  ^<7]/(CiWj)  =0  (12) 

-fu(f)  +  +  ^'^/(Ci/Zi)  =  0  (13) 

where  /  =  ( =  2fl  sin<*>)  is  the  Coriolis  parameter, 
H\  the  mean  ice  thickness,  the  ice  density,  V®  = 
[m(0,  v*0j  the  ice  velocity,  fa,i)  =  rL<’,')]  the 

wind  stress  on  the  ice,  and  riKi>  =  l4w,°,  the 
water  stress  on  the  ice.  The  air  and  water  stresses 
on  ice  are 


>■»  =  c.C,  i>> 

>•9  =  6w Cw  F®| 


(14) 


where  =  [k<0),  v(a>]  is  the  surface  air  velocity, 
F(w)  =  [m(w),  v<w)]  the  water  velocity  and  Cw  the  di¬ 
mensional  (m  s'1)  water  drag  coefficient. 


REDUCED  GRAVITY  OCEAN  MODEL 
WITH  RAYLEIGH  DISSIPATION 

The  ocean  model  is  a  two-level  model  in  which 
the  lower  layer  is  assumed  to  be  deep  enough  that 
motion  within  it  is  vanishingly  small.  Such  a 
model  is  referred  to  as  a  reduced  gravity  model. 
The  momentum  equations  are 

/?w<“)-/v(“>-g//wAew/(ewz?)aju(w)/a^ 


=  r<r>/(ew//w)  (15) 

R  v(M,)  +  =  r<;»/(ew/Z w)  (16) 

where  R  is  the  coefficient  of  Rayleigh  friction,  //w 
the  mean  thickness  of  the  surface  water  layer,  pw 
the  density  of  the  surface  water,  Apw  thejdensity 
difference  between  the  two  layers,  and  r',w>  = 
[r(x'>),  r(y- W)1  the  ice  stress  on  the  water.  has  the 
same  magnitude  but  the  opposite  direction  as 
i.e. 

>">  =  ->.0  =  -ewCw[F<w)-F<0].  (17) 


ICE  DIVERGENCE/CONVERGENCE 
CRITERIA 

The  surface  wind  is  the  forcing  term  of  the  ice- 
ocean  system.  Since  our  model  is  linear,  the  me¬ 
chanically  forced  ice  and  water  motion  V®  and 
V should  have  the  same  Fourier  component  as 
the  forcing  term  Ma),  i.e. 

uu\x)  =  u(i)  sinirx,  vw(x)  =  vro  sinxjc  (1.8) 

u<w)(x)  =  u(w)  sin  rx,  v(w\x)  =  v,w)  sinirx  (19) 

where  x  is  a  nondimensional  coordinate. 

Substituting  eq  5,  6,  14,  17,  18  and  19  intc  eq 
12, 13, 15  and  16  we  have  the  following  four  linear 
algebraic  equations  for  the  Fourier  coefficients 
«(*),  v(w),  mw  and  v(<): 
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(/?  +  ir2g//w  Sqw/(QwRL})  +  u',]u'w) 


-Zv'“,-ce)  u0)  =  0 

(20) 

Zu'“’)  +  (/?  +  W»)v'“,-a)l  v">  =  0 

(21) 

-u)j  U,W)  +  (J> 2  u(i)-f  v(i)  =  w,  u,a) 

(22) 

-U)2  v'*>  +  Z U(,)  +  U0)  =  0),  V(a) 

(23) 

Since  our  model  is  ^-independent,  the  ice  diver¬ 
gence/convergence  criteria  are  simply 

>  Ice  divergence 

u‘*  =0  No  ice  divergence/convergence  (30) 

<  Ice  convergence 

Using  eq  25,  the  criteria  (eq  30)  are  rewritten  as 


where  the  parameters  w,,  and  w,  are  defined  by 
s  eaCa/(ei//i) 


u,  =  ewcw/(eiHj) 


(24) 


O)}  —  Cw/ 


Solving  the  nonhomogeneuus  linear  algebraic 
equations  20-23  we  obtain  the  amplitudes  of  the 
ice  and  water  flow  fields: 

U,‘>  =  uJ,[a'2(l-«2)Mro,+/(l  +  a,u):)vfa>]/ 

[o)2(  1  -«, )( 1  -«j)  +/( 1  +  ajo>2)J]  (25) 
v">  =  «,|tf,(l-(rl)vw-/(l  +  a,o)2)u,a)]/ 

Md-a.Ml-aO+Zd+a,^)2]  (26) 
u(w)  =  ^[u^R  +  w^+fv1^/ 


[«„(/?  + CO,) +/2] 

v^»  =  u),[a„Vrt,-/wr',]/[ao(/?  +  a),)+/2] 

where 

a0  =  R  +  a)j  +  5r2gf/wAgw/(ew/?Z.2)  • 

Of  I  —  U)j(/?  +  O)j)/(c#o (/?  +  U!])  +  Z2] 

a;  =  W)a0/[ao(/?  +  ajj)  +  Z2] 
a,  3  aJ)/[«0(A’ +  ajj) +Z2]- 


(27) 

(28) 


(29) 


>  Ice  divergence 

Cr  =  |  v,a^/utt)  |  No  ice  divergence/convergence  (3 1 ) 
<  Ice  convergence 

where 

Cr  3  a>,d -«,)/( IZId +«>«!>] 

and  ;3  =  tan'1  |  v'a)/ufa)  j  is  the  icebreeze  deflection 
angle  due  to  the  earth’s  rotation.  In  the  present 
calculation  wc  fixed  aii  values  of  the  parameters 
except  for  (1)  the  Brunt-Vaisala  frequency  of  the 
atmosphere  N,  (2)  the  mean  ice  thickness  H„  and 
(3)  the  mean  depth  of  the  water  surface  layer  //w. 
Table  1  lists  the  values  of  the  fixed  parameters. 
Figure  3  describes  the  dependence  of  |  v<a)/u,a)  |  on 
R\  at  three  different  latitudes  (65°,  70°  and  75°). 
\vM/ufa)\  increases  with  increasing  latitude  and 
decreases  with  increasing  Rt.  At  65°  (N  and  S)  lati¬ 
tude  it  varies  from  0.9  (d  =  48°)  when  R,  = 
O.lOx  10~2  to  0.53  (/3  =  27.5°)  when  R{  = 
4.03  x  10‘2. 

We  also  computed  the  parameter  Cr  as  a  func¬ 
tion  of  Hj  and  //w  at  three  different  latitudes  (65°, 
70°  and  75°),  where  //j  varies  from  0  to  10  m  and 
Hw  from  0  to  50  m.  The  distribution  of  Cr  with  H\ 
and  Hw  is  shown  in  Figure  4.  High  values  of  Cr  are 
located  in  the  region  of  large  Hw  and  small  H\. 
The  smaller  the  H,  (or  the  larger  the  //w),  the  larg¬ 
er  the  Cr.  It  is  quite  understandable  that  the  water 
current  (*w)  is  very  weak  as  Hv  becomes  large,  be¬ 
cause  the  mechanical  forcing  terms  in  eq  20  and 
21,  u),u(l)  and  u),v(,),  become  small.  When  the  water 
current  F,w)  is  weak  (refer  to  large  //w)  the  turning 
angle  largely  depends  on  the  ice  thickness.  For 


Table  1.  The  standard  model  parameters. 


L  =  200  km 

v  =  5  m2  s'1 

Vh  =  10’ m2  s'1 

fi  =  0.7292x10-  s- 

g  =  9.81  m  s'2 

R  =  2.08x10-’  s  ' 

»bo  =  270  K 

Cw  =  1.76x  10"’  m  s'1 

Ca  =  3x  10‘2  m  s’1 

<?,  =  1.29  kg  m  ’ 

ew  =  10J  kg  nr’ 

e,  =  910  kg  nr* 

=  0.5  kg  m  ’ 

DT<,  =  5  K 

<t>  =  (65°,  70°,  75°) 
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Fig.  4  (cont’d).  Isolines  of  Cr  in  the  (H „ 
plane  for  three  different  latitudes. 


Hy(B> 

Fig.  5.  H,-Hw  curves  for  zero  ice  divergence/convergence  for  five  different  values  of  the  at¬ 
mospheric  Richardson  number:  (a)  R,  =  4  x  /O'1,  (b)  R ,•  =  3xl0~\  (c)  R,-  =  2xlO~2,  (d)  R,- 
=  10  *  and  (e)  R,  =  10  '. 


thin  ice,  eq  22  and^  23  show  that  deflects  by  a 
small  angle  from  therefore  the  ice  flow  is  di¬ 
vergent. 

Solving  the  algebraic  equation  23  for  H,  after 
setting  =  0,  we  get  the  dependence  of  Hx  on 
Hw  for  the  no  divergence  or  convergence  case,  as 
shown  in  Figure  5.  These  curves  separate  the  pa¬ 
rameter  plane  (H\,  H„)  into  ice  convergence  and 
ice  divergence  parts.  Figures  3-5  help  us  to  justify 
the  ice  divergence/convergence  in  the  MIZ,  and  to 
predict  the  formation  of  a  lead  or  polynya. 


AIR,  ICE  AND  WATER  FLOW  FIELDS 

We  compute  F(a),  and  V<w)  for  the  north  and 
south  polar  regions.  The  atmospheric  Brunt-Vai- 
sala  frequency  N  is  taken  as  10"1  s'1  (correspond¬ 
ing  to  a  Richardson  number  Rx  of  0.81  x  10'*).  The 
parameters  //w  and  H-t,  which  can  be  indicated  as 
points  in  the  plane  (Figure  5),  are  chosen  to 
represent  three  cases:  (1)  thin  ice-thick  surface 
water  layer  (a,  and  0,),  (2)  thin  ice-thin  surface 
water  layer  (a,  and  0,),  and  (3)  thick  ice-thin  sur- 
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face  water  layer  («,  and  *j>)-  The  points  <*,  (40  m, 
1.5  m),  ft.  (10  m,  1.5  m)  and  (10  m,  6  m)  are  for 
the  north  polar  region,  and  points  (40  m,  0.7 
m),  d:  (10  m,  0.7  m)  and  (10  m,  6  m)  are  for  the 
south  polar  region.  In  Figure  5,  points  <*,  and  ,i, 
are  located  in  the  ice  divergence  area,  <*.  and  in 
the  no  divergence  or  convergence  area,  and  <»,  and 
>i,  in  the  ice  convergence  area. 

Figures  6  and  7  show  F*a),  F*'1  and  F"*1  in  the 
north  and  south  polar  regions  (65°),  respectively. 


They  indicate  the  following  results: 

(1)  The  surface  wind  field  does  not  change  with 
«w  and  H\.  When  the  temperature  difference 
across  the  MIZ  is  I0°C  (i.e.  DT0  =  5°C),/V=  10  J 
s'1  and  <j>  =  65°,  the  maximum  surface  wind  speed 
can  reach  17  m  s  '. 

(2)  The  surface  wind  field  varies  from  north  to 
south  polar  regions.  In  the  north  (south)  polar  re¬ 
gion,  the  surface  wind  blows  from  the  ice  to  the 
water  with  a  deflection  to  the  right  (left). 
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a.  H,  =  1.5  m,  Hw  =  40  m. 


b.  H,  =  1.5  m,  Hw.  =  10  m. 


Fig.  6.  V3>,  Vu  and  V*'  in  the  north  polar  region  (<t>  =  65  °N,  N  =  10'1  s  '). 
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Fig.  6  (cont'd).  Va\  and  Yw<  in  the  north 
polar  region  (<t>  =  65  °N,  N  =  10~2  s~'). 


(3)  The  ice  flow  varies  with  //w  and  Hv  For  a 
thin  ice-thick  surface  water  layer  (points  a,  and  0, 
in  Figure  5),  the  ice  flow  reveals  evident  diver¬ 
gence,  as  shown  in  Figures  6a  and  7a.  The  maxi¬ 
mum  ice  flow  velocity  is  at  the  ice  edge.  It  reaches 
40  cm  s'1  (40.7  cm  s  ')  for  =  1.5  m,  =  40 
m,  and  <t>  =  65 °N  (H\  =  0.7  m,  //w  =  40  m,  <t>  = 
65  °S).  Such  steady-state  ice  divergence  will  break 
the  ice  and  cause  the  formation  of  a  lead  or  polyn- 
ya.  For  the  thick  ice-thin  surface  water  layer  case 


(points  a,  and  0,  in  Figure  5),  the  ice  flow  exhibits 
obvious  convergence  (Figures  6c  and  7c).  The 
maximum  ice  flow  velocity  is  also  at  the  ice  edge. 
It  reaches  41 .4  cm  s'1  for  H\  =  6  m,  //w  =  10  m, 
and  <t>  =  65  °N  (65  °S).  Such  steady-state  ice  con¬ 
vergence  will  prevent  the  ice  from  breaking.  For 
the  thin  ice-thin  surface  water  layer  case  (points  «2 
and  02  in  Figure  5),  the  ice  flow  displays  no  diver¬ 
gence/convergence  (Figures  6b  and  7b). 
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